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ISOELECTRIC SPECTRA AND MOBILITY SPECTRA: A NEW 
APPROACH TO ELECTROPHORETIC SEPARATION* 
By ALEXANDER KOLIN 
UNIVERSITY OF CHICAGO 
Communicated by R. S. Mulliken, January 31, 1955 


Introduction. Fluids which are of interest to the biochemist and microbiolo- 
gist are often complicated mixtures containing some components in solution and 
others in suspension. The separation and classification of these components con- 
stitute a problem of basic importance. Physical methods, ultracentrifugation and 
electrophoresis, are the most widely used general methods for accomplishment of 
these aims. The moving-boundary method employed in electrophoresis has the 
drawback that only the slowest and fastest ions can be isolated from the other 
components. Only a small fraction of the separable components can be isolated. 
The amounts of fluid required are relatively large, and the time required for analy- 
sis may be as long as 2 days.' The observation of the moving boundaries requires 
rather elaborate optical equipment, and special precautions must be taken to 
eliminate thermal convection. Filter-paper electrophoresis obviates some of the 
disadvantages of the moving-boundary method. The problem of thermal con- 
vection is eliminated, and the components of a mixture can, in principle, be sepa- 
rated from each other. Certain new drawbacks are introduced, however. For 
instance, there are proteins which are adsorbed by the filter-paper fibers and thus 
do not migrate and cannot be separated. The detection and quantitative deter- 
mination of colorless components is cumbersome, requiring staining. This is a 
disadvantage in preparative work. 

The time required in the analysis of a protein mixture by the method of filter- 
paper electrophoresis is approximately 12 hours.' For comparison, the method 
described below accomplishes a complete separation of a mixture, as a rule, in less 
than 10 minutes. In the standard methods the different components are character- 
ized by their electrical mobilities. The new method offers the prospect of making 
it possible to identify the individual protein components by their isoelectric points. 

The purpose of this paper is the presentation of a simple electrokinetic method of 
analysis which is applicable to suspensions as well as to ionized solutions. The 
components of a given mixture are separated from each other simultaneously and 
are sorted in a spatial arrangement to which we shall refer as an “electrophoretic 
“spectrum” is used here to designate a sorting in the 
” “energy spectra,”’ and “fre- 


line spectrum.” The term 
same sense in which one speaks of “mass spectra, 
quency spectra.” 
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We shall consider two kinds of electrophoretic line spectra: (1) ‘mobility spee- 
tra’ and (2) “isoelectric spectra.”” In ,the spectrum of the first kind the com- 
ponents are separated according to their differences in electrophoretic mobility, 
whereas in the spectrum of the second kind the components are separated accord- 
ing to differences in their isoelectric points. Both types of spectra may be super- 
imposed upon each other; in this case we shall speak quite generally of an “electro- 
phoretic line spectrum.” 

Mobility Spectra.—The aim of the following procedure is to condense the ions 
of the mixture into a thin layer (often » small fraction of a millimeter in thick- 
ness) whose ionized components, migrating with different velocities, finally separate 
into several layers. The side view of these layers offers the aspect of thin lines 
reminiscent of an optical line spectrum. Since, to achieve separation, two inter- 
penetrating layers do not have to move apart more than the small distance equal 
to the thickness of the layer, great speeds of separation can be achieved. 

Figure | shows the electrophoresis cell in which the electrophoretic line spectra 
are produced. The central U-tube section communicates at the top with the elec- 
trode compartments A and B. The bottom of the U-tube is filled with a buffer 
solution nearly saturated with sucrose to increase its density. The upper portions 
of the U-tube, as well as the electrode compartments A and B, contain a sucrose- 
free buffer solution. The layer VW (whose density is adjusted by addition of sucrose 
or glycerol so as to stabilize it between the light and heavy buffers) contains the mix- 
ture tobe analyzed. ‘The mixture is very dilute (about 0.01—0.1 per cent of dissolved 
or suspended solids), and its electrical conductivity is very small as compared to 
the conductivity of the buffer solution. 

Let us consider one ionic component migrating in the VM layer under the influence 
of the electric field. Due to the very low conductivity of this layer, the potential 
gradient in VW is very much larger than in the buffer, and hence the ions migrate in 
this region very much faster than they do in the buffer. Suppose that negative 
ions are moving from boundary / toward boundary 2 (Fig. 1, lower section) with a 
speed v,, uj /o» (Where j is the current density, u the ion mobility, and o,, the 
electrical conductivity of the M layer). The ion velocity will suffer a change when 
the ions enter the buffer, where v, uj/o,. When the last ions have reached 
boundary 2, the front of the ions has traversed in the buffer the distance d = rf, 
whereas during the same time the distance M covered in the M layer is M = vy. 
Hence 
T m 


M 


Tr 


Since om, <o,: d<M. Thus the ions, which originally fill the wide belt VM, are 
contracted into a thin layer of width d as they enter the buffer. If the zone VM con- 
tains a mixture of ions and of charged, suspended particles, practically one single 
common line is formed at the boundary of Mf. Ions of opposite sign are condensed 


at the opposite boundaries of M. 

After the formation of the common line, the current is continued until the mi- 
grating line splits into as many lines as there are separable components in the mix- 
ture, 

Figure 2, A, B, C, illustrates a mobility spectrum obtained with a mixture of three 
dyes, Figure 2, A, shows the VM layer before passage of current. (The marks / 
and 2 of Fig. | [lower section], indicate the positions of the boundaries of the M 
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Pig. |. -L, and Le: sources of light; £; and Ey: electrodes; A and B: channels filled with 
buffer: L and R: legs of the U-tube: 1/ mixture to be analyzed; m: boundary between light 
and heavy acid buffer; J and 2: marks indicating the positions of the boundaries of the mixture 
layer before separation, Upper section top view ower section front view. The dimensions 


are in mm 
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laver in all photographs.) Figure 2, B, shows the “‘common lines” formed at the 
M boundaries 2 minutes after the start, using a current of 13 ma at a voltage of 110 
volts across the electrodes. (Approximately the same voltage and current were 
used in all the experiments described below.) Figure 2, C, shows the splitting of 
the lower line, 12 minutes after starting the current, into two components which 
are easily identified by their color. Figure 2, D, shows, for comparison, the sepa- 
ration of two species of suspended green algae by the same method. The lines are 
much thinner and sharper. ©The time is 14 minutes. The same buffer of pH 7.3 


is used in all compartments (A, B, L, PR). 





Fic. 2.—A, B, ¢ mobility spectrum’’ of a mixture of three dves: methylene blue (a@), 
chrome hematoxylin (8), and bronphenol blue A> The M layer before separation. B: 
lons of opposite signs are condensed into “lines’’ at the 7 boundaries after 2 minutes. CC:  Sep- 
aration of 3 and 4 after 12 minutes. pH range, 7.7 —~ 7.7 The lines are broadened by diffusion. 

D; “mobility spectrum’’ of a mixture of two species of green algae, Ankistrodesmus (a) and 


Chlorella (8), after 14 minutes. pH range, 7 


” 
o> fed. 


It is possible to introduce a relatively thin layer (about 1 mm. thick) of a buffered 
mjxture (whose density has been appropriately adjusted by addition of glycerol) 
between two buffer columns of equal pH but of densities lying below and above 
the density of the mixture laver. This arrangement is a three-dimensional ana- 
logue to filter-paper electrophoresis. This method is not satisfactory, since the 
traveling ion layers are quickly erased by thermal convection. The performance 
ean be improved by creating a density gradient in the lower buffer column in which 
the ions migrate. The density gradient can be created by stirring the boundary 
between two buffer columns containing glycerol or sucrose in different concentra- 
tions or by superimposing many short columns whose densities vary in steps. 
The “ion condensation method” described above, however, is much more satisfac- 


tory, since the ion lavers can be made exceedingly thin and sharp, so that only 
short distances have to be traversed for complete separation. Thus the traveling 


ion lavers remain throughout the process sufficiently close to the steep density 
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gradient at the ./ boundary, where thermal convection is most effectively sup- 
pressed, 

Isoelectric Spectra-—In the separation of ampholytes, such as proteins, poly- 
peptides, and amino acids, special use can be made of the existence of an isoelectric 
point, in order to obtain particularly sharp and stable “line spectra.” The dis- 
tance which the ions must traverse in the process of separation can be limited by 
interposing the mixture between two buffers of suitable pH. The separation time 
can thus be greatly reduced. Thus far, the method has been applied to proteins 
only. 

The principle of the method is based on the following consideration: Let us 
imagine a column of electrolyte solution in which a pH gradient has been estab- 
lished. We shall assume that this solution also contains a protein at uniform 
concentration. The charge on the protein ions, being a function of the pH of ‘the 
solution, will vary throughout the pH gradient. For every ampholyte there is a 
characteristic pH at which the net charge of the molecule is zero. As we proceed 
from this point toward increasing pH values, the protein molecules acquire a 
negative charge, which increases with increasing pH. On the other hand, as we 
proceed toward the acid side of the gradient (decreasing pH), the protein molecules 
acquire a positive charge, which increases as the pH diminishes. An electric field 
of constant intensity will not move the protein ions in a constant direction with 
constant speed in our liquid column of varying pH. If the current flows in the di- 
rection of increasing pH, the protein ions which are on the acid side of the isoelectric 
point move with the current toward the isoelectric point, being positively charged. 
The protein ions on the basic side of the isoelectric point are negatively charged 
and hence move contrary to the direction of the current, i.e., likewise toward the 
isoelectric point. Thus the protein ions converge toward the isoelectric point, where 
all the protein can be eventually concentrated. Since different proteins have, as a 
rule, different isoelectric points, the idea suggests itself of separating different com- 
ponents of a protein mixture by concentrating them at their characteristic iso- 
electric points.’ We shall refer to a sorting pattern obtained on the basis of this 
principle as an ‘‘ isoelectric spectrum.”’ 

teversal of the current should produce the contrary effect, namely, a dilution 
of the protein solution at the isoelectric point due to migration of the protein ions 
away from this point. The predicted dilution effect, as well as the concentration 
and separation effects, are illustrated in Figures 3 and 4. 

The technique of producing the pH gradient for the isoelectric spectra, shown 
in the preceding photographs, is very similar to the method of obtaining a con- 
ductivity gradient for mobility spectra shown in Figure 1. The protein mixture 
of low electrolyte content is dissolved in a sucrose or glycerol solution and intro- 
duced at MV into the U-tube (Fig. 1, lower section). The dense sucrose buffer solu- 
tion filling the lower portion of the U-tube is an acid buffer of appropriate pH nearly 
saturated with sucrose. The buffer filling the leg R above mark m, as well as the 
chamber A, is the same acid buffer free of sucrose. A basic buffer of appropriate 
pH} fills the upper half of leg L and the chamber B. This arrangement establishes 
a pH gradient in the laver / by convection and diffusion. The pH varies from pH, 
at mark 7 to pH. at mark 2. To yield an isoelectric spectrum, the isoelectric 


points of the proteins must lie between these limits; otherwise, a mobility spec- 


trum of the proteins is obtained. 
Ideally, one would desire the pH distribution in MW not to be affected by the cur- 
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rent. WKohlrausch® and Dole® have shown for the simple case of a strong electrolyte 
column containing two species of ions whose concentration varies arbitrarily along 
the axis that the concentration distribution is not affected by the current. They 
also showed for strong electrolytes that the current does not affect the concentra- 
tion distribution in a solution containing many species of ions, provided that the 
concentration ratio of the ions does not vary along the current axis. A concentra- 
tion distribution satisfying this condition can be obtained by dilution of the elee- 
trolyte. These results are, however, not applicable with a sufficient degree of ap- 


proximation to solutions of weak electrolytes such as are encountered in buffer 


Fig. 3.--Eleetrophoretie concentration and dilution of hemo- 
globin in a pH gradtént. pH range (2.69.6). A: hemoglobin 
lavers before passage of current (the polarity is as indicated in 
Fig. 1) B concentration effect in left leg and dilution effect in 
right leg after 40 seconds. 

systems.’ In the absence of a satisfactory theory from which the required butter 
compositions for a stationary or nearly stationary pH distribution could be com- 
puted, a suitable combination of buffers was found empirically.* 

The “isoelectric lines” of proteins are much sharper and more stable than 


mo- 
bility lines” of proteins and of nonproteins. The greater sharpness of the isoelectric 
lines is due to the stable equilibrium of protein molecules in the isoelectric plane. 


Should a protein molecule be deviated from this position by such effects as diffusion 


or convection, it will acquire a charge of such sign as to be returned to its original 
position by the electric feld. As a rule, an isoelectric spectrum is formed more 
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rapidly than a mobility spectrum. It is confined mostly to the .W zone but may 
project somewhat beyond its boundaries. A mobility spectrum, on the other 
hand, lies entirely beyond the boundaries of 7. Figure 5 shows an isoelectric spec- 
trum of proteins. Four proteins were separated between buffers of pH 4.8 and 7.7. 
The proteins were cytochrome C, hemoglobin, catalase, and dyed collagen (‘ Azo- 
coll’). The cytochrome C line is not part of the isoelectric spectrum. The iso- 
electric point of eytochrome C lies at pH 10.6, outside the range of the pH gradient; 
thus eytochrome C forms a line at the upper 7 boundary due to retardation in the 
buffer. The isoelectric points of hemoglobin (pH 6.8), catalase (5.7), and collagen 
(5.3) he within the range of the pH gradient, so that the three lower lines form an 
isoelectric spectrum between the J/ boundaries. The separation shown in Figure 5 


y 


Was accomplished in 4 minutes. 


Ql 


oe 


ie ee 





Fia. 4.-Separation of eyvtochrome Fic. 5.—Electrophoretic line spee- 

C(a@) and hemoglobin (8) after 2'/« trum obtained in 4 minutes with four 
ida amg gO —e- ler : 

minutes, pH range, 3.3 — 7.7. The proteins: cytochrome C (a), hemo- 


hemoglobin line is. produced by “iso- globin (8), catalase (4), and dyed col- 

electric condensation, whereas the lagen (“Azocoll’’) (6). pH range, 

cytochrome C line is produced by re- 187.7. ais formed by retardation 

turdation at the buffer boundary at the buffer boundary, whereas 8, 

and 6 are formed by isoelectric 

condensation 
Electrophoretic Line Spectra. Since the arrangement for production of isoelec- 
tric spectra combines a conductivity gradient, such as is used for mobility spectra, 
with a pH gradient, it is evident that the above procedure will be effective in pro- 
ducing simultaneously a mobility spectrum as well as an isoelectric spectrum. The 
mixture may thus contain ampholytes as well as nonampholytes, and the iso- 
electric points of some of the proteins may lie outside the range of the pH gradient. 
This is the most general arrangement for production of “electrophoretic line 
spectra.”” In Figure 5 the eytochrome C line, for example, is not formed by tso- 
eleetric condensation, 

Whenever possible, it is advantageous to accomplish the separation by producing 
an isoelectric spectrum rather than a mobility spectrum. The lines are thinner, 





108 BIOCHEMISTRY: ALEXANDER KOLIN Proc, N. ALS. 


sharper, and more stable. The isoelectric spectrum is formed more rapidly, since 
the ions move in the .V laver, where the potential gradient is much higher than in 
the buffers in which the mobility spectrum is formed. A wider separation between 
the lines can be obtained in a short time. The distance traversed by the ions and 
hence the separation time can be made very short by making the height of the 
M layer small. The sharp density gradient at the adjacent 1 boundaries, com- 
bined with the gentle density gradient in the VW layer, accounts for the effective 
suppression of disturbances due to thermal convection in the W layer, where the 
isoelectric spectrum is formed. ‘The final prospective advantage of the isoelectric 
spectrum lies in the possibility that the pH could be measured by means of micro- 
electrodes at the sites of the individual isoelectric lines. This would permit one 
to identify the ampholytes by a rapid determination of their isoelectric points. 

In the preceding photographs colored substances have been used the lines of 
which can be seen in transmitted light provided by the source L, of Figure 1 (upper 
section). Lines due to colorless proteins can be made visible by means of dark- 
field illumination. A source of a nearly parallel beam of light (1, of Fig. 1) illumi- 
nates the laver \/ from the left side in reference to the observer. The Tyndall 
effect is usually quite strong, so that there is no difficulty in photographing the “line 
spectrum.” Figure 6 shows the colorless components of a commercial catalase 
preparation (Armour and Company). Figure 6, A, shows the dark catalase line 
in transmitted light, and Figure 6, B, shows the two colorless components on the 
acid side of the dark line as seen with dark-field illumination. 

Suspended particles may be separated with the same ease as protein molecules. 
This may be used for separation of microérganisms, various body cells, cell debris, 
ete. Figure 7 illustrates the separation of red blood cells from the two species of 
green algae Chlorella and Ankistrodesmus suspended in a 20 per cent sucrose solu- 
tion. 

Figure 8 shows the separation of tobacco mosaic virus from a preparation of Té6r-4 
virus. Figure 8, A, shows a single line obtained with the tobacco mosaic virus 
demonstrating the purity of the preparation. Figure 8, B, shows two lines ob- 
tained with the preparation of T6r+ virus. In Figure 8, C, one sees the three com- 
ponents into which a mixture of the T6r + virus preparation with the tobacco 
mosaic virus is resolved. The separated fractions of a mixture may be extracted by 
means of a micropipette 

The present paper is a preliminary communication. A detailed theoretical and 
experimental analysis of the method will be presented in a later publication. 


The author wishes to thank his colleagues Professor Martin E. Hanke and Dr. 


Nicholas Nicolaides for numerous biochemical preparations as well as for stimu- 
lating chemical discussions. The colored proteins provided by Dr. Nicholaides 
proved especially helpful. Illuminating theoretical discussions with Dr. Daniel 
Leenoy are much appreciated. The viruses have been kindly furnished by Mr. L. 
Barrington and the algae by Mrs. H. Gaffron and Drs. F. L. Allen, J. E. Brugger, 
and E. Kessler. Thanks are also due Mr. John Peters, of the chemistry shop, and 
Mr. Edward Wolowiec for their able technical assistance. 

The support of this work by generous grants from the Office of Naval Research, 
the Research Corporation, and the Smith, Kline, and French Foundation is grate- 
fully acknowledged. 
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Kia. 6 electrophoretic line spectrum of a commercial catalase prep- 
Company a: the dark fine photographed by 


aration (Armour ind 
transmitted light 38and >: color 
field illumination Separation time 


less components photographed by dark- 
O minutes. pH range, 4.8 — 7.7 





Fic. 7.—Electrophoretic line spec- 
trum of a mixture of the algae Ankis- 
trodesmus (a) and Chloreila (8), with 
human red blood cells (7). Separation 
time, 9 minutes. pH range, 6.2 — 7.7. 
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Fic. 8.-Electrophoretic line spectrum of a virus mixture (dark-field illumination pH range, 


$8— 7.7. Separation time, 8 minutes. A single line (a) obtained with tobacco mosaic virus, 
B: the two components (7; and rz) of a T6r+ virus preparation. C; three components (7 
p, 72) obtained with a mixture of tobacco mosaic virus with the T2r+ virus preparation. 


* These studies were aided by a contract between the Office of Naval Research, Department 
of the Navy, and the University of Chicago 

1H. J. Antweiler, Die quantitative Elektropho ese in der Medizin (Berlin: J Springer, 1952), 
p. 38. 

2A. Kolin, J. Chem. Phis., 23, 407, 1955. 

\. Kolin, thid., 22, 1628-1629, 1954 

‘ Buffers of different pH were prepared from two stock solutions which were mixed in varving 
proportions to obtain buffers of a desired pH. The composition of the basie buffer of pH 9.6 
was 19.43 gm. of sodium acetate (-+3H.O) and 20.43 gm. of sodium barbital per liter of solution 
The acid buffer of pH 2.6 was obtained from the basie buffer by adding | liter of the latter to 3 
liters of 0.1 V HCl and 0.6 liters of water. 

The pH values of buffer mixtures used in various experiments are indicated in the legends to 
the figures. For example, the statement “pHi range, 4.8 — 7.7” means that the pH of the acid 
buffer below the WV laver is 4.8 and of the basic buffer above the VW laver, 7.7 

> F, Kohlrausch, dann. Physik, N.F., 62, 209, 1897. 

6 V_P. Dole, J. Am. Chem. Soe., 67, 1119, 1945. 

7H. Svensson, Acta Chem. Scand., 2, 841, 1948 

§ The stability of the pH gradient could be improved appreciably by adding an appropriate 
amount of NaCl to the basic buffer (following a suggestion by D. Leenov). For the pH 7.7 buffer 
used in these experiments, | cc. of a5 N’ NaCl solution was added to 70 cc. of the buffer 

® While this paper was in press, a paper by R. R. Williams and R. Eb. Waterman:  ‘‘Electrodi- 
alysis as a Means of Characterizing Ampholytes,”’ Proc. Soc. Exp. Biol. & Med., 27, 56, 1929 (and 
subsequent related work) came to our attention through the courtesy of Dr. M. K. Brakke. 
Their apparatus consisted of 14 compartments filled with ampholytes in solutions of different pH 
separated by parchment membranes. After 60 hours the concentration was found to be greatest 
in compartments having a pH nearest the isoelectric pH of the ampholyte. Some of the later 
modifications of the apparatus use separate containers connected in series by electrolyte bridges 
This method may be considered as an intermediate step between standard electrophoretic tech- 


nique and the method of ‘isoelectric spectr 





A FURTHER STUDY OF THE APPARENT SYNTHESIS OF VITAMIN By 
BY MAMMARY CANCERS OF MICE 


By D. W. Woo.iey* 
ROCKEFELLER INSTITUTE FOR MEDICAL RESEARCH, NEW YORK 
Communicated January 13, 1955 


In a previous paper! evidence was reported which indicated that mice bearing 
spontaneous Mammary cancers synthesized vitamin By. Nontumorous animals, 
on the contrary, did not seem to form this vitamin in amounts sufficient for normal 
requirements. This finding of apparently greater synthesizing capacity in the 
cancerous animals ran counter to a prevailing belief, which is that cancers may 
differ from normal tissues by lacking some particular power of synthesis. Further- 
more, the finding of the apparent ability of mammary cancers to synthesize vita- 
min By (in contrast to the lack of such ability in the normal animal) provided a 
rational basis for the elaboration of drugs which would poison the tumor without, 
at the same time, harming the host animal.’ * * Such drugs have actually been 
produced.* > Consequently, because of the importance attributable to the finding 
about vitamin By synthesis in the cancerous mice, it became advisable to seek sup- 
porting evidence for such a property of the tumors 

The preceding evidence for the synthesis of vitamin By by mice bearing spon- 
taneous mammary tumors was obtained by animal assay. The young of cancerous 
mothers maintained on the By-deficient diet grew normally, indicating that the 
mothers had acquired and transferred to them enough vitamin to meet their needs 
for normal growth. This excess vitamin By probably was synthesized by the can- 
cer of the mother, or at least by the mother in response to her cancer, because she 
was subsisting on a B,-deficient ration. The evidence to be presented in this paper 
was obtained by microbiological assay. It has been found that tumorous mice 
maintained on a vitamin By-deficient diet contained more vitamin By than did 
similarly treated nontumorous mice. Thus, by direct assay of the tissues by a 
microbiological procedure, the cancerous animals were found to contain more of the 
vitamin than the noncancerous controls. This excess probably arose from synthe- 
sis by the tumor. 

Other explanations for the data can be proposed. One of these is that the 
tumorous animals stored the vitamin more than did the noncancerous controls. 
One may say that, although the diet contained no demonstrable vitamin By, it 
may have carried traces. These traces may have been filtered out and stored by 
the cancerous mouse more efficiently than by the control animal. No adequate 
evidence to disprove such an explanation was found, but some of the data were less 
readily understood on such a basis than on the idea of synthesis of the vitamin by 
the tumor. 


EXPERIMENTAL 


Diets and Animals.—Except where otherwise noted, the By-deficient ration of 
Zucker and Zucker was used. This ration was composed principally of a com- 
mercially available flour made from cottonseed meal and known as “Proflo.”’ It 
was found quite essential to use Proflo made by the old process. The product of the 
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same name made by the new process failed consistently to produce deficient ani- 
mals. Some experiments were carried out with the highly purified ration of Hart- 
man et al., in which alcohol-extracted casein constituted 45 per cent of the mix- 
ture.’ 

All the mice were of the Swiss strain maintained in this institute. These were 
not the Webster Swiss, which is the strain offered for sale by many dealers. As the 
females of the Swiss strain grew old, a rather large percentage of them developed 
mammary cancers. For the present experiments a large number of these old fe- 
males were kept on stock rations and examined weekly. As soon as a mouse 
showed a palpable mammary tumor, she was placed on the deficient diet and fed in 
this manner until the cancer had grown to a large size. The animals were caged in- 
dividually in screen-bottomed boxes and provided with food and water. 

Selection of a Method for Microbiological Assay.—The chrysomonad method was 
the only one found satisfactory for this investigation. Most of the determinations 
were carried out according to the directions of Hutner et al. and Barber ef al.’ It 
was found essential to elimisste completely from the solution to be assayed any 
trace of turbidity due to fat or to suspended particles which may have passed 
through the paper when the unknown was filtered. This was accomplished by 
centrifugation at 14,000 rpm. If this was not done, fictitiously high values were 
frequently found. When the revised procedure was published by Ford,® it was 
adopted for the remainder of the investigation because it reduced the dangers aris- 
ing from tiny particles or from emulsified fat. Before this was done, however, 
many of the samples which had been analyzed by the older procedure were re- 
examined by Ford’s modified method. The results on any sample tested by both 
procedures agreed well. 

The reason for the selection of the chrysomonad method was its greater specificity. 
This organism does not respond to pseudo—vitamins By which are found in the in- 
testines of animals but which show no vitamin By. potency for mammals and birds. 
Because whole mice were to be assayed in the present work, this was an important 
consideration. In the beginning, attempts were made to use the coli-mutant 
method.” To test its suitability, normal Swiss mice were fed the deficient ration 
for 3 days, and individual animals were taken for assay. The whole mouse was 
ground in a Waring Blendor with 10 volumes of water, and the homogenized 
sample was digested, autoclaved, and assayed (cf. Scheid et al.''). Values obtained 
on such samples by the coli method were usually much higher than by the chryso- 
monad method (see Table 1). Furthermore, removal of stomach and intestines 
from the mouse before assaying markedly lowered the value found by the coli 


TABLE ! 
COMPARISON OF B,y» CONTENTS OF NORMAL Mick, AS DETERMINED BY THE CHRYSO- 
MONAD AND Cott MreTHopDs 
Bis CONTENT By CONTENT 
(MILLIGAMMA PER GRAM) (MILLIGAMMA PER GRAM) 
Coli-Mutant Chrysomonad Coh-Mutant Chrysomonad 
Mouse minus stomach 
Mouse: and intestines: 
26 A23 
29 j A24 
65 : A25 
21 8 A27 


86 ( A28 
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method but did not materially change that obtained with the chrysomonad proce- 
dure. These data also are shown in Table 1. This latter finding tended to suggest 
that pseudo—vitamins By accounted for much of the extra vitamin found with the 
coli method. Because the coli values had been corrected for methionine by sepa- 
rate assay in the way recommended," the high results obtained by this method could 
not be attributed to methionine. These preliminary studies were taken as evidence 
that the chrysomonad method was the one which should be used, 

Many experiments were performed to test the reproducibility of the values ob- 
tained on any given sample by the chrysomonad procedure. ‘These indicated that 
the method was satisfactory from this standpoint. For example, mouse A9L was 
found to contain 10.7, 11.0, 11.0, 8.7, and 11.0 milligamma of vitamin By per gram 
in five separate assays. Similarly, recovery of vitamin B,, added to normal mice 
was tested and found to be satisfactory, in that 80-100 per cent of added vitamin 
could be recovered. The data of Table | will also indicate the variation among in- 
dividuals in a group and will show that most of the animals had about the same 
content of vitamin By. This was considered of importance in the evaluation of 
subsequent results. In all the analyses subsequently described in this paper, the 
stomach and intestines were removed from the mice before they were ground for the 
assay. 

Vitamin By Content of Nontumorous Mice Fed the Deficient Diet.——-Groups of three 
to eight adult female Swiss mice were fed the deficient ration in the manner in- 
dicated above, and, after various intervals of time, each animal was assayed. 
The data of Table 2 show that the average By. content was greater if the mice had 
been receiving the deficient ration for a short time but that, after they had eaten it 
for 2.5 weeks, no further significant change in By content was detectable. 

Similarly, the effect of age of the mouse at the start of the feeding of the deficient 
diet was studied. This was important, because most of the tumorous mice would 
be old and of varying age. It was found that the animals approached the same 
content of vitamin By when they were given the deficient ration, regardless of 
whether they were 2 months of age or 9 months of age when they were started. 
In other words, any normal adult mouse approached a vitamin By content of about 
9 milligamma per gram when it was fed this deficient diet for several weeks. 

Vitamin By Content of Mice Bearing Spontaneous Mammary Cancers when Fed 
the Deficient Ration.—Adult Swiss mice bearing small spontaneous mammary 
tumors were fed the deficient ration until the cancers were large. The time re- 


quired for this was usually about a month, but in some instances it was consider- 
ably longer. They were then analyzed in the way described. At the same time, 
normal Swiss mice of the same age were tested in the same way. The results, 


shown in Table 3, indicate that the cancerous mice contained about twice as much 
of this vitamin as did the controls. 

Vitamin By. Content of Normal, Pregnant Mice Fed the Deficient Diet.—Bevause 
many differences of cancerous mice from normals have, in the past, been shown to 
be due only to the resemblance of the cancer to embryonic tissue, it was considered 
advisable to determine whether this same criticism might be leveled against the 
difference with respect to vitamin By. A group of twelve normal Swiss mice in the 
first week of pregnancy was fed the deficient ration until the day of parturition or 
until a few days before this time. Each mouse was then analyzed for vitamin By. 
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The average content was 9.5 milligamma per gram. This value, when compared 
to 9.4 milligamma per gram for eight nonpregnant controls, showed that there was 
no evidence for synthesis of the vitamin by the embryos. 

Vitamin By, Content of Spontaneous Mammary Cancers in Comparison to the Vita- 
min By Content of the Host Mice from Which They Came.—When a cancer-bearing 
mouse was fed the deficient diet for several weeks, and its tumor was excised and 
assayed at the same time as the remainder of the carcass, the content in the carcass 
was usually found to be somewhat greater than that in the tumor. Thus in mouse 
275 the carcass contained 23 milligamma per gram, and the tumor 10 milligamma 
per gram. The corresponding average values from six such animals were 30 and 
15. 

Destruction of Vitamin By by Excised Spontaneous Mammary Cancers.—One 
reason for the lower By.-content of the cancers in comparison to the remainder of the 
carcass Was found to be that the cancer, when minced in the presence of the vitamin, 
destroyed some of it. Presumably, a similar destruction of endogenous vitamin 
occurred when the cancer was ground for assay. In a typical experiment, mouse 
AS4 was maintained on the deficient ration, and her cancer, weighing 2.6 gm., was 


TABLE 2 TABLE 3 
AVERAGE By. CONTENT OF NORMAL VivaMIn By Contrent or De- 
Mice Fep tHe Dericient RATION FICIENT Mice BEARING SpPon- 
FOR VARYING PERIODS OF TIME TANEOUS MAMMARY CANCERS, 
Peete A INCLUDING Dat A ON THE Er- 
Period on No Content FECT OF EXCLUSION OF THE 
Deficient Ration of Milligamma TUMOR FROM THE SAMPLE 
Mice per Gram 
2 Average Bu 
- a Content 
6 No (Milligamma 
~ g 4 of per 
6 Description of Mice Animals Gram) 


Greater than Nonearcerous 24 o3 
10 6 Q! Cancerous 13 17 
Cancerous, 
minus 
tumor 13 28 


excised. ‘Half of it was minced with scissors in 13 ec. of water and then ground 
with sand. The mixture was incubated for 1 hour at 37° C. The suspension was 
then autoclaved, and the clear supernatant liquid was assayed. ‘Two and one-half 
milligamma of vitamin By per gram of fresh tissue were found. The other half of 
the cancer was mixed with 20 milligamma of vitamin By per gram of fresh tissue 
before mincing, but other operations on it were the same. It was found to contain 
13 milligamma per gram, whereas it should have contained 22.5. Clearly, some of 
the added vitamin had not been recovered. 

Each tumor examined in this fashion was not found to destroy the same amount 
of vitamin By». Some destroyed very little (see Table 4). This variation among 
individual cancers may explain why the By» content of a few cancer-bearing individ- 
uals was found to be much less than the average for such mice. If a cancer was 
very active in destruction of the vitamin, the value for By content found after 
homogenization of the mouse plus cancer would be unusually low. 

Vitamin By Content of Cancerous Mice Analyzed without Their Cancers.—The ex- 
periment on the By. content of cancerous mice maintained on the deficient ration 


was repeated, except that the cancer was removed from each mouse before it was 
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ground for analysis. The data in Table 3 show that under these conditions the 
animals were found to contain considerably more of the vitamin than if their tu- 


mors had not been removed. The results of this experiment indicated that the 
cancer tissue was destroying some of the vitamin By when carcass plus cancer were 


ground together before the analysis. Nevertheless, even with this technique of dis- 
carding the cancer before grinding, the By content of a cancerous mouse might be 
as high as 55 milligamma per gram, while others might show only 7 milligamma per 
gram. It may be that various strains of cancer differ quantitatively in ability to 
synthesize vitamin By as well as to destroy it. 


EXPERIENCES WITH TRANSPLANTED MAMMARY CANCERS 


a) Vitamin By Synthesis (as Suggested by Content) of Mice Bearing Cancer Strain 
W137.—As explained more fully elsewhere,> transplanted mammary cancer strain 
W137 was obtained from one of the Swiss mice of the present study and was trans- 
planted in Swiss mice. The experiments on By synthesis were performed with the 
fourth to eighth passages of this strain. This strain. was of interest because it was 
the most susceptible to antimetabolites of dimethyldiaminobenzene and should, 
therefore, show evidence of By synthesis if the general working hypothesis * has 


TABLE 4 
DesTRUCTION OF VITAMIN By BY SPONTANEOUS CANCERS IN ViTRO* 
Bie CONTENT 


MILLIGAMMA PER GRAM PeRCENTAG® 
Mouse Tumor Tumor and Br DESTRUCTION 


AS4 2.5 13.0 50 
146 19.0 37.6 7 
275 10.0 21.0 45 
277 21.0 35.0 30 
* The tumor tissue was incubated with 20 milligamma of vitamin By per gram of fresh tissue for 
one hour at 37 
validity. Three mice bearing established transplants of W137 were fed the de- 
ficient ration for one week. The tumors were discarded, and the carcasses were 
analyzed for vitamin By. The average value was 30 milligamma per gram, which 
was much more than in the controls (ef. Table 1). 

A second set of experiments was performed in this way, except that the Proflo 
ration was replaced by the By-deficient diet of Hartman et al.’ This was a semi- 
synthetic diet rich in aleohol-extracted casein. The nontumorous control mice had 
an average By content of 18 milligamma per gram. The mice bearing transplanted 
cancer strain W137 contained 24 milligamma per gram. Evidently this ration did 
not render mice so deficient in vitamin By as did the Proflo diet and therefore did 
not magnify so much the difference in the two kinds of mice. 

b) Synthesis versus Storage in Strain W137.—One might contend that by using 
well-established cancer transplants one had given opportunity to the host plus 
cancer to store extra vitamin By before the deficient ration was fed. The resultant 
greater content of the vitamin in the animal thus might represent not synthesis 
but, rather, merely more tenacious storage. This point was tested in the following 
way. 

Swiss mice were fed the Proflo-containing deficient diet for two weeks. Half of 
them were then implanted with strain W137, and the others were maintained as 
controls. Both groups were continued on the deficient ration, and, when the tu- 
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mors were well developed (2-4 weeks), analyses were performed. The average By 
content of the controls was 10 milligamma per gram, and of the cancerous mice, 
15.5 milligamma per gram. This result tended to suggest that the tumor strain 
W 137 transplanted into animals depleted in vitamin By led to accumulation of this 
vitamin in the carcass, presumably by synthesis. However, a garnering and stor- 
age of traces of the vitamin which may have been in the diet cannot be entirely ex- 
cluded as an explanation. 

c) Failure of Transplanted Strain W113 To Gwe Any Evidence for Synthesis.—In 
the preceding paper on the evidence for By synthesis by spontaneous mammary 
cancers,' the failure to find any such synthesis with a transplanted strain was re- 
corded.'? In the present study, another strain of transplanted mammary cancer 
of Swiss mice has likewise been found to give no evidence for synthesis of vitamin 
By. Clearly, then, some transplanted strains do and some do not give evidence for 
synthesis. 

Swiss mice were implanted with strain W113,5 and 1-2 weeks: later, when the 
tumors were established, the animals were fed the deficient diet. Two weeks 
later (when the tumors were large) the animals were killed and their tumors re- 
moved and discarded. Analyses of the carcasses of six of these mice showed them 
to contain, on the average, 5 milligamma of vitamin By per gram. Therefore, the 
vitamin content was somewhat less than in nontumorous controls (cf. Table 1). 

d) Active Destruction of Vitamin By by Transplanted Tumor Strain W113.—The 
strain which failed to give evidence for synthesis was found to destroy vitamin By 
actively. Eight mice implanted with strain W113 were treated exactly as described 
in the preceding section, except that the tumors were not removed from the animals 
before the analyses were performed. The average By content of these mice was 2.1 


miligamma per gram. Because the nontumorous controls had more than three 


times this amount, it was clear that the cancer had destroyed much of the By in 
the sample when the whole mouse was ground. This experience with strain W113 
gave the most clear-cut example of By destruction which has been seen in this 
work and led to the trials on By destru¢tion by spontaneous cancers which were 
described earlier in this paper. In the case of the spontaneous cancers, the de- 
struction of the vitamin was not found to be as extensive as with this transplanted 
strain, 

Direct proof of destruction of vitamin By by the excised cancer was obtained in 
the way described for the similar experiments with the excised spontaneous tumors. 
With strain W113, 1 gm. of minced cancer tissue was incubated 1 hour at 37° with 
10 milligamma of By. Analysis showed that only 2.9 milligamma of the vitamin 
remained. Since the tissue without added By» gave a value of 1.2 milligamma per 
gram, 74 per cent destruction of added vitamin had taken place. 

Four strains of transplanted mammary cancers were examined in the way just 
described for their abilities to synthesize vitamin B,. while growing in Swiss mice 
and to destroy it when incubated in vitro. Evidence for synthesis has been ob- 
tained for three of the strains, and evidence for destruction (at least to some extent) 
with all four. The most active destroyer was W113, just described, for which no 
evidence for synthesis was found. 
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DISCUSSION 

Although the data of this paper serve as additional evidence for the idea that 
spontaneous mammary cancers of mice synthesize vitamin By, they do not establish 
this idea as a fact beyond question. It has not been possible to eliminate the alter- 
native that the tumorous animals trap and store traces of this vitamin which may be 
present in the deficient ration. The difference from normal animals would, if this 
were the true explanation, then reside in a greater capacity of the cancerous ani- 
mals for trapping and storage of vitamin By. However, the present evidence ob- 
tained with microbiological assays corroborates that found previously by animal 
assay in showing that the cancer-bearing mice either synthesize or store more By 
than do normal controls. 

The destruction of vitamin By» by certain mammary cancers is noteworthy. 
Some enzymic system was possibly responsible for this destruction, since the ac- 
tivity was lost by heating the tissues. Prior to this finding, no evidence for a By- 
destroying system has come to our attention. To be sure, there are preparations 
from normal gastric juice and from other sources which combine with, and thus 
inactivate the microbiological activity of, this vitamin. However, this activity 
may be regained by heating the reaction mixture. In the present case, however, 
the activity could not be regenerated by heating and may represent true destruction 
of the molecule. 

Some may wonder whether the destroying system is anything except the By-syn- 
thesizing system acting in reverse. This may well prove to be the case. However, 
present evidence does not permit any conclusions on this point. 

The variations among individual tumorous mice both in the amount of apparent 
synthesis of vitamin By and in By-destroying capacity were striking. Thus, al- 
though control animals maintained on the deficient diet each had approximately 
the same By content, the variation among cancerous mice was rather large. Some 
showed five times the average nontumorous content; one had less than the average 
of the controls. The majority of the cancerous individuals did, however, cluster 
around a By content which was twice that of the controls. This failure of a few of 
the cancerous individuals to show evidence for synthesis was found not only by the 
methods of the present study. Careful examination of the records of the older in- 
vestigation! showed that there, too, was an individual which gave no evidence for 
synthesis. Clearly, then, there are some mice which, when invaded by the cancer, 
do not show evidence for By formation. In the case in which the cancer was 
examined in a transplanted condition (i.e., strain W113), failure to find any evi- 
dence for synthesis was associated with the greatest capacity for destruction of the 
vitamin. It remains to be determined whether causal relationship existed. 


SUMMARY 
A microbiological method for the determination of vitamin By-content of mice 
was sought. The coli-mutant method was found to be unsatisfactory in comparison 
to the chrysomonad method. When the chrysomonad method was employed, 
analyses of nontumorous mice fed a By-deficient ration showed them to contain on 
the average, 9 milligamma of the vitamin per gram of live weight. The correspond- 
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ing value for similar mice bearing spontaneous mammary cancers was 17 milli- 
gamma. The extra By was thought to have arisen from synthesis by the tumor. 
Other explanations were also considered. A transplanted strain of these cancers 
was found to give similar evidence for vitamin By synthesis. A second trans- 
planted strain showed no evidence for such synthesis but, instead, was demon- 
strated to bring about active destruction of the vitamin in vitro. Several individ- 
ual spontaneous mammary cancers were likewise shown to destroy this vitamin, 
but not so actively as did the one transplanted strain. This destruction was shown 
to be different from mere combination of a protein with the vitamin, such as is 
known to occur in gastric juice. 
* With the technical assistance of V. Armbrust and G. Schaffner 
' 1D. W. Woolley, these PRocEEpINGs, 39, 6, 1953 
2D. W. Woolley, J. Exptl. Wed., 93, 13, 1951 
D. W. Woolley, Proc. Second Intern. Biochem. Congr. (Paris), Julv, 1952 
D. W. Woolley, Cancer Research, 13, 327, 1953 
*D. W. Woolley and G. Schaffner, Cancer Research, 14, 802, 1954 
6, M. Zucker and T. F. Zucker, Arch. Biochem., 16, 115, 1948. 
A. M. Hartman, L. P. Dryden, and C. A. Cary, Arch. Biochem., 23, 165, 1949. 
‘S$. H. Hutner, L. Provasoli, and J. Filfus, Ann. V.Y. Acad. Set., 56, 852, 1953; F. W. Barber, 
): L. Baile, C. B. Troescher, and C. N. Huhtanen, Ann. N.Y. Acad. Sei., 56, 863, 1953 
* J. E. Ford, Brit. J. Nutrition, 6, 324, 1952 
B. D. Davis and E. 8S. Mingioli, J. Bacteriol., 60, 17, 1950 
‘tH. BE. Scheid, M. M. Andrews, and B.S. Schweigert, J. Nutrition, 47, 601, 1952 
‘2 This, in fact, was used as an argument against the idea that the apparent synthesis in te 
spontaneously tumorous animals was due to bacterial infection or necrosis of tissues 


A GENERALIZED THEORY OF PLASTICITY INVOLVING THE VIRIAI 
THEOREM* 
By Henry EyrinG ano Tarkyvk REE 
UNIVERSITY OF UTAH 
Communicated December 22, 1954 
The theory of plastic deformation proceeds by considering a series of relative 


displacements of patches lying on the two sides of each shear surface. Displace- 
ments along such shear surfaces occur as relaxations obeying a generalized absolute 


reaction rate equation. A general statistical formulation has been achieved which 
fits the usual cases of shear. Application of these results to various types of cases 


will appear elsewhere. A generalization of the absolute reaction rate equation is 
also presented. 

A calculation of the probability, ?,, that a molecular system, A, will have the 
total energy, ¢,, distributed in some unique way among its degrees of freedom. 
where A is in equilibrium with a large number of oscillators whose mean available 
energy is y, yields the result! 


P, ce eT . (1) 


The appropriate value to take for y depends on the circumstances. There are a 
yariety of important cases involving different types of degrees of freedom. 
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Case 1: For degrees of freedom at equilibrium ¥ will have the value &T, where 
7 is the temperature of the system as a whole. 

Case ll: If some irreversible rate process is in operation, it may be that for cer- 
tain degrees of freedom the local statistical y may be well approximated by kT), 
where the subseript / indicates the local value. Thus the appropriate 7, may be 
considerably different from the over-all measured temperature 7’. 

Case II]: If we have a frozen-in equilibrium with respect to certain degrees of 
freedom, tt may be convenient to write for them y, = kT, Here 7’, is the tem- 
perature at which the equilibrium was frozen in and may have nothing to do with 
the ambient temperature, 

Case 1V: If, after an equilibrium is frozen in, certain components are arbitrarily 
removed or modified, we may be obliged to give up the frozen-in equilibrium con- 
cept and replace it by a detailed accounting for the individual species present. 

Case V: Inthe same way it may be that there is no local temperature, 7;, which, 
combined with equilibrium theory, will properly describe all the relative concen- 
trations in certain groups of states. In this case, the group of states must be 
broken into subgroups until such a local temperature can be properly applied to 
each of the subgroups. This final subdividing into equilibrated subgroups may 
end with a single quantum state in a subgroup. 

Now if, beginning with equation (1), we go through the usual quasi-equilibrium 
derivation of the specific rate constant, k’, we obtain 
Kyo ~ Qo (art 

€ 


Ik’ 
h 


Here yo is the mean available energy of those oscillators equilibrating with the 
translational degree of freedom corresponding to crossing the potential barrier, and 
y, is the same mean available energy of those oscillators equilibrating with the 7th 
group of activated complex degrees of freedom requiring the activatiqnal free 
energy AF, ft. If there are classical degrees of freedom equilibrating with all the 
degrees of freedom of the activated complex, each y can be replaced by kT, and we 
obtain the familiar rate expression, 


k’ ; saath (3) 


Various kinds of incomplete equilibration of activated complex degrees of free- 
dom are observed in plastic flow. Consider that shear occurs along sets of parallel 
planes. The distance between these planes is indicated by the symbol A,. The 
shear force per unit area is indicated by f. A patch of atoms or molecules whose 
cross-sectional area is AsA; shift or jump as a unit on either side of the shear plane. 
We obtain, then, for the rate of shear S, when the system is in equilibrium with an 


external temperature 7’, the expression® * 


ee ee fr\rr rd KT . 2 Aa. fro r _ 
¢ 2 sinh —— a K exp 4 — es! ae 2 sinh =a ( 
: 27 hy kV OR ee ee 
Here k’ is the frequency of displacement of the flow patch per second in the for- 
ward and in the backward flow direction when f = 0. The distance the patch 


moves per jump Is A, 
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The generalized equation when local degrees of freedom fix y becomes 


‘2 sinh — 
ae y, I 27; 


. — 


N Ko oe Agee) 2. hodAS 
KY a ‘ % 2 { J 


ry h 


Equation (5) reduces to equation (4) when yo = y; = y, = kT. When this is not 
true, the problem remains as to the appropriate values to take for the y’s. 
The Simple Virial Theorem Applied to Viscous Flow.—-Consider the usual simple 
equation 
d 
(mx) == m&x + mx?*. (6) 
c 
Here m, x, X, and ¥ are the mass, distance, velocity, and-acceleration, respectively, 
associated with a given degree of freedom. If both sides of identity (6) are in- 
tegrated from the initial time ¢, to some final time ¢,, there is obtained the result 


(mXx)ty — (mex) = (mXx + mx*) (t, — ty) (7) 


l 


, \(mx)., — (mx), = mtx + mx? = 0. (8) 
PS ae 


The final zero value follows from looking at the material to the left of the first 
equality in equation (8). This must be zero, since the difference inside the braces 
fluctuates about zero, whereas | /(t; — ¢,) approaches zero for averages over long 
time intervals. Since m*¥ = X, where X js the force in the zx direction, we write the 
usual virial expression‘ 

Xr + mx? = 0. (9) 


Here the average is over a very long time. If equation (9) be applied to the reac- 
tion co-ordinate for plastic flow, the following identifications can be made: 


X = fod2ds (10) 


and 


x = gy, (11) 


where g is the ratio of the length of free vibration, 2, of the flowing patch whose area 
is AAs to the distance jumped, \. Thus g should approximate! the ratio of sound 
velocity in the'solid to that in the gas—a value in the neighborhood of 7. Here f, 
is the local microstress to be distinguished from the over-all macrostress, f. Fur- 


ther, 
mx? = ¥;. (12) 
We now define AF,* as the mean work required to displace unit area along the 


slip p'ane a distance \/2. Whence the free energy, AF,*, to move the patch of area 
AeA; the distance \/2, i.e., to the top of the potential barrier, is 


ae 


AF,* = AF,*naAs. (13) 
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Combining formulas (9), (10), (11), and (12) yields 
fodAodsgr — Vue 


Substituting equations (13) and (14) in equation (5) gives 


1 +f » 
= AF,” AP. : 
x - » 2 sinh —. (15) 
1 v1 3 =forgS 2afo 


As far as the available evidence indicates, it is usually satisfactory to write yo = 
kT. However, to the approximation that yo can be identified with the average 
energy of an oscillator, we have 


Yo _ l + 4 hy 
h hi 2 hy 


é — | 
AT 
at low temperatures this has the value ' » y and at high temperatures the usual 
value, kT h. Rate processes involving electrons and He'' seem at least in quali- 
tative accord with the low-temperature form of equation (16), since there is no 
evidence that processes cease at absolute zero in the way indicated by the factor 
kT h. Because of quantum effects at low temperatures, the transmission co- 
efficient, x, including barrier leakage, will be important, and Kyo/h will be best 
treated asaunit. Finally, we can write 


1 + I wos ’ 
— AF, AF, . AF, Pa 
pi tore eee, (17) 

i=1 ¥; RT, R17 

Here AF,* is the free energy of activation for these degrees of freedom when they 
were frozen in at 7’,, while the AF,* is the free energy of activation for the degrees 
of freedom which are in equilibrium under the conditions of the experiment. 

This leads us to the general equation for shear rate, 
. \. sF sah ars f 


K exp \— 2 sinh 
Ay h ' 


RT, RT forg | 2ofo 
If we write AF,? AH,* — TAS,*, we see that (18) can also he written as 


VN KT (As Ak: f 

. 0 
= exp — ~~~ > 2 sinh ; 
ha ( R RT \ gfe 

where we have for AS,;,, the apparent activation entropy, the value 

e re ee AF; RAF, 

ASin49 = 483+ Rink - —— -— ’ 

T', JoAg 


or, even more simply, we can define an effective activation entropy 


b a N 
ASi, = ASi,, + Rin 7 
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and write for the rate of shear 


kT AS;i, AH,’ f 
exp oe ly sinh ——. (22) 


R RT‘ ~ 2Gfo 


t 


The quantity 1/g A/a (see eqs. [11] and [22]) has the value 6.32, according to 
Tobolsky and Eyring’s analysis® of Trumpler’s data for stress relaxation of steel.° 


According to the theory, 1/g should approximate, as it does, the relative velocity 
of sound in the solid and the hypothetical vapor, provided that the distance, A, 
which a patch jumps approximates a lattice distance. In this case, for steel, fo is the 
initial stress of the steel. Further, Tobolsky and Eyring found that AF? decreases 
with increasing stress fy in accordance with the term AF,* /fuAg of equation (18). 
These authors also noted that 1g is a constant for the creep of rubber at constant 
load. The form of the argument of the hyperbolic sine in equation (15) is of con- 
siderable interest. The interpretation of g has been given. In stress relaxation the 
local microstress, fo, has been identified as the initial macrostress. The independ- 
ence of temperature of the argument of the hyperbolic sine has not previously 
been explained but has been widely observed.* *7 It has been found that equa- 
tion (19), with AH,* 0, gives an excellent account of the plasticity of many 
metals at temperatures far below their melting points.’ 

* This work was supported in part by the Office of Naval Research. 

1H. Eyring, 8. Walter, and G. FE. Kimball, Quantum Chemistry (New York: John Wiley & Sons, 
1944), pp. 296-298. 

2H. Evring, J. Chem. Phys., 4, 283, 1936 

* A. Tobolsky and H. Eyring, J. Chem. Phys., 11, 125, 1943. 

‘Cf Hy ring, Walter, and Kimball, op. cit., pp. 355-358. 

’W.E, Trumpler, Jr.,./. Appl. Phys., 12, 248, 1941 

®S.S. Kistler, /. App!. Phys., 11, 769, 1940; W. Kauzmann, Trans. Am. Inst. Mining Met. Engrs., 
143, 57, 1941; A. V. Tobolsky and R. D. Andrews, J. Chem. Phys., 13, 3, 1945; T. Ree and H 
Eyring, University of Utah Tech. Rept. (ONR), 28, December 1, 1952; 0. D. Sherby and J. E. Dorn, 
Trans. Am. Inst. Mining Met. Engrs., (J. Metals), 194, 959, 1952; J. Metals (Trans.), 5, 324 


1953. 
7 J. W. Frederickson and H. Eyring, Am. Jnst. Mining Met. Engrs. Tech. Rept., 2423, 1948. 


INVESTIGATIONS ON LIGNIN AND LIGNIFICATION. XV. 
HETEROGENEITY OF NATIVE AND ENZYMATICALLY LIBERATED 
LIGNINS AS’ ESTABLISHED BY ELECTROPHORESIS AND PAPER 

CHROMATOGRAPHY 
By WALTER J. ScHUBERT* AND F. F. Norp 
DEPARTMENT OF ORGANIC CHEMISTRY AND ENZYMOLOGY, T FORDHAM UNIVERSITY, NEW YORK 


Communicated by M.L. Wolfrom, December 14, 1954 


As it occurs in the natural state in wood, lignin is a complex polyphenolic com- 
pound, and it is generally assumed that it is formed through the polymerization 
of some fundamental unit, such as the phenylpropane moiety ‘This polymeriza- 
tion may occur through various possible linkages whose natures are as yet not es- 


tablished. ! 
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However, it is certain that the lignins isolated from various species of wood or 
other lignified tissues possess quite different properties. Thus chemical and phys- 
ical investigations from this laboratory have revealed significant individual dif- 
ferences among the native lignins of softwoods,? hardwoods,’ cork,* and bagasse.® ° 
Other native lignin preparations were investigated by Schuerch,’ who reported 
native spruce and aspen lignins to be almost completely soluble in weak alkali 
and therefore not subject to partition under the conditions of countercurrent dis- 
tribution. He concluded, therefore, that his samples were homogeneous. 

However, it is still not inconceivable that the lignin of any one species of wood 
alone is not homogeneous but may present minor differences with respect, for ex- 
ample, to the nature of the linking of the phenylpropane groups to one another. 
Furthermore, differences may also arise from a varying degree of polymerization 
of the lignin; that is, the number of aromatic units actually present in each indi- 
vidual lignin particle may vary. 

The term “lignin” thus might imply a mixture of polymers, all possessing a 
very similar constitution but with the possibility of small differences in molecular 
size and chemical structure.® 

The study of the chemistry of such a substance is greatly handicapped by the 
difficulty of separating the various fractions and the uncertainty as to whether or 
not these fractions are real chemical entities, especially if the compounds are 
polymeric. 

Current chemical methods of separation, such as fractional crystallization, dis- 
tillation, chemical reactivity, or countercurrent distribution, are of but limited 
applicability to such conglomerations, because of the close similarity of the mem- 
bers involved; hence we have had recourse to certain physicochemical techniques, 
such as electrophoresis and partition chromatography. 

Experimental.—-The native lignins studied were obtained by thoroughly re- 
moving “extractives” from the ground, air-dried wood species with water and ether. 
The sawdust was extracted at room temperature with ethyl alcohol, until the alco- 
holie solution no longer responded positively to the phloroglucinol-HC! color test 
for lignin. After removal of the alcohol by distillation at reduced pressure, a resi- 


hous material was obtained. This was purified by dissolving in dioxane and pre- 


cipitating alternately into water and ether until a constant methoxyl value re- 
sulted.® 

Enzymatically liberated lignins were obtained by decaying ground wood (pre- 
viously freed of its native lignin) with various “brown rot” fungi, which preferen- 
tially attack the wood cellulose. After a certain period of decay, the mycelia of the 
fungi were removed, the wood collected and dried, and the liberated lignin extracted 
and purified according to the procedure for native lignin." The percentage yield 
of enzymatically liberated lignin obtainable is directly related to the duration of the 
fungal decay. The liberation of the total amount of lignin in a wood sample is there- 
fore dependent only on the duration of the enzymatic\degradation of the cellulose 
associated with it. 

For the paper-chromatographic investigations | per cent solutions in dioxane 
were prepared of each of the following native and enzymatically liberated lignins: 
white Scots pine, oak, birch, maple, and bagasse. 

Strips of Whatman No. | filter paper were spotted with approximately 0.02 ml. 
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of one of these solutions, and the spots were developed with the solvent mixture 
isobutyl alcohol: benzene: water, 5:45:50. The strips were dried, sprayed with 
dilute HCI, redried, and then resprayed with an alcoholic solution of phloroglucinol. 
The red-violet coloration typical for lignin developed at the origin and at points along 
the chromatogram. The FR; values of these spots were measured in the customary 
way. 

Results and Discussion.—The chemical compositions of the native and enzy- 
matically liberated lignins investigated are presented in Table 1. The differences in 
elementary composition and methoxy! content of the lignins from varying species 
are apparent from these data. 


TABLE 1 
CHEMICAL COMPOSITIONS OF LIGNINS (PER CENT) 


Lignins Cc “Hs Lignins 
White Scots pine N.L.* 64.0 4 Birch E.L. 
White Scots pine k.L.f 64.2 Maple N.L. 
Oak N.L. 58.6 6.3 Maple E.L. 
Oak E.L. 58.4 ae ) Bagasse N.L. 
Birch N.L. 61.4 5.! ( Bagasse E..L. 


* N.L. = Native lignin 

t E.L. = Enzymatically liberated lignin. 

On electrophoretic analysis, some of these lignin preparations appeared per- 
fectly homogeneous, whereas other samples presented a small boundary of lower 
mobility next to the main boundary. For example, homogeneous material is ob- 
tained from the native maple and cork lignins. On the other hand, native white 
Scots. pine lignin reprecipitated several times presents a well-marked trailing boun- 
dary." 

It is worthy of note that the same trailing boundary is also observed in samples 
of enzymatically liberated white Scots pine lignin. Native bagasse lignin also 


appears to be somewhat inhomogeneous, presenting a detectable second boundary. 


Enaymatically liberated bagasse lignin behaves similarly.'! , 
The results of certain of these electrophoretic analyses are summarized in Table 
2. Here the mobilities of the main boundaries only are presented. It is evident 


that under similar conditions of pH and buffer composition all the lignin samples 


present comparable mobilities. ‘ 


TABLE 2 
ELECTROPHORETIC Mosiuity or Various LicNins In GiycinE-NaCl Burrer, pH 10.7 


Mosiritry (Cmi/V Sree X 10° 
LIGNINS Ascending Descending 


White Scots pine N.L. —10.2 —10.1 
Bagasse N.L. —9.§ —i-% 
Maple N.L. —9.. —8.3 
Oak N.L. —s =$.2 
Cork N.L. —9 ! =e 


In view of the complexity of lignin and the possibility of variations in its struc- 
ture, it was significant to find that the majority of our lignins gave patterns char- 
acteristic of electrophoretically homogeneous compounds, with the white Scots pine 
lignins being exceptional. 

It was of particular interest to observe that this lignin, whether native or en- 
zymatically liberated, presented an identical pattern, and, on mixing the samples, 
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the electrophoretic patterns did not change. Therefore, it must be assumed that 
this lignin is a mixture of electrophoretically distinct components, whose ease of 
extraction is the same both before and after enzymatic decay. 

In view of the electrophoretic inhomogeneity of the native and enzymatically 
liberated white Scots pine and bagasse lignins, we attempted to apply the paper- 
chromatographic technique of separation to all our lignin preparations. 

The R, values of the spots developed are recorded in Table 3. The results indi- 
cate a definite inhomogeneity with respect to the native and enzymatically liber- 
ated lignins of white Scots pine, maple, and bagasse. Oak and birch lignins gave 
no indication of inhomogeneity. 


TABLE 3 
PAPER-CHROMATOGRAPHIC ANALYSIS OF LIGNINS 
Ry Vatues Re VaLves 


Immobile Mobile Immobile Mobile 
LiGNIN® Spot Spot LIGgNINS Spot Spot 


White Scots pine N.L. 0.0 0.81 Birch E.L. Nil 
White Scots pine EL 0.75 Maple N.L. 0 
Oak N.L. Nil Maple EL. 0 
Oak E.L. Nil Bagasse N.L. 0.7 
Birch N.L. Nil Bagasse E.L. 0 


Applying these R, values, it was possible to cut additional spotted, developed, 
but unsprayed strips into small portions, containing the mobile and immobile 
spots. In addition, a third portion of the strips was taken to serve as a blank for 
any extractable material which might be present in the paper. The mobile, im- 
mobile, and “blank” portions of several strips were then eluted separately by ex- 
traction with ethy! alcohol in Soxhlet extractors. The solutions were dried over 
anhydrous sodium sulfate and filtered. 

The ultraviolet absorption spectra of the solutions were then determined, with 
the “blank” solution serving as the standard solvent. The results for the two frac- 
tions from native white Scots pine lignin are shown in Figure 1. The immobile spot 
of this preparation presents a spectrum which is evidently the characteristic lignin 
curve. However, the mobile spot from this lignin shows only a shoulder at about 
255 Mu. 

The spectra of the two fractions from the enzymatically liberated white Scots 


pine lignin reveal essentially the same situation (Fig. 2). The proximity of the 
curves of the mobile spots of the native and enzymatically liberated lignins suggests 
a definite similarity between the components of the two ‘‘mobile spots.” 

The native and enzymatically liberated maple and bagasse lignins behaved in 
a similar way. That is, in all cases, the immobile spots present curves character- 
isticaliy associated with the lignins proper,® while the mobile spots showed curves 


only with end absorption. 

Commenis.—A comparison of the results obtained by the electrophoretic and 
chromatographic methods of analysis reveals a satisfactory corroboration. Thus 
both the native and enzymatically liberated samples of oak and birch lignin behave 
as homogeneous preparations in the two methods, while the native and enzymati- 
cally liberated samples of white Scots pine and bagasse lignin give evidence of 
secondary components in both methods. However, the fact that any one method 
of analysis alone might fail to detect minor components in certain cases is indicated 
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Fig. 1.—U.V. spectra of mobile and immobile spots of native white Scots pine lignin 
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by the homogeneous behavior of native and enzymatically liberated maple lignins 
in electrophoresis, as contrasted with the indication of a secondary component in 


both these lignins by the paper-chromatographic technique, 

As a result of this study, lignins should be regarded as mixtures of components, all 
possessing similar structures but with the possibility of certain chemical differ- 
ences, as discussed above. It is for such mixtures of substances which may differ 
both in structure and in molecular size that Staudinger introduced the term Grup- 
penstoffe."* 

* Procter and Gamble Fellow. 

t Communication No. 293. Presented before the Division of Biological Chemistry of the 
American Chemical Society, New York, New York, September 13, 1954. This work was sup- 
ported by the National Science Foundation and the Office of Naval Research. For Communica- 
tions Nos. XII, XIII, and XIV of this series see notes 5, 8, and 11. 
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THE COSMIC ABUNDANCES OF POTASSIUM, URANIUM, AND 
THORIUM AND THE HEAT BALANCES OF THE EARTH, 
THE MOON, AND MARS* 


By Harouip C. Urey 
INSTITUTE FOR NUCLEAR STUDIES AND DEPARTMENT OF CHEMISTRY, UNIVERSITY OF CHICAGO 
Communicated January 24, 1956 


In a discussion of the abundances of the elements the writer' did not select 
values for the abundances of uranium and thorium because of the great variability 
in the older data, but especially because the more recent values for uranium and 
thorium contents of meteoritic material seemed to be much too high to permit an 
understanding of the heat balances of the earth and moon. Chackett, Golden, 
Mercer, Paneth, and Reasbeck? found 0.106 and 0.335 p.p.m. for the mean content 
of uranium and thorium in the Beddgelert meteorite, and Davis* found somewhat 
smaller values in achondrites. Such values would require that ‘much more heat 
has been generated in the moon and earth than seems likely to me, and I have been 
unable to think of any reasonable process to account for less uranium in the moon 
and earth than in the meteorites.”! The difficulties represented by the chemical 
homogeneity of Mars‘ should have been included, for it seems most probable that 
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Mars would have melted during geologic time if it had contained the reported 
amounts of these elements and hence would have formed a core, and this appears 
not to be the case on the basis of present data. Since this paper was written, 
many other data have accumulated which indicate that the general level of the ura- 
nium abundances in the Beddgelert meteorite is indeed representative of the chon- 
dritic meteorites. The amounts of uranium in the iron meteorites as reported by 
Dalton, Golden, Martin, Mercer, and Thompson® and by Davis’ are of the order 
of a few per cent of this value for Beddgelert. Since uranium is very electroposi- 
tive,® the element is probably not dissolved in the metal but occurs in occlusions 
which are known to be present in these objects. A few per cent of such otclusions 
seems probable, but some 10 or 20 per cent, as would be required if the abundance 
of uranium were one-tenth that reported for Beddgelert, i.e., about that reported 
for the Modoc meteorite by Patterson et al.,7 1.1 & 107° gm/gm, would be an 
impossible assumption. Either uranium is dissolved by metallic iron for reasons 
that are not understandable on the basis of its chemical properties, or the analytical 
data on the iron meteorites are incorrect for reasons that are far from obvious, or 
some material containing uranium in about the concentrations indicated for 
Beddgelert are occluded in iron meteorites to the extent of a few per cent. Recent 
results of Patterson® on the uranium and lead contents of several chondritic meteor- 
ites indicate that their uranium contents are about the same as that secured for 
Beddgelert. It seems that we must accept the conclusion that there is much more 
uranium and thorium in these meteorites than appeared to be the case on the basis 
of the data reported for Modoc. 

The Heat Balance of the Earth.—The rates of heat generation by the radioactive 
elements have been calculated previously, but because of the improvement in data® 
it is possible to recalculate these values. By adding the energies of transitions, 
correcting the alpha-particle energies for recoil, and using 35 per cent of the maxi- 
mum beta-ray energies as the non-neutrino energy, values of the heat energy per 
gram for U2, U2?®, and Th**? have been calculated and are given in Table 1. There 


TABLE 1 
Heat Enercy Constants or K, U, anp Th 
HEAT Heat Eneray E . dE/dt 


ENERGY Eros/Gm (Cat/Gn Ergs/Gm See 
NUCLIDE Mev) x 10°") x 10°%) x 10°7 Cal/Gm Yr 


kK“ 071 0.171 0.410 0.935 0.223 
Th?*? 39.9 1.66 3.96 0 828 0.198 
[235 44.2 1.82 4.34 17 64 1 21 

[238 17.1 1.91 1.56 2. 935 0 701 


may be some doubt in regard to the choice of some data. In the case of Ac?® the 
energy was calculated as that of a beta ray of 1.7 mev plus a gamma ray of 0.80 
mev; for Pb?’ the energy was taken as a beta ray of 0.331 mev and a gamma ray 
of 0.0238 mev; for the Bi?!? to Pb transitions a total heat energy of 9.82 mev is 
used, which is the heat energy for the path Bi?!? —~ TI? — Pb®® and which is a 
reasonable estimate for the heat energy by the more involved Bi?!? + Po?!” 
path. The energies for other transitions seem to be obvious. The results are 
probably good within a few per cent. Bullard’ secured 43.7, 47.6, and 40.1 mev for 
U2, U5 and Th?*?, respectively, and these values are not greatly different from 
those used here. Birch'! secures almost identical values for U?* and Th**? but 


us 


—> Pb? 
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gives a l-mev larger value for U*®. The mean energy for K*® has been reviewed 
by Birch,'® who secures 0.708 mev. It seems difficult to improve on this value, 
even though recent evidence indicates a different branching ratio.'* The half- 
lives, A, are taken as K®, 1.27 & 10°; Th**?, 13.9 « 10°; U?®, 7.13 x 108; and 
Us 4.51 & 10°, all in years. The energy released per gram and per second given 
in the last two columns is then calculated from the equation dE /dt = E(0.693/X), 
where £ is the total heat energy per gram and dF ‘dt the rate of generation of heat 


energy. Then energies are given both in ergs per second and in calories per year. 


The abundances of the important radioactive elements in meteorites are not 
well known, but Paneth’s value for the Beddgelert meteorite is consistent with 
Patterson’s values, and for the present calculations we shall use 0.106 and 0.335 
p.p.m. for uranium and thorium, respectively, as the mean for the meteorites. 
Arguments have been given previously that the chondrites represent an average 
value for meteoritic abundances, and hence Davis’ values for the achondrites are 
not averaged with the Beddgelert chondrite. For potassium we shall use 880 
p.p.m., which is the average given by Edwards and Urey" and is closely the same 
as that of Ahrens, Pinson, and Kearns.’ The former authors have investigated 
the potassium content of the carbonaceous chondrites and found it somewhat less 
than that of the chondrites and more variable. These objects appear to be chon- 
dritic material of the high iron group’ which has been subjected to the action of 
water, and we do not know to what extent they should be included in a proper 
average. The potassium content of the achondrites is variable, and again it is 
difficult to estimate what their contribution to average material should be. We 
shall take Nier’s value of 0.000119 for the fraction of K“® in potassium at the 
present time. 

Table 2 gives values for the present instantaneous rate of heat generation and 
for the heat generated during the last 5 & 10°, 4.5 & 10°, 5.0 K 10°, and 5.5 & 10° 


TABLE 2 
Herat Propuction IN CHonpritic METEORITES 


dE/dt —Torat E 1~ Jovies rer GRaM or 
AT PRESEN1 CHonpritic MerTreortte iN Last 


Eroes/Yr Ga 5 xX 106 $5 « 10° 5 x 109 5.5 &® 10% 
NUCLIDE or METEORITE) Years Years Years Years 


kK“ 56 1,910 2,570 3,430 
Th? 2 139 1,400 1,570 750 
U3 0 9 1,090 1,770 2,890 
[238 3 162 2,020 2,370 2,690 
Total 7 366 6, 420 8, 280 760 

years in average chondritic meteoritic matter for each of the four important nu- 
clides K®, Th®*?, U?®, and U2*. The heat of fusion of such material is about 400 
joules per gram. If the earth’s mantle was made of this material and if it was at 
the melting point at the beginning, enough heat has been generated to melt nearly 
the entire mantle during the last 5 & 10° years. The rate of heat loss from the 
earth is estimated to be 9 XK 10” ergs per year," and this value is believed to be 
accurate within 20 per cent. In 5 X 10° years this would amount to 110 joules 
per gram, taking the mass of the earth’s mantle as 4.1 X 10” gm. Still, there 
should have been enough heat to melt some 65 per cent of the earth’s mantle. The 
large values for the heat generated during the last 4.5 & 10° years make it very 
probable that it would have been at\the melting point at the beginning of this 
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period, and in fact melting should have occurred early in the earth’s history, and 
the mantle should have remained at the melting point or above through all geologic 
time. It is impossible to believe that this has been true. We must conclude that 
the earth’s mantle does not contain as much of these radioactive «'ements as do the 
meteorites. Since the writer believes that serious difficulties wiil arise in under- 
standing the structures of the moon and of Mars unless the concentrations of these 
elements in the planets are less than in the meteorites, we must explore the possi- 
bility that abundances of these elements in the meteorites do not represent the 
cosmic abundances but that these elements have been concentrated in the meteor- 
ites. The alternative is that these elenients have been depleted in the earth, the 
moon, and Mars. This possibility will be explored later in this paper. 

At the Rome meeting of the U.G.G.L. (September, 1954) Dr. Roger Revelle sug- 
gested to me that the rate of heat loss of the earth is now and perhaps always has 
been equal to the rate of generation of heat in the earth. He suggested that a rise 
in temperature of the mantle will cause more rapid convections and a more rapid 
lissipation of heat, and a fall in temperature will produce the opposite effects until 
equilibrium is attained again. The assumption of near equality of heat loss and 
heat generation is made implicitly by Bullard” and by Birch,'! though they both 
believe that the earth is cooling and that the heat loss exceeds the rate of heat 
generation. If the temperature of the earth is falling with time, then the rate of 
loss will be higher than the instantaneous rate of heat production, and. if it is rising, 
the reverse is true. This statement of equality of rate of generation and loss of heat 
is accepted as true for the purposes of this paper, even though the writer favors a 
lower concentration of radioactive elements in the earth than is required by this 
principle. 

In accordance with Revelle’s suggestion, we conclude that the abundances of 
Kk, U, and Th in meteorites are too high for the earth. It is not necessary that 
these abundances shall be too high by exactly the same factor, but there is no evi- 
dent criterion for distinguishing between them, and hence the same factor will be 
applied. In order that the rate of heat generation shall equal the rate of loss, we 
must decrease the abundances of these elements by the ratio (7.00 & 4.1 & 10°) 
(9 X 10”), or 3.19. Applying this ratio to all the radioactive elements gives 276, 
0.0332, and 0.105 p.p.m. for K, U, and Th, respectively. In these calculations any 
contribution of heat from the earth’s core is neglected, since it is difficult to under- 
stand how these very electropositive elements could be present in any important 
amounts in a metailic iron core.6 However, these elements should follow the sili- 


cates in any separation of the iron and silicate fractions, and the larger amount of 


iron in solar material as compared with the earth’s mantle would produce a lower 
concentration of these elements in objects containing these larger amounts of 
elemental or oxidized iron. This could increase the factor to 4 instead of 3.19 for 
such objects. We shall use the more conservative value of 3.19, but remember that 
the values for heat production so calculated may be too large. Patterson con- 
cludes on the basis of his data that the factor should be 4, as compared to 3.19, and 
has noted this discrepancy between the rate of heat generation in meteorites and 
the rate of heat loss by the earth. The use of his data and factor would not change 
the results of this paper significantly. 

The abundance curve of the elements of odd mass number in the region between 
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the maxima at Mg” or Al” and at Fe’ appears to be irregular and neither smooth 
nor similar to the curve for the even mass elements in the same region, if the usual 
abundance of potassium, i.e., based on 0.088 per cent by weight, is used. How- 
ever, if the potassium abundance is lowered by a factor of 3.19, an irregular curve 
is still secured, but it is strikingly similar in its irregularity to the curve for the even 
mass elements. Such similarity of these curves may be expected according to 
the suggestions of Suess.'* Other elements, particularly Rb, Cs, Ba, Sr(?), and 
the rare earths, might behave similarly to potassium and hence be too abundant 
in the meteorites. The abundance curve in the neighborhood of these elements 
is not known sufficiently well to permit a reasonable decision as to whether lower 
abundances are indicated, but certainly there is no evident difficulty in selecting 
such lower abundances. Astronomical data seem to favor a higher abundance of 
potassium, but the variation in abundances for the elements in general in different 
astronomical sources or in the same source as studied by different investigators 
suggests that the revision suggested here is not outside probable errors of these 
measurements. !* 

The amount of potassium, uranium, and thorium assumed must not be less than 
that contained in the crust of the earth. Potassium constitutes about 2.6 per cent 
of the surface rocks of the earth, and, assuming 15-20 km. for the mean thickness 
of the éarth’s crust, we find that 50 per cent of the potassium is in the earth’s 
crust if its mean concentration in the mantle and crust is 276 p.p.m. Similar values 
are secured for uranium and thorium. Since these elements are more concentrated 
in granites, and since granites are presumably more characteristic of the surface 
rocks than of the deeper rocks of the crust, no serious contradiction exists between 
our assumed abundances and what we know of the earth’s composition. 

It is assumed that the concentrations of these elements in the chondritic meteor- 
ites are higher than the cosmic abundance, and this increased concentration must 
be explained by some reasonable process. These elements are concentrated at the 
earth’s surface by partial melting processes which are taking place at the present 
time. Potassium, uranium, and thorium are markedly concentrated in granites 


as compared to basalts and in both of these as compared to ultrabasic rocks. If 


these have been concentrated in the chondritic meteorites by processes similar to 
those occurring on earth, we must assume that melting processes such as those 
which have produced granites and basalts have been effective in the formation of 
some of the constituents of chondrites. Such constituents may well be the achon- 
dritic meteorites. These are of two principal kinds: those of approximately the 
composition of the silicates of chondrites, the chondritie (Urey and Craig) or the 
low-caleium (Prior) achondrites, and the basaltic (Urey and Craig) or the high- 
ealeium (Prior) achondrites. The former might be material from which some 
liquid has been separated (Johnstown is very low in sodium and p>tassium!* '), 
or they might be melted chondrites, and hence potassium would nut be low; this 
is true in some cases. Their uranium content is generally low, as would be ex- 
pected; this has been shown by Davis and Tilton’ in the case of Norton County, 
Cumberland, and Shaw. However, the second group, which might be the liquid 
fraction of the melting process, does not contain high concentrations of sodium and 
potassium.'* The uranium contents of this group have not been investigated. 
We do not know of any achondritic material which would be exactly the kind of 





132 CHEMISTRY: H. C. UREY Proc. N. ALS. 


material required to increase the concentration of potassium, uranium, and thorium 
in the chondritic meteorites. Further investigations are highly desired. Never- 
theless, we shall assume that such material existed and was incorporated in chon- 
dritic material. The over-all process of concentration was very inefficient com- 
pared with terrestrial processes, for an increase by a factor of only 3.19 is indicated 
for these elements in the meteorites, instead of about 100, as in the case of the 
earth’s crust. 

The Melting of Asteroidal-T ype Bodies.— Kuiper” has suggested that the asteroids 
and the moon have been melted by radioactive heating due to K, U, and Th, He 
assumes that these bodies melted internally, until only a thin outer layer of solid 
material remained, and then solidified by convection and cooled off by conduction, 
and that this left an outer layer with increased concentrations of these elements. 
The writer does not believe that this is a reasonable description of such a melting 
process. Chemical fractionation between phases, separation of the phases, and 
convection would occur during melting as well as during solidification. Kuiper 
assumes 2 completely quiet melting with a surface remaining intact. If the surface 
laver had the same composition as the melted liquid, it would surely sink into the 
liquid below. Moreover, we cannot assume that the liquid is more dense than the 
solid surface in some cases and not in others. If the liquid rises only to some dis- 
tance below the solid surface in the Hyginus Cleft and the Alpine Valley of the 
moon, as he assumes, this must mean that the liquid is more dense than the solid, as 
in the case of liquid and solid water. On the other hand, he supposes that liquid 
flows from the moon’s interior completely to the top of the Wargentin crater, 
which is a collision crater completely filled with what appears to have been liquid 
lava; this implies a less dense liquid than the solid surface of the moon, as is true of 
terrestrial lavas relative to the solid crust of the earth. Gilbert?! considered these 
same surface lunar features and concluded that great splashes of melted silicates 
produced by the collisions that formed the circular maria partially filled the Hyginus 
Cleft, completely filled the Wargentin crater, and produced other well-known 
flooded regions of the moon. It seems very doubtful that a moon with a com- 
pletely molten interior is necessary or at all possible as an explanation of any of 
the moon’s features. It is necessary to reconsider the process of melting of the 
moon and asteroidal bodies by radioactive heating. 

We assume that planetesimals accumulated from primordial dust by some proc- 
ess at present unknown.*? These objects were at some temperature at least low 
enough so that some water condensed, since it is retained both in the earth and in 


the meteorites. The iron was presumably in a completely oxidized condition, 


since elementary iron is unstable at low temperatures in a cosmic mixture of hydro- 
gen and water. They may have contained some ammonia and compounds of 
carbon. Since they accumulated at some time as grains of low thermal conduc- 
tivity, radioactive heat would slowly raise their temperatures. 

In general, melted terrestrial silicates are less dense than the solid silicates with 
which they are in chemical equilibrium. Wahl,** however, has argued that this 
condition may not be true. Bowen and Schairer*‘ have studied the MgO-FeO-SiO, 
system, and their results show that a liquid may have a higher density than the 
solid with which it is in equilibrium, provided that the FeO content is high. Such 
high concentrations of FeO do not occur generally in the earth, so far as we know, 
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but they might well occur generally in the moon, which, because of its low density, 
matter in which the iron is completely 
or nearly completely oxidized to FeO and FeS. The abundance value for sulfur 
as secured from the meteorites is probably too low because of the escape of hydrogen 
sulfide during their formation. Suess estimates its atomic abundance as 2,510, 
using silicon equal to 10,000. Iron combines with sulfur rather than with oxygen, 
and thus sulfur removes iron from the silicate fraction, since the iron sulfide and 
silicate are only slightly soluble in each other. By using Urey and Craig’s'* low- 
iron-group abundances, which they show give material of about lunar composition 
as judged from densities, and using Suess’s value for the sulfur abundance, per- 
centages of minerals of the cosmic silicate fraction with the iron completely oxidized 
to FeS and FeO have been calculated and are given in Table 3. This neglects all 
minor constituents. The weight fraction of FeO in the silicate is 25 per cent. 
The calculations could be modified in the direction of lower percentages of iron 
oxide and of the olivines by assuming higher sulfur and lower iron abundances. 
(The writer believes that a somewhat lower abundance of iron is probable.) It 
seems certain that liquid high in the feldspar constituents, with some water content 
and thus of low density, would be produced first in a melting process. Such a 
liquid would rise in a gravitational field if sufficient time were available for the 
separation of the phases. This liquid would carry the radioactive elements K, U, 


’ 


probably consists of cosmic ‘nonvolatile’ 


TABLE 3 
Per Cent Per Cent 
NasO+ ALO,;- 6810 9.7 FeO-SiO, 6.7 
CaO+ ALO - 28102 a4 Mg0-Si0, 10.3 
CaQ- SiO, 3.9 2FeO-SiO, 29.3 
2Mg0-SiO, 37.4 


and Th and other elements not important to the present discussion from the in- 
terior toward the surface, just as has occurred in the general surface regions of the 
earth. This happens because the potassium goes with the feldspar fraction, and 
uranium and thorium are insoluble in the solid silicates. The densities would vary 
with the ferrous iron content, but the solubilities for these elements would not be 
changed markedly by the substitution of ferrous iron for magnesium ion in the 
melted silicates. The removal of these elements from the interior by the partial 
melting process would hardly be complete. Thus heating in the interior would be 
reduced but not completely eliminated. After this process had occurred, it might 
well be that any liquid produced would be as dense as, or more dense than, the 
solids because of the concentration of the iron in the liquid, and little or no sepa- 
ration of liquid and solid would occur. If the amounts of the radioactive elements 
characteristic of chondritie meteorites had been originally present, enough of them 
would probably have remained in the interior to produce complete melting, and, if 
this had been the case, the deep interior of the moon would be melted now. If the 
first liquid to melt was more dense than the solid, it should have transported the 
radioactive elements toward the interior. Also, if this had been the case, the sur- 
face regions of a completely molten moon would have solidified first and the radio- 
active elements would have concentrated in the residue and would again have been 
concentrated in the interior. These conditions would have resulted in a molten 
interior of the moon throughout all geologic time and are inconsistent with the 
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nonisostatic shape of the moon. For the purposes of the following discussion, we 
assume that the first liquid to melt will be less dense than the solid phases in con- 
tact with it. 

The character of the heating process within cosmic objects of varying sizes de- 
pends on many factors. If an object were small, it would not heat up at all; if 
somewhat larger, it would dry out by the diffusion of water vapor to the surface; 
if still larger, liquid water with some substances dissolved in it might flow slowly 
to regions near its surface and evaporate; if larger yet, or in longer times, the lower 
melting silicates containing some water would flow toward the surface, break 
through, and spread on the surface and solidify. These latter melts would carry 
the radioactive elements toward the surface, and thus the rate of heating in the 
interior would decrease, and it would increase at regions to which the melt moved. 
The basalts of the earth are produced by some such process, and lava flows are 
. in temperature. However, if the objects contained little 


never above 1,200° ¢ 
water, the temperatures might become higher, and, if water was present, they 
could be considerably lower than 1,200° C. The processes just outlined for the 
largest objects might apply to the asteroids and the moon. It is not at all certain 
that the object would ever melt completely. Even the surface regions would 
probably never be melted except locally and temporarily, and subsurface 
temperatures would remain somewhere near the melting point of hydrous silicates. 

Using the abundances of K, U, and Th given above, we can calculate the amounts 
of heat which would be generated by radioactive substances. Since the date for 
the formation of the solar system seems to be some 4.5 & 10° years ago,”: * and 
some time is needed for radioactive heating, we shall calculate these heat energies 
for the last 4.5 & 10° and 5.5 & 10° vears and for the interval between 5.5 and 4.5 
billion years ago. The results are given in Table 4. The heat of fusion of silicates 

TABLE 4 


Hear EnerGy in Joutes per Gram GENERATED IN INDICATED TIMES 
(K = 276 p.p.m., U = 0.0332 p.p.m., Th = 0.105 p.p.m.) 


4.5 & 10° Years 5.5 * 10° Years 5.5 & 10° Years to 
Nuclide to Present to Present 15 & 10° Years 


K' 600 180 
Th232 440 5 110 
U2 340 570 
L238 630 210 
Total 2,010 : 1,370 
of cosmic composition is about 400 joules per gram, and the heat capacity is about 
1 joule per gram at low temperatures. The heat of fusion of water is 333 joules 
and of vaporization 2,260 joules per gram; hence any water in these objects would 
require considerable heat for its melting and evaporation. Some 2,000 joules of 
heat are required to raise the temperature of silicates from ordinary temperatures to 


their melting points and to melt them even if no free water is present. The 1,370 
joules of total heat per gram given in the last column of Table 4 is not sufficient, 
and yet some 10° years are required for its production. It is difficult to meet the 


requirements of no net heating or cooling of the earth at the present time and at 
the same time melt the moon or asteroids by radioactive heat. 

The Chondritic M eteorites.—Wasserburg and Hayden* have determined the ages 
of the chondritic meteorites by the K“-A* method and find 4.7 & 10° and 4.8 X 
10° years; and Patterson ef al.,% from the ratio of Pb®* to Pb*” in the chondrites, 
assuming that the lead of Canyon Diablo and Henbury iron meteorites gives the 
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primitive isotopic ratio, find 4.5 K 10° years. Since the half-lives of U?® and U?* 
are more reliably known than the half-life and branching ratio of K®, Patterson’s 
date is probably the more reliable. The K“-A“® dates are particularly informa- 
tive. They mean that A“ has not been lost from the materials of these meteor- 
ites for this length of time, and this means that all melting processes which produced 
the minerals of these objects occurred more than 4.5 X 10° years ago. The chon- 
drules of these objects are best explained as having been produced as molten drop- 
lets which have only partly crystallized even at the present time. The process 
which broke up the crystals and produced the melted chondrules was one of con- 
siderable violence and would result in the loss of A® from these objects. It seems 
most likely that this process also occurred before or at 4.5-4.8 X 10° years ago. 
The fragments could not have accumulated much later by falling on an asteroidal- 
type body, because small objects some millimeters and smaller, such as the chon- 
drules, would have been moved toward the sun, due to the Poynting-Robertson 
effect, and would no longer remain in the neighborhood of the asteroids after the 
sun became a bright star.” The chondritic meteorites contain bits of plessite, 
kamacite, and taenite, the iron-nickel alloys of the Widmanstatten figures of the 
iron meteorites. Hence this differentiation had taken place before these dates. 
This process of producing the chondrites required as a minimum (1) the reduction 
of oxidized iron to metal, the melting of the iron, its accumulation into massive 
masses, and its slow cooling to form the iron-nickel alloys; (2) the melting of sili- 
cates and their crystallization to produce the minerals of the chondrites and achon- 
drites; (3) the violent brecciation of these materials; and (4) the accumulation of 
these materials into the chondrites. The processes were partially repeated, since 
pieces of chondrites and achondrites of particular distinctive types are often im- 
bedded in stone meteorites of other types.** No ages of these polymict meteorites 
have been measured. The main process of meteorite formation was evidently 
completed 4.5 X 10° or more years ago. It is difficult to account by radioactive 
heating for the melting which produced the materials of the meteorites and at the 
same time to account for the heat balance of the earth at the present time. A 
somewhat larger amount of radioactive elements or an age for asteroidal-type ob- 
jects even greater than 5.5 X 10° years would be required. One wonders where 
these objects were and what evolutionary processes of the solar system were occur- 
ring during the long time between the accumulation of such objects and the period 
of violent collisions which formed the chondrites some 4.5 10° years ago. 

Yet we have assumed that the radioactive elements have been concentrated in 
the meteorites by a factor of 3.19 by a melting process. We shall try to outline a 
possible radioactive heating process by which this could have occurred in spite of 
the obvious difficulties noted above. As the temperature of an asteroidal body 
slowly increased, consolidation of its material would occur and partial melting of 
its lower melting constituents would take place. The radioactive elements K, U, 
and Th would dissolve in the first liquid, and this liquid would flow toward the sur- 
face. This would decrease the heating on the interior and increase that toward the 
surface. The exact events would defy any confident quantitative description; 
they depend on the strength of the gravitational field, the water content (which 
affects the viscosity and melting points), and the thermal conductivity of the object. 
This melting process is apparently occurring on the earth at the present time, and 


even in this case no quantitative description has been possible. However, the 
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differentiation of the surface materials required in order to produce increased con- 
centrations of the radioactive elements in the meteorites does not appear to be 
impossible. Some iron may even have melted in regions where the radioactive ele- 
ments were concentrated. If the massive metallic iron-nickel was produced in 
this way, then reducing agents such as carbon compounds must have been present 
to reduce the oxidized iron present in the primitive materials. Such reduction 
processes would consume heat, as would the vaporization of volatiles generally. 
Also, the different chemical composition of the high and low iron groups" and a 
mechanism for the intense brecciation of the chondrites must be explained. It 
seems most likely that this chemical fractionation is related to the general frac- 
tionation between the silicate and metallic iron fractions required to explain the 
varying densities of the planets.” 

The carbonaceous chondrites contain appreciable quantities of water and other 
volatiles, and chemical analyses which will be published shortly by H. Wiik show 
that these objects belong in the Urey and Craig high iron group."* X-ray diffraction 
studies made by Miss Ann Pletinger and the writer show that these objects are 
probably ordinary chondrites which have been subjected to water solutions. The 
Orgueil and Ivuna meteorites, for example, show the X-ray diffraction patterns 
of magnetite with only very faint indications of any other pattern, and magnetite 
is the stable oxide produced by the action of water on metallic iron. Other car- 
bonaceous chondrites indicate that hydrated silicates are present in addition. 
Though we may be able to account for the increased concentrations of the radio- 
active elements, and possibly for the fusion of iron and silicates by radioactive 


heating due to K, U, and Th over a long period of time, many difficult points re- 
main, and other possibilities should be considered. 

However, one point must be made definite. The discussion of this paper is con- 
sistent only if the chondritic meteorites originated from the surface regions of an 
asteroidal body. This requires revision of previously expressed views of many 


people, including those of the writer. 

The Melting of the Moon.—The moon has differences in moments of inertia cal- 
culated from its dynamical motions which indicate a projection toward and away 
from the earth about 1 km. in height, as compared to the polar radius, whereas the 
equilibrium height should be only about one-seventeenth of this.* The writer has 
maintained that this indicates great strength in the moon, that it is not now molten 
on the interior, and that, if it had ever been molten, even moderate amounts of 
radioactive heating would have maintained its temperature near the melting points 
of the silicates. This calculation has been repeated using more recent data, and 
the previous conclusion seems to be justified. The writer suggested that the 
irregular shape might be due to the fall of the large planetesimals on the surface 
of the moon. Regardless of the time and mode of formation of this projection, it 
must be true that the moon has had a very rigid structure since it was formed. By 
contrast, Kuiper® concludes that this projection is not due to a frozen tidal wave 
but that it was formed by the addition of mass to localized areas by the objects 
which formed the maria and that these events occurred at or near the time when 
the moon was totally or largely in a melted state. Kuiper’s assumptions appear 
to be inconsistent. Surely the moon would take an equilibrium shape when com- 
pletely melted, regardless of whether objects fell on it or not. This nonequilibrium 
shape of the moon is an important datum in considering its thermal history. 
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The writer has regarded the moon as a primitive object like the planetesimals, 
which he found necessary for the understanding of the chemical problems relating 
to the origin of the solar system. The primary argument for this view depends 
on the moon’s density. As is well known, this density is 3.34, or about 3.31 cor- 
rected for pressure. The mean temperature of the moon can hardly be 300° K. 
If it is 1,200° C., the density of its material would be about 3.45 at room temper- 
ature and zero pressure, and this is about the density expected for cosmic nonvolatile 
material if all the iron has been oxidized to FeO and FeS."* The figures can be 
manipulated to some extent by assuming the presence of some water in the moon’s 
substance, which would lower the density, and of some metallic iron, which would 
raise its density. This composition of the moon requires that it have been sepa- 
rated from gases of cosmic composition at a low temperature, for otherwise the iron 
would be substantially completely reduced to the metallic state.** 

But some radioactive heating of the moon must have occurred. | This has been 
considered previously by the writer.** The abundance of potassium in the meteor- 
ites was taken as 0.2 per cent, whereas it is now considered to be 0.088 per cent. 
The abundance of uranium was taken as 1.3 X 10~° gm/gm, while now it appears 
to be about eight times as great; and the age was taken as 3 X 10° years, while now 
it is taken as 4.5 X 10° years. Needless to say, the previous calculations need 
revision. A simple substitution of the new values for the elements K, U, and Th in 
chondritic meteorites leads to a highly melted moon and earth, as explained before. 
However, the lower abundances assumed in this paper lead to an acceptable con- 
clusion. A heating process such as that described for an asteroidal body can be 
applied qualitatively without change to the moon. Due to greater size and gravi- 
tational field, the process would be quantitatively different. The increased field 
would favor the separation of liquid and solid silicates, but viscosities might be 
increased as well by higher pressures. The transfer of radioactive elements to 
tie surface or toward the surface would enable the radioactive heat to be lost more 
readily than was previously caleulated by the writer. It is possible that great ri- 
gidity would be preserved in the moon since the time of formation of the lunar pro- 
jection toward and away from the earth. A little consideration of the data of 
Table 4 shows that it is only necessary that about one-half of the original concen- 
tration of kK, U, and Th should have been transported to the outer parts of the 
moon, from which the heat is readily lost by conduction, in order that the tem- 


perature of the moon should stay below the melting point of its silicates, if the 
moon accumulated some 5.5 X 10° years ago. The figures can be manipulated 
in various ways, again, because of uncertainties regarding initial temperatures, 
thermal conductivity and the consequent loss of heat, amount of heat generated 
as determined by amounts of radioactive materials and age, ete.; but one sees that 
the nonequilibrium shape may well have been formed some 4.5 X 10° years ago 


and still remain today on the moon. 

But, if heating processes other than radioactive heating by U, Th, and K (see 
below) produced the melting in the meteorites, then the moon and asteroids may 
have formed some 4.5 & 10° instead of 5.5 & 10° years ago. In this case a slightly 
smaller concentration of these elements in the moon than is assumed in this paper 
would mean that no large transfer of these elements to the surface is required to 
avoid melting of the moon, provided that the moon accumulated at low tempera- 


tures. 





138 CHEMISTRY: H.C. UREY Proc. N. A. 8. 


Heat Balance of the arth and Mars.—The earth is a derived object and not a 
primitive one. This is shown by the essential absence of the high-atomic-weight 
inert gases. These gases must have escaped from the earth’s protoplanet before 
the earth had accumulated and while its mass was distributed over a large volume. 
The present material of the earth must have been in solid or liquid form, for other- 
wise the inert gases could not have been separated from the volatilized “nonvolatile” 
fraction, and it must have been at low temperatures and hence solid, for otherwise 
the more volatile constituents of the “nonvolatile” fraction would have escaped,** 
which is not the case. The writer has described this stage as a miniature globular- 
cluster type of object. When the inert gases and most of the hydrogen, water, 
methane, and ammonia were gone, only then could the earth have accumulated. 
It is a derived object, and younger than the moon. (The writer would not try to 
state how much younger.) The total heat generated in the last 4.5 billion years 
was thus 2,000 joules per gram, and at the present rate of loss of heat the earth 
would have lost 1,000 joules of heat per gram. It seems not unreasonable that 
the rate of loss of heat may have been greater in the past. Hence the earth can quite 
easily have remained at temperatures comparable to those of the present time 
even if the energy of the formation of the core is included. In fact, it is difficult 
to account for the high values of the internal temperatures generally recognized 
as existing within the earth,** unless the energy of the formation of the core is in- 
cluded if it accumulated at low temperature, as the writer has argued on evi- 
dence that has not been criticized.* The energy of core foymation has been estimated 
as 1.67 X 10* ergs, with an error which is difficult to estimate. The radioactive 
heating during 4.5 & 10° vears is calculated to be 2 X 10" ergs per gram of the 
mantle, or a total of 8.2 & 10” ergs. At the present rate of heat loss of 9 X 10” 
ergs per year a total loss of 4.05 & 10” has occurred. This gives a total of 20.85 X 
10* ergs, or 3,500 joules per gram of the earth. If the loss of heat was larger in the 
past than it is now, which seems probable, no insuperable difficulties in regard to 
the heat balance of the earth exist. In fact, according to Revelle, the earth has 
always been able to dissipate its heat approximately as rapidly as it is generated. 
Proponents of an initially molten earth will secure this additional heat from the 
“original” heat rather than from the energy of formation of the core. 

The Clairaut constant for Mars is 1.22, using Trumpler’s value for its radius, or 
1.19 using Kuiper’s value.” These values are so close to the value for a planet of 
uniform density, namely, 1.25, that they indicate that even the difference between 
them is due to increased densities in the interior produced by compression only. 
The density of Mars is 4.18, which gives 3.96 estimated for zero pressure. This 
indicates that the planet contains about 30 per cent by weight of metallic iron. 
These conclusions can be changed by a redetermination of the radius of Mars, but 
for the present the conclusion in regard to Mars is that, though it contains approxi- 
mately 30 per cent of metallic iron, it has no core or only a small one. This con- 
clusion is inconsistent with the high rate of generation of heat calculated in Table 
2, for the planet would have melted and formed a core regardless of its temper- 
ature of formation. 

However, the 2,000 joules per gram of cosmic matter generated in 4.5 & 10° 


vears (Table 4) is a manageable figure, at least. Some loss of heat has occurred. 


Though there seem to be no mountains on Mars, lava flows may have occurred. 
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Solid iron would hardly sink through solid silicates, and at ordinary pressures 
silicates melt below the melting point of iron-nickel of meteoritic composition and, 
in fact, if some water is present, some 500° below the melting point of iron. Hence 
the partial transfer of radioactive elements from the interior of Mars to its surface 
may have occurred. McLaughlin® has suggested that volcanoes do occur on Mars. 
With a moderate temperature of formation some 4.5 X 10° years ago, the low 
concentrations of the radioactive elements assumed here, and their possible transfer 
to the surface to a reasonable degree, the low temperature of the planet as indi- 
cated by its lack of a core can be understood. 

Effect of Other Permissible Abundances of K, U, and Th.—Up to this point it has 
been assumed that the particular choice of abundances for K, U, and Th, namely, 
276, 0.0332, and 0.105 p.p.m., respectively, is rigorously correct. This, of course, 
is only a working hypothesis. The total heat generated may be larger if the earth 
is heating up or lower if the earth is cooling down. Because of the marked chemi- 
cal similarity of UO. and ThOs, these two elements should accompany each other 
and should be concentrated in the chondrites to the same degree. On the other 
hand, potassium may not be equally concentrated in the chondrites as compared 
to uranium and thorium. Thus it is possible that the cosmic abundance of po- 
tassium is higher or lower and that of the uranium and thorium lower or higher, 
respectively, than has been assumed. If the potassium abundance is higher than 
has been assumed, somewhat more heat would be generated during geologic time and 
somewhat more heating would be available in earlier times, and if it is lower some- 
what less heating would be available. Such amounts of heat may be consistent 
with the present high rigidity of the moon and the absence of a core in Mars. As 
stated previously, it seems more likely to the writer that the total amounts of these 
elements are, in fact, lower than has been assumed in this paper. In this paper an 
attempt is made to show that these data on the moon and Mars and the radioactive 
heat required by the assumption of equality of the rate of loss and generation 
of heat in the earth can be consistent with each other. 

As was pointed out previously, an alternative assumption would be that the 
earth, the moon, and Mars had lost part of their radioactive elements during their 
formation, while the meteorites retained them in more nearly their proper cosmic 
abundances. Recently Dr. H. E. Suess outlined a possible course of events to the 
writer which might have accomplished this and which is presented here with his 
permission. 

After the elements were synthesized at high temperatures, the mass cooled down 
and condensation of the elements occurred, with considerable segregation of the 
elements in different solid particles such as metallic iron, the silicates, and oxides. 
These crystals had different densities and, during the formation of the planets, be- 
came sorted mechanically, with variable amounts of the different chemical sub- 
stances being retained by the different planets, thus producing the different densi- 


ties of the planets. In this process the radioactive elements were lost preferentially 
relative to the silicate fraction, while the silicate fraction itself was being lost rela- 


tive to the metallic fraction. Von Weizsicker concludes that an extended system 
such as a protoplanet having high angular momentum would separate into a cen- 
tral body, the planet, and a fraction which is lost to space and which carries the 
excess angular momentum. However, the details of the process are too complex 
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for exact mathematical and physical description. Such a process might explain 
Mars and the earth, but some special assumptions must be made for the moon, 
which did not lose much of its silicate fraction as compared to its metallic frac- 
tion, that is, if it has any iron fraction at all. 

The suggestion does not seem to be an impossible one at all in view of the uncer- 
tainty of the process of planet formation. However, it is difficult to say what the 
physical differences between the fraction which carried the radioactive elements 
and the other fractions were, since these elements are so definitely oxyphile in 
character and are, in fact, among the most electropositive elements. These ele- 
ments should have condensed as part of some silicate fraction, and this fraction 
could hardly have had a density very different from that of any other silicate 
fraction. What is required is a separation of the one fraction from the other, since 
these elements—potassium, uranium and thorium—will not be present in the 
metallic phase. Perhaps some difference in particle size may have made a sepa- 
ration possible. 

It is still necessary to secure an accumulation of the bodies which produced the 
meteorites and to account for the melting processes which produced the silicate 
minerals and the massive iron of the meteorites. Also, we must account for the 
massive objects of high tensile strength which plowed out the grooves in the center 
of the moon’s disk.** 

If larger amounts of the radioactive elements were present in the primordial 
materials, greater heating effects could be produced. Yet the data of Table 2 show 
that some 5 X 10° years, ie., from 5 & 10° to 4.5 & 10° years ago, would be re- 


quired to produce some 1,860 joules of heat per gram, which is somewhat less than 


the amount required to raise meteoritic material from low temperatures to the 
melting point and meit it, but sufficient to partially melt such material in com- 
pletely insulated bodies. This is a long period of time, and one wonders at the 
required series of events, i.e., an accumulation of bodies large enough to be ther- 
mally insulated, a quiet period of 5 & 10° years for radioactive heating, a violent 
brecciation process to produce chondrites, and then no processes sufficiently violent 
to remove A“ for the next 4.5 X 10° years! 

A Reconsideration of Radioactive and Other Sources of Heat.—The arguments of 
this paper have been presented as favorably to the radioactive heating postulate as 
possible, but it is very difficult to believe that heating by K, U, and Th is a feasible 
explanation for the high-temperature stage required to produce the meteorites and 
to produce the fractionation of the silicate and metallic phases in the planets gen- 
erally.*° Urey® has considered three sources for these high temperatures, the 
required chemical separations, and the collisional processes for the formation of the 
meteorites: (1) Evaporation of silicates by a high-temperature sun early in its 
history was suggested. There appears to be no evidence for such behavior of a 
new star, but we do not know everything about the process of forming stars. (2) 
An attempt was made to modify the Eucken-Kuiper condensation process in pro- 
toplanets in such a way as to avoid fractional loss of the more volatile elements. 
It is difficult to secure high temperatures by the adiabatic compression of gases 
and at the same time to lose those gases, but the process has not been described 
quantitatively because of its great complexity. (3) It was postulated that the 
planetesimals of the asteroidal type antedated the formation of the sun and that 








a 

















Vou. 41, 1955 CHEMISTRY: H.C. UREY 141 


the solar nebula was formed during a collision of the new sun with a gas, dust, and 
planetesimal cloud, with the heating and collisional processes occurring during this 
collision. It seems possible chat any of these heating processes may have produced 
an increase in the K, U, and Th concentrations on the surfaces of the asteroids by a 
factor of 3.19 or more, as required by the present calculations, though it is easier to 
describe plausible models if radioactive heating is assumed. 

The age of the solar system now appears to be so great that we may question 
whether the accretion of substantial objects did not occur so soon after the syn- 
thesis of the elements that some of the radioelements with half-lives of some million 
years may not have completely decayed when this occurred. Blaauw*! has shown 
that the stars in the Perseus cluster must have originated 1.5 X 10® years ago and 
hence that some very rapid mechanism for the formation of stars exists. Without 


a model for the formation of these planetesimals, it seems impossible to conclude 
that they could not also have formed in such a period of time. There are eight 
known radioactive nuclides that may be important as sources of heat within 10% 
10’ years after the synthesis of the elements. These are listed in Table 5, together 


TABLE 5 
Heat ENERGIES OF SOME RADIOACTIVE NUCLIDES 
Heat Abundance of Energy Energy 


Energy Half-Life Neighboring Nuclide (Ergs/Gm (Joules/Gm 
Nuclide (Mev) Years) (P.p.m.) of Nuclide) of Meteorite) 


Als ° 6.3 10° ‘5 Mg** 7500 1.2 X10" 9 xX 10° 
Cl 0.285 x10 Cl# 120 6 «K 10" 8 xk 104 
Zr 0.021 5 xX 10°* 1/,Nb* 0.2 2 x 104 
Pd 0.035 ¢ 108 1/, Agi? 0.67 l « 10'4 
['=9 0.052 1 x 10° J !37* 1.25 6 x< 10" 

0.038 


0 090 
C's!3 0.070 2.9 xX 106 Ba'® 0.19 5 «x 10'4 10 
[236 1 58 2.39 X 10 [235+ 0 0666 1 87 x 10" 125 
Np*” 44.01 Z.a KW (2/3)U2f 0.0222 1 x 10" 394 
* The abundance of I'* is probably too high by a factor of 5 


+ One-third of the abundance of U?" 4.5 & 10% years ago is used as an estimate of the primordial abundance of 
the odd mass-odd atomic number element Np?”, and that of U** as equal to that of U®* at the same time. 


with their heat energies of disintegration (35 per cent of the maximum beta-ray 
energies are used again), their half-lives, the abundance of a similar stable element 
which is taken as an estimate of their primordial abundance, and the energy per 
gram of nuclide and per gram of meteorite produced by their complete disintegra- 
tion. The calculations are very approximate, but it seems doubtful whether 
any nuclide except Al’ and Cl*, or some other at present unknown nu- 
clides, need be considered, and Cl*®* has such a short half-life that its effect 
would be negligible as compared to Al*. Hence only this recently discovered 
nuclide Al** need be considered at all.*? Even after the lapse of some 15 half- 
lives, the amount of remaining heat from the Al** would be 3,000 joules per gram. 
This element also is concentrated at the earth’s surface by partial fusion processes. 
If this heat was available, a great variety of assumptions can be made, since the 
amount of heat depends so markedly on the time of formation of insulated bodies. 
The possibility of radioactive heating as a source of these heating processes has 
been considered by the writer previously, but heating by such means has seemed 
to him inconsistent with other facts. Uranium, thorium, and potassium either re- 
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quire such long periods of time to produce the early heating processes, i.e., some 
10° years, that they seemed highly improbable, or their abundances must be so 
high that the irregular shape of the moon and the absence of a core in Mars could 
not be understood. Radioactive elements in the mass range above iron could not 
furnish a sufficient amount of heat to be important, as is also true of U2" and Np*, 
unless other, as yet unknown radioactive nuclides exist. Cl has such a short 
half-life that very rapid accumulation of the planetesimals would be necessary, and 
this seemed improbable. However, Al** perhaps gives a real possibility of ex- 
plaining these heating processes. Its energy is large, and perhaps accumulation of 
planetesimals would be possible within some 15 X 10® years, that is, not too quickly 
or too slowly, and hence some proper amount of heating is possible. 

The order of events becomes the synthesis of the elements during some major 
process such as that which formed the Perseus cluster of stars, the accumulation of 
the sun and planetesimals of the asteroidal and lunar type within some 10’ years, 
heating due to Al**, a fractionation process which partially separated the silicate 
and metallic fractions and separated these in varying amounts in different planets, 
and accumulation of the planets during some fairly long time. These processes 
included some intense collisional mechanism during which the brecciated structure 
of the meteorites was formed. All these processes may well have occurred some 
1.5-5.0 X 10° years ago, during a time span of some 10’ years. 

It seems most probable that some combination of these various processes occurred 
during the formation of the solar system. 

Conclusions. Reasons have been given for believing that the abundances of the 
elements K, U, and Th in the chondritic meteorites are higher than the mean 
abundances of these elements in planetary bodies by about a factor of 3.19. Using 
decreased abundances of these elements, it is shown that the thermal properties of 
the earth, the permanent nonequilibrium shape of the moon, and the absence of a 
core in the planet Mars can be understood. The lower abundance of potassium 
gives a more regular curve for the abundance of elements in this region of the 
periodic system, though it disagrees with the astronomical data in this respect. 
[t is difficult to secure sufficient heat for this purpose on the basis of abundances of 
K, U, and Th as used in this paper. It is suggested that Al?* may have supplied 
the radioactive heat for this purpose and that other heating processes previously 
suggested may have been important as well. As an alternative, it is also pessible 
that depletion of these radioactive elements in the earth, the moon, and Mars dur- 
ing their formation may account for the lower abundances of these elements re- 
quired to secure reasonable thermal histories for these planets. 


The writer is indebted to Professor Roger Revelle for his statement of the con- 
clusion that the rate of heat less from the earth has always been approximately 
equal to the rate of heat generation. This idea gave the key to a solution of prob- 
lems that have troubled me for several years. Thanks are also due to Dr. C. 
Patterson for a discussion of his analytical results on lead and uranium in meteor- 
ites before they were published, to Professor Julian Goldschmidt and Dr. Harmon 
Craig for their discussion of the melting processes in silicates, to Dr. Maria Goeppert 
Mayer for a very constructive criticism of this paper, and to Dr. H. E. Suess for 
his presentation of alternative possibilities noted in the text. Apologies are due 
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Professor F, A. Paneth and his co-workers for the writer’s skepticism in regard to 
their determinations of the uranium and thorium contents of meteorites for these 
several years. 


* This work was partially supported under contract AT(11-1)-101 with the U.S. Atomic 
Energy Commission. 
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SOLUTION OF A PROCESS OF RANDOM GENETIC DRIFT WITH A 
CONTINUOUS MODEL* 


By Moroo Kimura 
DEPARTMENT OF GENETICS, T UNIVERSITY OF WISCONSIN 


Communicated bu Sewall Wr ight, November 15, 1954 


The problem of random genetic drift in finite populations due to random sampling 
of gametes in reproduction was first treated mathematically by R. A. Fisher,! 
using a differential equation. Fisher’s general method was appropriate, but owing 
to the omission of a term in the equation, his result for the rate of decay of variance 
was only half large enough. The correct solution for the state of steady decay was 
first supplied by 8S. Wright,? using the method of path coefficients and an integral 
equation. 

Later Fisher® corrected his results and also elaborated the terminal part of the 
distribution in the statistical equilibrium by his method of functional equations, 

In all these works, however, it was assumed that a state of steady decay had been 
attained, but nothing was known about the complete solution which might show how 
the process finally leads to the state of steady decay. The present writer, by cal- 
culating the moments of the distribution and with the help of the Fokker-Planck 


equation, obtained a solution which assumed an infinite series under the continuous 
model, showing that the process approaches asymptotically the state of steady 
decay. At that time, however, only the first few coefficients in the terms of the 
series could be determined. Pursuing the problem further, he arrived at the com- 
plete solution, which will be reported here. After obtaining these results, the write: 
recently discovered the work of 8S. Goldberg.’ In his unpublished thesis Goldberg 
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solved the diffusion equation for the gene-frequency distribution in a finite popu- 
lation when recurrent mutations occur. His solutions have a direct connection 
with the frequency distribution of unfixed classes in the case of pure random 
genetic drift, as will be seen below, 

Consider a random mating population of N diploid parents. Let A and A’ bea 
pair of alleles with frequencies z and | — z, respectively. In order to single out the 
effect of random drift, we shall assume an idealized situation in which selection, 
migration, and mutation are absent and generations do not overlap. The process 
of the change in gene frequencies is most adequately described by giving the fre- 
quency distribution f(z, t) at the tth generation, where z takes on a series of dis- 
crete values: 0,1/2N,2/2N,...,1—1/2N,1. For fairly large V, however, x can 
be treated as a continuous variable without serious error. 

First, we shall derive the moment formula which is useful heuristically. Let 
x, be the gene frequency in the (th generation, and let 6a, be the amount of change 
due to random sampling of gametes per generation, such that 


Dena = 2, + G2; (1) 


rt +1) 


Let pu’), = K(x...) be the nth moment of the distribution about 0 in the 
(¢ + 1)th generation. We write the expectation of x7,; in terms of (2, + 6z,)" in 
two steps: (1) taking expectation for the random change 6z,, which will be de- 
noted by £;, and (2) taking the expectation for the existing distribution, which will 
be denoted by Fg. 

Noting that £,(éz,) = 0, F(6x,)? = xr,(1 — 2,)/2N, ete., 


pit) = E(2, + &,)* 
Eyi2i + (pzit "BY, (8x) + (329 ~7E, (52,2? +... 


ae nin ~ 1) 9 %Al — 2;) 1 \i 
= E, ) om 2 Mom ON +0 rare 


The intrinsic assumption in the continuous model is that the effective siz V is 
sufficiently large so that terms of order 1/N? and higher can be omitted without 
serious error. 
Thus 
nin —1)) wa n(n — 1) 
ie, ee = er 4N 


r(t) 


Mn-l- 


For large NV, the moments change very slowly per generation, and we can replace 
the above equation by the system of differential equations: 


#(t) ql 
dy epee n(n pie 1) f v(t) pi i} (n = 1 9 3 
: q Oy cae 


dt an. t#e (3) 


If the population starts from the gene frequency p (0 < p < 1), u’\” = p", and we 
can obtain the nth moment as the solution of (3), 


w'? = pt i + 1)pq(-1)'F(1 — 3,2 + 2, 2, p) X 
1=] 


(n— 1)...(n —12) 
(n+ 1)...(n + 2) 


—[i(i+1)/4N]t 
’ 
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where F(1 — 7, 4 + 2, 2, p) is the hypergeometric function and q = 1 — p. For 
finite n the series is finite. Putting n = 1, 2, 3, and 4, it will be seen that, for large 
N, the resulting formulas give a very good approximation to the exact moment 
formulas obtained by A. Robertson® (p. 205). 
The probability f(1, 4) of the gene A being fixed in the population by the éth 
generation can be obtained by using the relation 
fai,d) = lim > 2f(z,t) = lim p’™. 


—— 
n—-> © z= 0 n> «a 


The resulting series, 


fii +1)/4N]t 


fa,0 =p+ ¥ (i+ ipq(—1)'F(1 — 7, i + 2, 2, pe! , 3 
1 I 


is now an infinite series whose convergence must be examined. It is convenient 
. . . om . ° 
here to introduce the Gegenbauer polynomial 7° ;(z) which is related to the hyper- 


geometric function by 


T) (2) = ua + I) ve 5 
i~1(2) i+ a6 Fae oF 


The properties of this function have been thoroughly studied (see Morse and 
Feshbach’ {pp. 782-783]). Using this relation and putting p = (1 — r)/2, (-I1< 


r< 1), we obtain 
: ; . } rl 
fa, t io: ee F —— 2) | (rye 


where 7,'(r) = 1, T7;*(r) 3r, Tor 3/,(5r? — 1), 7'3'(r) 
Here, if we use the recurrence relation: (27 + 1). — r2)T7;),(r) = i + LD) 
ia + 1)P,.,(r), the above formula becomes 
. (—1)' 
f{ - > 
JC, t) p+ 2 = 


i l - 


’ 
}Pealr) — Paalr) fe 


where P(r) represents a Legendre polynomial: Pp») = 1, P: = r, P2 = '/2(8r? — 1), 
P; = | ,(5r? — 3r), ete. For ¢ = 0, the partial sum of the first n terms of equation 
(7) is (—1)"""(P,4 — P,)/2 (n > 3). By using a proper integral expression (see 
later part), it can easily be shown that, if —1 <r < 1, the partial sum tends to zero 
as n goes to infinity. For ¢ > 0, the series (7) is uniformly convergent and ténds to 
past— o. ‘The probability of the gene A’ being fixed in the population by the 
tth generation /(0, t) can be obtained by replacing p with g and r with —r. If we 
note that P,(—r) = (—1)"P,(r), the frequency of the fixed classes is seen to be 
f0,0 +f/0,0 =1=— ¥ [Pelr) — Pasele)} eg SUT OGITOANK —°*“8) 
j=0 
which is 0 when ¢t = 0 and tends to l ast{—~> ©, 

Let us now consider the probability distribution of unfixed classes. Let $(z, t) 
be the probability density that the gene frequency in the (th generation is between 
xandx+dzr (0 <2 < 1). It has been shown that, under the assumption of the 
continuous model, ¢(z, ¢) satisfies the following partial differential equation: 
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oO SF 
= ” Ww ax? {x(1 — xg}, (9) 


which is a Fokker-Planck equation for the case of the random drift.4* This 
equation has singularities at the boundaries, and no arbitrary conditions can be 


imposed there. But the moment formula which can be obtained from equations 


(4) and (5) by calculating, 


rit) 


wn — 1f(l, b) S o'x"o(z, i) dz, (10) 


suggests that equation (9) must have the solution of the form 


- , 4 ly 
+ 2 CA lee * elec 


— 


where C, are constants and X,(x) are functions of x only. In order to solve equa- 
tion (9), if we put @ « X,(x) exp } —i(¢ + Ie iN} (¢ = 1,2,3, 


.), we obtain the 
hypergeometric equation 


r(l x) + 2(1 — 2z) 


ax, dX, 
dx? j 


— (1 — wt(¢+ 2)X, = 0, 


or, putting = (1 — z)/2 such that z 


z< 1), we obtain the Gegen- 
bauer equation: 


d?X, dX ; 
(s* =~ 3) + 42 


1. 42 — (4 — 1)(¢ + 2)X 0. (11) 
dz* dz 


From the comparison of the results obtained from equation (10), it can be found 


that a Gegenbauer polynomial X, 7; 


,(z) is a pertinent solution. 
Thus 


+ ml 
T; 1\<)e 


(12) 


The coefficients C, can be determined from the initial condition that the popula- 
tion starts from the gene frequency p. 


Mathematically, 


A ~ 1 mt 
(a—-p)= L CTA), 


: l 


(13) 

where 6(x) represents the delta function. Multiplying (1 — 2?)T7}_,(z) on both 
sides of equation (13) and using the orthogonal property, 

asl :, 

+ =e ae 2(1 + 1)t 

JS (L — 22) Tn '(z) 7; 1(z) dz = bn. 1 . . 

21 + 1) 


where m in Kronecker’s notation represents zero or a positive integer, we obtain 


(14) 


242 + lt 


2{1 — (1 — 2p)*}7i_1(1 — 2p) = C, 
~} \ =p ‘ i—!] «p P 
. 21 + 1) 


(2i + 1) ” 


C, = 4na— ; 
Pl +) 


(1 — 2p). 
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Thus the formal solution is 
> (21+ 1) —-— r? 


) Wt. z 
: G+ D T)_,(r)T'_(2)e7 HOT DANE (15) 
= ] 


d(x, t) 
or, in terms of the hypergeometric function, 


o(z, t) > pit + 1) (224+ IF - t,2 + 2,2, p) X 
l 


F(1 — 6,68 4+ 2,2, ze OTDM 015)! 


’ 


For ¢ > 0, the series is uniformly convergent for x and p, since the exponential term 
approachs zero rapidly. It is interesting that this solution agrees with the “absorb- 
ing barrier solution” of Goldberg® by putting a = 8 = 0 in his formula. 

The probability that both A and A’ coexist in the population in the ¢th generation 
(Q,) is easily obtained from equation (15), by noting that dP,(z)/dz = T}_,(z) and 
P,(1) = 1: 


*! l 
} d(x, t) dx = i (zx, t) 
0 e 1 


: (4m — 1) (1 — r*) 19 2 / 
— Te. alt) l)2m Nye 


ae (2m — 1)2m 


« 


(16) 

m 
For t > 0, the series is easily seen to be convergent, and, ast—> ©, 2, goes to zero. 
For t = 0, we must show that this series converges to 1. Let Q9, , be a partial sum 
of the first n terms; then, by the recurrence relation (4m — 1)(1 — r?®)Tb,—2(r) 
(2m — 1)2m = Pom_o(r) — Pom(r), we have Q, , = 1 — P(r). By using an in- 
tegral expression of P,, ie., P,(z) = (1/) Ys 2+ V22—- 1 cos t}" dt, we can 
show that, for |r| < 1 Ps,(r) ~ 0 asn—~ o. For 


] ® 
r+Vr? — leost!*" di 
0 


[ {r? + (1 — r*) cos? t}" dt +0 (n—> ©), 
0 


T e 
Furthermore, from (16), 


sce lidee ect (17) 
Consequently, it can be seen that, from equations (8) and (17), f(l, 4) + Q, + 
f(0, 0) = 1, as it should be. 

The processes of the change in the distribution of the unfixed classes when the 
population starts from p = 0.5 and p = 0.1 are illustrated in Figures 1 and 2, respec- 
tively. In Figure 1 it will be seen that after 2N generations the distribution curve 
becomes almost flat, and the genes are still unfixed in about 50 per cent of the 
eases. In Figure 2, the initial gene frequency is assumed to be 10 per cent, and it 
takes 4 or 5N generations before the distribution curve becomes practically flat. 
By that time, however, the genes are fixed in more than 90 per cent of the cases, and 
the simplest asymptotic formula Ce~‘'/?4” may not be useful as in the case of p = 
0.5. 
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Fics. 1-2.—The processes of the change in the probability distribution of heterallelic classes, 
due to random sampling of gametes in reproduction. It is assumed that the population starts 
from the gene frequency 0.5 in Fig. 1 (left) and 0.1 in Fig. 2 (right). ¢ = time in genera- 
tion; N = effective size of the population; abscissa is gene frequency; ordinate is probability 
density. 


The probability of heterozygosis is calculated by equation (15): 


~4 = (22 + 1), : 
H, = Jo 2x(1 — x)o(z, 0) dr > 1 T;-1 (1 — 2p) X 
‘ 1 v(t + 1) 


l 
i (1 — 2%) T)_, (2) eB + DAP ge. 
a 1 


By virtue of equation (14) (put m = 0), the last integral is 0 except for i = 1. 
Hence 

H, = pe a l JF aan = 2pqe~ "1/2" = Hye 1/2 (18) 

t Pq 2 3 <Pq of , 

showing that the heterozygosis decreases exactly at the rate of 1/(2N) per generation. 
This is readily confirmed by a simple calculation: Let p be the frequency of A in 
the population, where the frequency of the heterozygotes is 2p(1 — p). The 
amount of heterozygosis to be expected after one generation of random sampling of 
the gametes is 


E{2(p + dp) (1 — p — bp)} = 2p(l — p) — 2Ep)? = 
' i 
2nA1 — p) — 2 ON = (1 — oN 2p(1 — p), 


as was to be shown. 
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AN X-RAY ANALYSIS OF CHROMOSOME DUPLICATION* 
By Kart Sax AND Epwarpb D. KING 


ARNOLD ARBORETUM AND BUSSEY INSTITUTION, HARVARD UNIVERSITY, JAMAICA PLAIN, 
MASSACHUSETTS 


Communicated January 5, 1955 


The induction of half-chromatid aberrations by X-rays provides evidence re- 


garding the multiple nature of the chromosome and may be of considerable signifi- 


cance in the analysis of chromosome duplication. The occurrence of half-chro- 
matid exchanges was first described by Swanson! in 1943. These aberrations, 
found in pollen-tube mitoses of Tradescantia, were of sporadic occurrence. Re- 
cently Crouse* has found that practically all the X-ray-induced aberrations induced 
at the first meiotic metaphase stage in Lilium involve half-chromatid breaks and 
exchanges. Breaks in one of the half-chromatids of each of two chromatids are 
followed by reciprocal translocation to produce a chromatid bridge at anaphase. 
The close association of the coiled half-chromatids prevents the terminal separation 
of the two anaphase chromosomes (cf. Fig. 1). 

Occasional half-chromatid aberrations have been found at anaphase in the divi- 
sion of the microspore nucleus of Tradescantia following X-irradiation. The 
dosages commonly used (50-150 r) produced considerable stickiness of the chromo- 
somes, so that accurate analysis of chromosomal aberrations could not be made 
until about 6 hours after raying. By reducing the dosage to 25 r and keeping 
the inflorescences at 3° C. during irradiation, it was possible to obtain clear figures 
of mitoses and induced aberrations as early as 3 hours after raying. Irradiation at 
3° C. was done to compensate for the low dosage, since it has been shown that the 
chromosome aberration frequency can be increased by raying at low temperatures. * 

Chromosomal aberrations induced at 4—6 hours before anaphase consisted almost 
entirely of half-chromatid bridges at anaphase. Presumably, these chromosomes 
were irradiated at prometaphase or metaphase. At these stages the two .ister 
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chromatids have developed into sister chromosomes, which later pass to opposite 
poles at anaphase. The visual evidence provides little or no indication that each 
of the sister chromosomes is bipartite, but the X-ray evidence proves that each 
consists of two half-chromatids. Breaks in a half-chromatid of each sister chromo- 
some are followed by a reciprocal interchange or a lateral fusion between a half- 
chromatid of one chromosome and a half-chromatid of the other. It has not been 
possible to determine which type of aberration is produced, but, judging from the 
types of chromatid aberrations induced at prophase, both occur. In either case 
a bridge is formed at anaphase with no free fragment. If there were a reciprocal 
exchange of half-chromatids, the daughter chromosomes at anaphase would not be 
able to separate because of the close association of the half-chromatids in a common 
coil (see Fig. 1). If a dicentric half-chromatid were formed, it would produce a 
bridge at anaphase and the acentric half-chromatid fragment would be held in 
close association with the unbroken sister half-chromatids. 

As shown in Table 1, the transition between half-chromatid aberrations and 
chromatid aberrations appearing at anaphase occurs between 6 and 8 hours after 
irradiation. After about 8 hours only chromatid aberrations appear at metaphase 
and anaphase. The chromatid aberrations are induced at prophase when the 


TABLE | 


TRANSITION FROM Hatr-CuromMatTip TO CHROMATID ABERRATIONS IN T'radescantia M1icROSPORES 
FoLLOWING X-RAY Dose or 25 r at 3° C. 


Tyres AND FREQUENCIES OF ABERRATIONS 


HocuRS APTER Torat Half-Chromatid ~—-Chromatid 

EXPOSURE CHROMOSOMES Bridge Exch Per Cent Iso Exch Ye Per Cent 
j SSS 67 6 8.2 0 0 0 0 
5 366 20 0 5.6 0 0 0 0 
6 108 22 0 54 0 0 I 0.2 
7 786 33 2 $5 4 0 0 1.0 
Ss 1386 6 0 0.4 52 6 23 10.2 
10 324 0 0 0 14 5 } Cx 
12 {SO 0 0 0 16 l 5 + 6 


chromosomes are clearly differentiated into two sister chromatids. One or both 
may be broken at any given locus to produce deletions of one or both chromatids, 
or broken ends of chromatids in different chromosomes may unite to produce re- 
ciprocal translocations or dicentric chromatids. Only chromatid aberrations are 
found in microspores fixed at 10-27 hours after irradiation. At 27-30 hours some 
chromosome aberrations are found, and the transition from chromatid to chromo- 
some aberrations continues for about 10 hours, varying with radiation dosage and 
environmental factors such as temperature and light. After the transition only 
chromosome aberrations are found.‘ 

“The “chromosome aberrations” are induced during the resting stage. The breaks 
and unions induced at this time appear at metaphase and anaphase as chromosome 
deletions in the form of double rods or rings, centric ring chromosomes, and dicen- 
tric chromosomes. ‘These aberrations involve the chromosome as a single unit, and 
the X-ray evidence would indicate no bipartite structure at the resting stage. 
The chromosome aberrations continue for 5-8 days after raying the microspores, 
depending upon environmental conditions. Aberrations induced at the preceding 
meiosis are rarely recovered in microspore nucleus divisions because of the lethal 


effect of the deletions. 
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Fic. 1.—Types of chromosomal aberrations induced by X-rays at various stages 
of the nuclear cycle of Tradescantia microspores. The breaks induced at the resting 
stage involve both sister chromatics at any one locus and may lead to the formation of 
rod or ring deletions or ‘ov exchanges between chromosomes or chromosome arms to 
produce centric ring ana dicentric chromosomes. .Breaks induced at prophase may 
involve one or both of the sister chromatids at a given locus, which may lead to the 
deletion of a single chromatid or to the deletion of both chromatids. The isochromatid 
deletion is usually accompanied by lateral fusion to form a chromatid bridge and fused 
deletions at anaphase. Exchanges between chromatids of different chromosomes of 
chromosome arms results in chromatid rings and chromatid dicentrics, with the accom- 
panying deletions. Aberrations induced at prometaphase involve only one of the half- 
chromatids of each sister chromosome and result in bridges with no free fragment at 
anaphase. 

Camera lucida drawings of some of these types of aberrations are shown in sections 
1~7 of the figure. Chromosome aberrations: /, deletion; 2, dicentric. Chromatid 
aberrations: 3, isochromatid; 4, chromatid dicentric. Half-chromatid aberrations: 


~ 


5, 6, ?—bridges found at anaphase by exchange or lateral fusion of half-chromatids. 
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Thus the X-ray-induced aberrations would indicate that the chromosomes of the 
Tradescantia microspore are single threads during the resting stage, double at 
prophase, and four-partite at prometaphase and metaphase. But, if the chromo- 
some is four-partite at metaphase, each of the daughter chromosomes must be two- 
partite, and the chromosome must go into the resting stage as a two-strand thread. 

There are several lines of evidence to show that the X-ray effects are not suffi- 
ciently restricted in space to serve as a precise tool for analyzing the finer details 
of chromosome structure. A single “hit” can break one or both of two sister 
chromatids at anaphase when they are microscopically differentiated into two sepa- 
rate threads and are separated by at least 0.2 u. During early prophase, when 
the sister chromatids are close together, both are likely to be broken at any one 
locus, but as prophase progresses and the sister chromatids become more separated, 
the frequency of isolocus breaks decreases and the frequency of aberrations in- 
volving a break in only one of the two sister chromatids increases. Lea® has esti- 
mated that 15-20 ionizations are needed to produce a chromosome break and that 
a single effective “hit” covers a diameter of 0.9-1.3 u. 

The bipartite structure of the resting-stage chromosome is indicated by indirect 
cytological evidence. Most cytologists are agreed that the anaphase chromosome 
consists of at least two chromatids which are intertwined in a common coiled 
structure. The X-ray evidence shows the metaphase chromosome to be four- 
partite, and thus the anaphase chromosome must be bipartite. Yet the resting- 
stage chromosome responds to X-ray breakage as a single unit. It is not surprising 
that the resting-stage nucleus should react to X-rays as a single unit in view of the 
relatively large area effected by a single “hit.” But why do each of the sister 
chromatids at metaphase react as a bipartite structure when the half-chromatids of 
each chromatid must be very closely associated? This association of half-chroma- 
tids at metaphase certainly is closer than that of sister chromatids at prophase, 
when X-rays can break one or both of the threads at the same locus. The half- 
chromatids at metaphase would be expected to be more closely associated than 
are the sister threads in the resting-stage chromosome; yet in the latter stage both 
threads are invariably involved in X-ray-induced aberrations. 

The metaphase chromosome is four-partite, but only one of each of the sister 
half-chromatids is involved in the X-ray-induced aberrations. If both were 
broken at the same locus, we should get simple deletions; but such aberrations are 
very seldom found. It is possible that such aberrations do occur but that the 
deletion is held in place by the chromosome pellicle and does not become a free 
fragment until the chromosome sheath disintegrates at telophase or interphase. 
At metaphase each of the two sister chromosomes presumably is inclosed by its 
own pellicle, vet half-chromated exchanges do occur between sister chromosomes. 
Exchanges between both half-chromatids of the two sister chromosomes might 
be expected to occur if both are broken, but these types of aberrations are not 


found. Chromatid types of aberrations do not appear at metaphase and anaphase 
until 7 or 8 hours after irradiation. 


The peculiar behavior of the metaphase chromosomes in response to X-irradia- 
tion may be related to the time of chromosome duplication. The evidence from 
the types of X-ray-induced aberrations would indicate that the duplication of the 
chromosome occurs shortly before metaphase. According to Swift,? the DNA 
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content of the T'radescantia microspore nucleus is doubled just prior to nuclear 
division. However, autoradiographic analyses of Tradescantia microspores made 
by Taylor*® show that P*? is incorporated in the chromosomes at about. 26-32 hours 
before metaphase. This is the time when the X-ray-induced aberrations shift 
from chromosome types to chromatid types. Howard and Pelc® have found that 
P**-labeled DNA is incorporated in the chromosomes of Vicia faba root tips during 
this transition period, which in their material was 10-20 hours before metaphase. 
It is possible that chromosome synthesis begins at the onset of this transition period 
but that chromatid replication is not complete until some later stage. 

The behavior of the meiotic chromosomes in response to X-irradiation is similar 
to that in postmeiotic mitosis. Crouse found chromatid bridges and free fragments 
at the first meiotic anaphase following irradiation of pachytene stages. Each of 
the two homologous chromosomes paired at pachytene consists of two chromatids, 
a stage comparable to prophase of mitosis. Irradiation of the first meiotic meta- 
phase resulted in half-chromatid bridges at anaphase, indicating the four-partite 
nature of each of the two homologous chromosomes. Meiotic metaphase is the stage 
of chromosome duplication as detected by X-ray breakage and is comparable to 
prometaphase in mitosis. We have also obtained half-chromatid bridges by ir- 
radiating meiotic chromosomes at, or just prior to, the first meiotic metaphase in 
Tradescantia. Aberrations induced at prometaphase and metaphase of root-tip 
chromosomes of Tradescantia also consist almost exclusively of half-chromatid 
bridges, indicating the similarity of metaphase chromosomes in the somatic, meiotic, 
and haploid cells. 

DNA synthesis during meiosis shows Jittle relation to chromosome duplication 
as indicated by the nature of X-ray-induced aberrations. The indirect cytological] 
and radiation evidence indicates that the leptotene chromosome is comparable with 
the resting-stage somatic chromosome, that the pachytene stage in meiosis is 
comparable to prophase in mitosis, and that the first meiotic metaphase is com- 
parable to prometaphase in mitosis. Autoradiographs of meiotic chromosomes of 
Tradescantia show an increase of DNA at pre-leptotene but no further increase 
during the entire meiotic cycle, according to Taylor.’ Plaut,'® working with 
Lilium, found no evidence from P*? autoradiographs that any new DNA was in- 


corporated in the meiotic chromosomes between leptotene and the end of the two 
meiotic divisions. Sparrow,'! however, found an increase in DNA synthesis during 
pachytene and diplotene in Trillium—stages of meiosis comparable to prophase 


in mitosis. 

If chromosome replication involves identical developmental stages of both sister 
chromatids, it is difficult to explain why only one is broken by irradiation at meta- 
phase. If, however, the original chromatid produces the daughter chromatid in a 
sequence of chemical changes, the newly found chromatid might be particularly 
resistant or susceptible to radiation damage at certain stages of development. 
Perhaps the evidence from DNA synthesis might be of some aid in the analysis of 
chromosome duplication, but the differences in magnitude between molecular syn- 
thesis and gross chromosome duplication would appear to be too great to justify a 
direct comparison. 

If, as many competent cytologists maintain,'? the somatic anaphase chromosome 
is four-partite, any relationship between chromosome duplication and breakage by 
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X-rays must. be secondary. It is possible that chromosome duplication occurs 
one cell generation before the two components can be differentiated by the induction 
of half-chromatid breaks and two cell generations before they separate at anaphase. 
Perhaps the actual duplication could be pushed back even further in chromosome 
synthesis, until DNA duplication is synchronized with chromosome duplication. 
If so, we would have to consider the chromosome as a multistranded unit at all 
stages in the nuclear cycle, emerging as a discrete unit only after many replications 
at the molecular level. 

Summary.—In response to X-ray-induced aberrations, the Tradescantia micro- 
spore chromosomes react as single threads when irradiated during the resting stage, 
as double threads when rayed at prophase, and as four-partite strands when rayed 
at prometaphase and metaphase. The breaking of only one of the half-chromatids 
at metaphase, but of both at the resting stage, may be related to the nature of the 
half-chromatids at the time of duplication. The time of chromosome duplication 
as indicated by the types of X-ray-induced aberrations is compared with autoradio- 
graphs of DNA synthesis in microspore development and in the meiotic cycle. 

* This work was supported in part by contract with the United States Atomic Energy Com- 
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BREAKAGE OF CHROMOSOMES TIN TRADESCANTIA WITH A CALCIUM 
DEFICIENCY* 
By DALE STEFFENSEN 
BIOLOGY DEPARTMENT, BROOKHAVEN NATIONAL LABORATORY, UPTON, LONG ISLAND, NEW YORK 


“Spontaneous” chromosome breakage in T'radescantia paludosa increases when 
plants are grown in decreasing amounts of magnesium.' Since calcium and mag- 
nesium show certain biochemical similarities and are components of chromosomes 
and nucleoprotein,® it seemed of interest to determine the effects of a calcium de- 
ficiency on chromosome stability. The present paper deals with attempts to 
demonstrate the importance of calcium in the structure and stability of chromo- 
somes, since the motivating hypothesis of this work states that bivalent cations 
fulfil a structural requirement within nucleoprotein. 

Materials and Methods.—The cultural methods were the same as those described 
previously.' Clonal divisions from a single plant of Tradescantia paludosa (clone 
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5 of Sax) weighing about 5 gm. each were placed in solution | of Hoagland and 
Arnon.’ The concentration of calcium in the cultural solutions was adjusted to 
three different levels. The first solution contained 2.5 p.p.m. of calcium. Ten 
clones were grown in this solution for three months, after which the plants were 
transferred to solutions lacking calcium altogether. The second (suboptimal) 
solution, again with ten clones, contained 10.0 p.p.m. of calcium for the first three 
months and 2.5 p.p.m. thereafter. The third solution contained 250 p.p.m. of 
calcium, an optimal supply according to Hoagland and Arnon. Fourteen plants 
were started in this solution. For ten of these the pH of the initial cultural solu- 
tion was 5.8. For the other four plants grown in optimal caleium, the pH of the 
solution was adjusted weekly to 5.5 and the solution replaced every two weeks. 
The other solutions were replaced every six or seven weeks. During the interim 
between replacements, none of the pH readings of these solutions exceeded 7.6. 

Buds were collected during the sixth and seventh months of incubation. ‘These 
were fixed in a 3:1 solution of ethanol and glacial acetic acid, and subsequently the 
anthers were stained with propiono-carmine. 

Results.—Anthers from plants grown in the lowest calcium concentration con- 
tained only abortive cells, none of which were undergoing microspore division. 
Chromosoma! alterations could not be analyzed in these anthers. Anthers from 
plants grown in optimal and suboptimal calcium solutions showed complete divi- 
sion cycles. From these, the frequency of chromosomal aberrations was deter- 
mined. The frequency was much higher in the cells of the plants grown in the sub- 
optimal level of calcium. 

Variations from the normal chromosome constitutions of six haploid chromo- 
somes of nearly equal length with medial or submedial centromeres were found 
in anthers from plants grown both in optimal and in suboptimal calcium. The 


types of variation and their respective frequencies are shown in Table |. The 


constitutions of the spore types revealed much concerning their relationship in time 


TABLE 1 


FREQUENCY OF CHROMOSOME ABERRATIONS AT METAPHAS® OF THE First Micro 
SPORE Division IN 7T'radescantia 
SUBOPTIMAL 
CALCIUM OrpTiMaAL CALCIUM SUBOPTIMAL 
No No OprTimat 
Obs. Per Cent Obs Per Cent CALCIUM 
Total cells* 105 1,390 
Six normal chromosomes + 
centric or acentric fragments 17 , 2 
Chromosome fragments (5 
chrs. with centric chr. 4 
frag.) 
Dicentrics, inversions, and 
acentric rings 
Chromatid exchanges 
Chromatid deletions 
Six normal chromosomes + 
micronucleus 
Seven normal chromosomes 


Total cells with aberra- 
tions 10 9.88 14 1.01 9. 7T 
* Nine slides (buds) examined in each case 


t A comparison between the total number of aberrations of optimal and low calcium cells gave a 
chi-square value of 84.59, with DI 1 P <0.001 
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of occurrence. Types 1, 6, and 7 originated prior to microspore formation; types 
f and 5 probably arose during development of the microspore; types 2 and 3 could 
not be determined. 

The frequency of abnormal chromosomal! alterations was higher in anthers from 
plants grown in suboptimal calcium. This was especially true for spore types 1, 
2, and 5; also, types 4 and 7 occurred only in plants grown in suboptimal calcium. 
The chi-square test indicated a highly significant difference between the two 
groups of plants (? < .001). 

The number of micronuclei was determined in microspores at a stage Just prior 
to prophase. The results of these observations are shown in Table 2. Spores from 
plants grown in suboptimal caleium contained 4.42 times as many micronuclei as 
spores from plants grown in optimal calcium. The presence of micronuclei in the 
microspores reflected the occurrence of anomalous chromosome behavior at a pre- 
vious spindle figure. Either chromosomal! aberrations resulting in acentric chroma- 
tid formation or delayed chromosome separation could have accounted for the ob- 
served event. The variability in size of the micronuclei, together with the evidence 
presented in Table 1, suggested that structural alteration within chromosomes 
rather than faulty chromosomal separation was responsible for the origin of the 
majority of micronuclei. 


TABLE 2 
FREQUENCY OF Muicronuc Er (Mten.) IN Microspores PrRiIoR TO PROPHASE 


Orptimant CaLcium SUBOPTIMAL CALCIUM 
Per 100 Total per Per 100 Total per 
Cells 100 Cells ‘ Cells 100 Cells 


Total no. buds 
Microspores without 
micn, 9, 680 2,954 
One micn. per cell 82 0.840) 03 3 O4) 
0.881 , 
Two mien. per cell 2 0 0208 13 0.424 
Aborted microspores 303 2.02 760 19.90 


3.89 


Totals 10,067 3,820 


$42 = Difference between optimal and suboptimal calcium 


A third set of observations was made to determine the frequency of appearance 


of micronuclei in postdivision microspores or pollen. The number of pollen grains 
with 0, 1, 2, and 3 micronuclei from both groups of plants is recorded in Table 3. 
Here, again, differences between the two treatments were noted. Pollen grains 
containing micronuclei were 17.8 times more frequent in anthers from plants grown 
in suboptimal calcium than in anthers from plants grown in optimal caleium 
(Table 3). Abortive pollen grains were present in all anthers. 

From Table 3 it can be seen that the number of pollen grains with micronuclei 
and the number of abortive grains varied from bud to bud. There is an obvious 
correlation between the frequency of these two types of cell within individual buds 
from plants grown in suboptimal calcium. Statistical tests show that no such cor- 
relation exists in buds from plants grown in optimal calcium. The relationship 
between the numbers of micronuclei and pollen abortion in buds from suboptimal 
calcium is presented graphically in Figure 1. Using the standard regression 
method, the equation for the line in Figure 1 is y = 26.7 + 2.24(2 — 12.4), and the 
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TABLE 3 
FREQUENCY OF MiIcRONUCLEL IN EARLY BtNUCLEATE POLLEN 
CELLS 


Wirnot MickONUCLEL PER CeLt CELLS 
FLOWER Bubp l , ; ABORTEI 
Optimal 


Calcium: 


Totals [2.4 
Per 100 ¢ 


Suboptimal 
Calcium 


5O6 
000 , : 6 
1,021: x 0 
613 j 3 0 
1,154 : ) 0 
1,220 2 0 
Totals 6,380 7 , 
Per 100 cells Q 62 0.812 0 098 
Total per 100 cells 11 54 28 |‘ 


7.28 = Difference between frequencies of micronuclei in optimal and suboptimal 
calcium cells 


REGRESSION OF MICRONUCLE! AND POLLEN ABORTION 


IN CALCIUM DEFICIENT TRADESCANTIA 
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slope of this line is 0.46, which differs significantly from one and zero. Chromatin 
duplication in one pollen grain of a tetrad would presumably produce a corre- 
sponding chromatin deficiency in another grain, causing it to abort. This could 
explain the relationship between the numbers of micronuclei in pollen and of pollen 
abortion. 

Discussion.—The salient point of these experiments is that the frequency of 
chromosome fragmentation and the presence of micronuclei are functions of cal- 
cium concentration. This contention is further supported by unpublished experi- 
ments which involve the study of several different calcium levels. 

Further evidence indicates that plants grown in decreasing concentrations of 
calcium are increasingly more sensitive to X-rays. Of particular interest to the 
case in point is a differential effect between types of chromosomal aberration ob- 
served following X-irradiation of buds grown in low calcium. The proportional 
increase of chromosome fragments and interstitial deletions is higher than the in- 
crease of dicentries and centric rings. These observations were twice repeated, 
and each time differences were significant. These unpublished data support the 
conclusion that calcium is directly involved in the chromosome-breakage process in 
Tradescantia. 

Bernstein and Mazia‘ report that chelating or sequestering agents, such as citrate 
and sodium phosphate, following treatment with distilled water produce a deoxy- 
ribonucleoprotein particle approximately 4300 A by 250-300 A in size. The ions 
removed by sequestering presumably hold the particles within the chromosome to- 
gether by chemical bonds. In macromolecular states, it is indicated that caletum 
and magnesium form bonds with deoxyribonucleic acid (DNA)® and with protein.® 
Other studies show that cations provide stability for viruses’ and for DNA.* 
The exact nature of bonding in DNA, nucleoprotein, or chromosomes is not known. 
But a priori bonding could occur both within and between the key groups of mac- 
romolecular constituents. Important is the fact that calcium and magnesium pre- 
fer to form complexes with oxygen, in particular with the oxygen of phosphate 
groups. ® 

The above biochemical evidence, together with the finding that a deficiency of 
calcium can significantly alter both spontaneous and X-ray-induced chromosomal 
aberrations, leads to the following hypothesis: Bivalent cations, e.g., caletum 
and magnesium, are involved with complexing and binding of chromosomal nucleo- 
protein and are thereby responsible for a part of the chromosome’s structural con- 
stitution. 

Summary.—Tradescantia plants from the same clone were grown in different 
calcium concentrations. Observations from microspores and pollen showed that 
plants grown in suboptimal calcium concentrations exhibited at least seventeen 
times more chromosomal aberrations and micronuclei than did plants grown in 
optimal calcium. Coupling these data with the biochemical evidence presented, 
it was concluded that bivalent cations may bind together the macromolecular 
constituents of the chromosome and thereby provide structural stability. 


* Research carried out at Brookhaven National Laboratory under the auspices of the United 
States Atomic Energy Commission. 

1D. Steffensen, these PROCEEDINGS, 39, 613-620, 1953 

? £. S. Horning, in G. Bourne (ed.), Cytology and Cell Physiology (London: Oxford University 
Press, 1942), pp. 160-88; L. G. Allgén, Acta Physiol. Scand., Vol. 22 (Suppl. 76), 1950. 
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3D. R. Hoagland and D. I. Arnon, Calif. Agr. Exptl. Sta. Circ., No. 347, 1950. 
*M. H. Bernstein and D. Muazia, Biochim. et biophys. acta, 10, 600-606, 1953; thid., 11, 59-68, 


1953. 
’G. Jungner, Acta Chem. Scand., 5, 168-174, 1951. 
6 W. G. Crewther, Australian J. Biol. Set., 6, 597-616, 1954. 
7 Karl G. Lark and Mark H. Adams, Cold Spring Harbor Symposia Quant. Biol., 18, 171-183, 
1953. 
8S. Zamenhof, G. Griboff, and N. Marullo, Biochim. et Biophys. Acta, 13, 459-470, 1954. 
*R. J. P. Williams, Biol. Rev., 28, 381-415, 1953. 


ON THE NICKEL, COBALT, AND COPPER CONTENTS OF DEEP-SEA 
SEDIMENTS* 


By G. E. Hurcutison, R. J. Benorr,t W. B. Correr, ano P. J. WANGERSKY 
OSBORN ZOOLOGICAL LABORATORY, YALE UNIVERSITY, NEW HAVEN, CONNECTICUT 
Communicated January 9, 1955 


In two recent papers Pettersson and Rotschi,' using a magnificent series of 

analyses of deep-sea cores collected by the Albatross expedition and analyzed by 

totschi,? have suggested that the nickel content of Pacific red clay is mainly of 
extraterrestrial origin. They observed considerable differences between typical 
Pacific red clays and the series of Atlantic samples from a single core that they 
analyzed, the Pacific material ordinarily containing at least five times as much 
nickel as the Atlantic. They attribute this to the greater lithogenous sedimentation 
in the Atlantic, which masks the effect of the extraterrestrial contribution. 

Goldberg,’ in a stimulating study of the role of electrochemical precipitation in 
the sea, concludes that the nickel content of the sediments analyzed by Rotschi 
is in a broad way dependent on their manganese content. He supposes that, 
during the electrochemical precipitation of manganese oxides, positively charged 
nickel is adsorbed on the precipitate. He further points out that in manganese 
nodules cobalt and copper as well as nickel are enriched; the ratios of these elements 
to nickel in such nodules are much greater than in meteoritic material in which the 
Ni: Co ratio is at least 10:1 and the copper content much less than that of nickel.‘ 
If an extraterrestrial origin for the nickel is postulated, special hypotheses for the 
enrichment of the other two elements must also be entertained. 

Since we have been engaged in the chemical study of certain deep-sea cores, dur- 
ing which a number of analyses of nickel, cobalt, and copper have been made, 
it has seemed proper to complete Goldberg’s argument by publication of our data 
for these elements. This is particularly desirable since Goldberg’s original analyses 
relate solely to manganese nodules, while Rotschi’s great body of data does not 
include determinations of cobalt and copper. 

Three cores have been available, two from the Atlantic and one from the Pacific 
ocean. A 160-8 is for the most part a typical red clay, from 5,030 meters at lati- 
tude 16°54’ N., longitude 59°31’ W. It is 870 cm. long; twenty samples were ana- 
lyzed, but, in presenting the data, two are omitted. One of these is from a sharply 
defined band of voleanic ash and is obviously irrelevant; the other is from a man- 
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ganiferous zone which probably represents a discontinuity in the process of ordinary 
or lithogenous sedimentation. A 152-134 is a chemically typical red clay, curiously 
pale in color, from 5,080 meters at latitude 35°54’ N., longitude 62°17’ W. The 
core is 350 em. long; seven samples were analyzed. MP 35.2, our only Pacific 
core, is from 4,850 meters at latitude 10°25’ N., longitude 174°36’ W. It isa typical 
red clay. The @ore is 376 cm. long; nine samples were analyzed. 

Nickel and cobalt were determined by the method of Perry and Serfass,® using 
nitrososalicylie acid; copper was determined by the standard spectrophotometric 
method with sodium diethyldithiocarbamate. The ranges and means of our de- 
terminations, together with those of Rotschi for nickel, are given in Table 1. 


TABLE 1 
Copper, NICKEL, AND CoBALT (MG/KG) 
Ocean Core Copper Nickel Cobalt Ni:Co 
Atlan- <A 160-8 0-610; mean 109 25-122; mean 63 < 10-60; mean 25 2.6:1 
tic A 152-134 17-63; mean 41 57-142; mean 92 14-26; mean 19 ‘3:3 
Rotschi 273 ; 40-110; mean 62 ; 
Pacific MP 35.2 218-447; mean 320 114-535; mean 364 93-310; mean 194 1.9:1 
Rotschi, all cores 30-786; mean 267 e 
While the latitudinal variation in manganese noted for the Atlantic by Hutch- 
inson® and for the Pacific by Goldberg’ may make detailed comparisons difficult, it is 
reasonable to conclude in the light of all the data that Pacific red clays are con- 
siderably richer in all three elements than are Atlantic. The only probable ex- 
planation is that the ratio of the rate of hydrogenous or chemical sedimentation to 
ordinary or lithogenous sedimentation is much greater in the Pacifie than in the 
Atlantic. 
The analyses of MP 35.2 in Table 2 indicate that the upper part of the core is 
poorer in the three elements than is the lower. There is also definitely more 


TABLE 2 
Copper, NICKEL, AND CoBALT IN CorE MP 35.2 (mG/ka) 

Cu Ni Co 

0-2.5 cm. 225 247 107 

46-53 232 214 121 

100-110 218 259 O4 

155-161 310 235 129 

199-215 447 410 253 

259-268 326 503 218 

305-313 373 535 239 

363-374 380 516 310 

374-376 372 526 274 
manganese’ and calcium in the lower half. These differences presumably reflect a 
change in the ratio of lithogenous to other types of sedimentation; some of Rot- 
schi’s cores show similar and some much more irregular variations. Although the 
attractive hypothesis of the cosmic origin of the nickel in the Pacific red clay must 
evidently be abandoned, the variations in composition must mean something, and 
their elucidation is likely to be of immense importance to paleoceanography. It is 
reasonably certain that all the values in Table 2 represent an enrichment over any 

probable mixture of source materials in the accessible lithosphere. 

It is finally desirable to mention that the copper contents of A 160-8 are distrib- 
uted most irregularly. The single high value of 610 mg/kg might conceivably be due 
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to contamination of a short section of the core after collection. There are, how- 
ever, three levels at which no copper could be found; this was confirmed by replicate 
analysis and apparently must be valid. Such findings seem to suggest that the 
hydrogenous sedimentation of copper depends on rather delicately balanced condi- 


tions. 


We are greatly endebted to Professor Maurice Ewing, of the Lamont Laboratory, 
Columbia University, for the material from cores A 160-8 and A 152-134 and to Dr. 
J. D. Frautschy of the Scripps Institution for Oceanography for that from core 
MP 35.2. 

* This study was supported by ONR Contract No. SAR/NONR-260(00). It forms part of the 
program of the NRC Committee on the Chemical Composition of Sediments. 

t Present address: Department of Limnology, Academy of Natural Sciences, Philadelphia, 
Pennsylvania. 

1H. Pettersson and H. Rotschi, Nature, 166, 308, 1950; Geochim. et cosmochim, acta, 2, 81-90, 
1952. 

2H. Rotschi, Cahiers du Centre recherche etude oceanog. (Paris), No. 4, pp. 2-22, 1951. 

3 E. Galdberg, J. Geol., 62, 249-265, 1954. 

*K. Rankama and T. G. Sahama, Geochemistry (Chicago: University of Chicago Press, 1950). 

5M. H. Perry and E. J. Serfass, Anal. Chem., 22, 565-567, 1950. 

* G. FE. Hutchinson, Ecol. Monographs, 17, 299-315, 1947. 

’ The manganese and other determinations of the major constituents, made by M. Tavola and 
H. B. Wiik, of the Geological Survey of Finland, will be published in a later paper. 


BOUNDS FOR THE RATIOS OF THE CO-ORDINATES OF THE 
CHARACTERISTIC VECTORS OF A MATRIX 


By ALFRED BRAUER 


UNIVERSITY OF NORTH CAROLINA, CHAPEL HILL 


Communicated by H. S. Vandiver, December 15, 1954 


During the last fifty years the prob'em of finding bounds for the characteristic 
roots of a matrix has been studied by many authors. I myself obtained such 
bounds in a number of papers.! Now I can improve many of these theorems fur- 
ther by obtaining bounds for the ratios of the co-ordinates of the characteristic 
vectors. Since my former results are used in different fields of application, I want 
to give here a short account of the new method even before the investigations are 
finished. I intend to publish my results in detail later in the Mathematische Zeit- 
schrift in a volume dedicated to the memory of I. Schur. 

Let A = (a,,) be an arbitrary square matrix of order n. Let w be one of its char- 
acteristic roots and ry = (x), %,...,2,) a characteristic vector belonging tow. As- 
sume that 


= max /2, = Secvcc ee t,o. - 13) 


We set 


~ { 
> lay 


A=] 
Ate 
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It is well known that w lies in the interior or on the boundary of the circle 
~ 4,| S f,. (2) 


I proved that w lies more exactly in the intersection of the circle (2) with the oval of 
Cassini 


2 — Ge! & PP, (3) 


Moreover, I proved that in certain cases a smaller oval than (3) must contain the 
root w. 

It will be shown in this paper that the circle (2) can be replaced by a smaller con- 
centric circle if at least one of the elements of the rth row, say, a,,, is different from 
0 and if the circle 


(4) 


has no point in common with (2). Similarly, the ovals can be replaced by smaller 
ones. The new results make it sometimes possible to compute easily a character- 
istic root as accurately as we please. 

THEOREM 1. Let R be the circle (2) or a closed subregion of it which contains w. 
Forwx = 1,2,..., n, let d, be the minimum distance of a,, from R. We set 


$O=1 f d.=0 “ef SE ot 


P 
t,,() = min (1, ‘, 
d, 


and, for d, > 0, 


l - ' 
t,,° = min ) i. ( ag) + Dd laalta’™’ ) (» = 2, 3, 
{ d, A=1 f 


Ache 


Then, for each v, 


ee ae POE Re a 2 sag Me ee 


« 7 


Proof: Inequality (8) is trivial for those values of « and » for which 4, = 1 
because |z,) S |2z,| by relation (1). Since r is a characteristic vector belonging 
to the root w, it follows that x, re, ..., 2, are a solution of the system of linear 
equations 


n 


(w — @)r, = p dF, fy es Py See 
p=l 
ueer 


From the «th of these equations for x # r we obtain 
w— a,,| |2z,| S P,\2,|; 


hence, since w is a point of PR, 


diz | a 


« + 7 


By formulas (6) and (10) our theorem is proved for » = 1. 
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We now assume that the theorem is already proved for y — 1; hence, by relation 

(8), 
i a (A= 1,2,..., 9; #1). (11) 


Substituting relation (11) in the «th of equations (9), we obtain 


' 1) 
an g*” Ir|, 


n 
< la,||z,| + > 
A= 
d 


Hx, 


' 
1 
r 


n 
| (v—1 | | (») ‘ 
ry ( Oy, oa me ay tyr : )= I, . {xk = a 2. 
A=1 
At«r 


and the theorem is proved by induction. 
If we substitute relation (8) in the rth of equations (9), then we have 


n 


Ory |Z) 
(12) 


Let us denote thd circle (12) belonging to R by G(R). Then we have proved 

THEOREM 2. Each characteristic root of A whose rth co-ordinate of a characteristic 
vector is the absolute greatest and which lies in R must lie in the intersection of R and 
G(R). If this intersection is empty, then R does not contain such a root. 

Theorem 2 can sometimes be used to compute a characteristic root as accurately 
as needed by excluding more and more regions which cannot contain characteristic 
roots. 

Let R be the circle (2). Assume that at least one of the elements of the rth row 
of A, say, a,,, is different from 0 and that the cirele (4) has no point in common with 
(2). Then 


a,, — a,,| > P, + P,, d, = |a — a,,| — P, > P,. 


AK 


Hence t,,'”. < 1 by relations (6) and (7), and G(R) is smaller than R. 

If we denote the circle (4) by C, (x = 1,2,...,), then each characteristic root 
of A lies in the union of G(C)), G(C2), .. . , G(C,), hence in the union of G{G(C))}, 
GiG(Cy)},..., @IG(C,)}, and so on. 


! Alfred Brauer, ‘Limits for the Characteristic Roots of a Matrix,’’ Duke Math. J., 13, 387-395 
1946; II, ibid., 14, 21-26, 1947; III, ibid., 15, 871-877, 1948; IV, zbid., 19, 77-91, 1952; V, ibid., 
19, 553-562, 1952; VI (with H. T. LaBorde), to be published zbid., Vol. 22, 1955; “Matrices with 
All Their Characteristic Roots in the Interior of the Unit Circle,’ J. Elisha Mitchell Sci. Soc., 68, 
180-183, 1952; ‘Bounds for Characteristic Roots of Matrices,’’ Naél. Bur. Standards Appl. Math. 
Ser., 29, 101-106, 1953; “Uber die Lage der charakteristischen Wurzeln einer Matrix,” J. reine u 
angew. Math., 192, 113-116, 1953. 





AN EXTENSION OF THE MAXWELL THEORY 
By R. L. INGRAHAM 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF CONNECTICUT, STORRS, CONNECTICUT 
Communicated by Oswald Veblen, January 6, 1955 


1. Jntroduction.—This note, summarizing an article to appear elsewhere, con- 
cerns the elimination of infinities from classical electromagnetic theory. In the 
past this has been accomplished by noninvariant devices like giving the particle 
a finite extension in 3-space! or by the admixture of unphysical elements like ad- 
vanced fields.2, The present method is essentially different from these. It consists 
of the extension of the conventional theory to the whole of a “natural”? domain of 
definition suggested by the conformal form-invariance of the Maxwell laws. This 
introduces a new degree of freedom, A, alongside space-time. In this preliminary 
note the field of a particle in the extended theory is obtained and its formal proper- 
ties briefly investigated, with especial attention to the question of infinities. Specu- 
lation as to the possible physical meaning of the new co-ordinate \ is postponed. 
For \ + 0 the field of the particle is shown to be everywhere finite (A operating as 
a Lorentz invariant cut-off). The special subspace \ = 0 (precisely the space of 
the conventional theory) is then excluded as representing an unattained limiting 
situation in nature. Also touched on are the inertial versus noninertial character 
of motion in the electromagnetic field, the self-force problem, energy-momentum, 
and motion in the Coulomb field as they appear in the extended theory. 

2. The Extended Field Equations.—The Maxwell laws in empty space are known! 
to admit a larger class of equivalent observers than that of the Lorentz observers; 
that is, they are form-invariant not only under Lorentz transformations but under 
all conformal transformations #? = f?(x") of space-time.‘ The conformal group is, 
in facet, the largest group for whose equivalent frames light-velocity is constant; 
it has the physical interpretation® of uniform relative accelerations of cartesian 
3-frames (the Lorentz group forming a subgroup of zero relative accelerations). 
Now the conformal group is a group of sphere transformations: it takes space- 
time spheres 


/ 


— r)? — c(t’ — {)? = —j? (2.1) 


(r 


into other space-time spheres and nullspheres (A = 0) into nullspheres,* and its ac- 
tion on points of space-time is obtained by identifying these with the nullspheres of 
which they are the centers. Let‘ the 5-ple 2*, where x" = the real Lorentz co- 
ordinates of the center, and x° = \, the (directed) radius, of the sphere (2.1), be 
called natural co-ordinates for the sphere. (Spacelike spheres are excluded for 
physical reasons; hence the minus prefixing A? in (2.1).) Then the Maxwell equa- 
tions have the aspect of a set of laws form-invariant under conformal sphere trans- 
formations #* = f%(x*) but defined only over the subspace \ = 0 of nullspheres. 
We extend the theory by extending its domain of definition to all spheres in space- 
time. Then the fields will be functions of 2” and i, and a 5-tensor formalism is 
necessary to write equations covariant against transformations of the five x*. The 
prescription for writing down the extended laws is that they be formally identical 
with the old laws in terms of the corresponding 5-tensors; there is no arbitrariness 
in this prescription. 


165 
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The angle dé between neighboring spheres x* and x* + dx* turns out to be 
dO? = Ya, dx* dx® = ——*(dr? — c? dt® + dd*). 


d@ is the (dimensionless) invariant of spherespace which replaces the length in- 
variant ds of special relativity; the metric’ y,, of spherespace is defined by (2.2). 
Then, if Pas = —/'gq is the 5-force, the laws: Fs is a curl, and the divergence 


. . 2 « > ° ® 
yo 2 Oy! ty y™ F;.) vanishes, become, in natural co-ordinates 


(a) Chad al + Oo, F; ” +- OF nn = 0; (b) Om!’ sn ists 0, bm = OsF mn; 
(c) Ok", = — Osl'sn + A7'Fs,; (d) Ons" = 0. (2.3) 


Here 0, means 0/Oxr* and F™, = g”’Fy,, Where gm, is the Lorentz metric. The 
effect of the extension has been to adjoin a 4-vector F;, to the 4-force F,,,, and to 
make both these fields functions of (r", A). On the nullspheres (A = 0) we take 
F’;, = 0; then all the terms involving an index 5 in (2.3) fall away, and we get back 
the old theory for F,,,.(2”). 

3. The Field of a Particle.—It can be verified that Fg = Of 3 — Osl,, where 


€ Um 


tr Rv,’ 


oe Ff, = Q; R°R, = —\» (3.1) 
is a solution of (2.3). Here (y"(w), 0) is the path of the source, w any invariant 
running parameter (N.B.: the path points are all nullspheres); the velocity 
vy” == dy"/dw, R? = x? — y’(w), indices are raised and lowered with the Lorentz 
metric and —e is the charge on the source (rationalized units). The field is 


evaluated at the sphere (2”, \), and the last equation fixes the retarded time as a 
function of (2”, A). On the nullspheres we get just the classical Liénard-Wiechert 
solution for a point particle in arbitrary motion, so (3.1) is the extension of this 
solution to the whole spherespace.’ Incidentally, the retarded time condition has 
the invariant meaning that the field sphere (2”, A) and the source nullsphere 
(y™(w), 0) be orthogonal at y* = the retarded time. When dA = 0, this degenerates 
into the old condition that R? be a null 4-vector. The force is computeu io be 


( ” A a 
Fi 2 me (a. _ rn), F,, = («, _ rs) (3.2) 
lr R? | ie 4nr R? R 


where the acceleration’ a” = d°y"/dw?, R= R’v,, A = Ra, — v’v,,and [| means 
the antisymmetric part, e.g., 2R;..0,) = Rava — Rivm. The main result of this 
theory now follows immediately. The force (3.2) could only be infinite where 
R = Rv, = 0. But, since the source velocity v? must be taken timelike, the 
orthogonal vector R? would have to be spacelike or a zero vector; thus R’R, > 0. 
But this contradicts the last equation of (3.1) for \ + 0. Therefore, the field of 
the particle is infinity-free on the nonnull spheres.'° There are good reasons for 
now excluding the nullspheres from the observable domain of the field, whence the 
force is everywhere finite." 

In virtue of the motion equations adopted (ef. (4.3)), the components of the 4- 
force F,,,, continue to have their usual identifications in terms of the electric and 
magnetic fields. In addition, there is the 4-vector G,, = F;,, with no classical. 
counterpart (if one ignores a suggestive formal analogy with the 4-current; see 
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(2.3). The energy-momentum density tensor of the field, whose divergence 
vanishes in virtue of the field equations, is 


Tig = FUP ip — “erg F (3.3) 


ag 


This leads in the usual 
way to an expression for the energy-momentum 4-vector of the field as an integral 
over a hypersurface in spherespace. The (time-constant) energy of the field at con- 


where indices here are raised with the angle metric 7 


stant A comes out to be 


—> —, 


ils Di (ED te H + G + GG) dx dy dz, (3.4) 


where the limits on the variables are — © to ©. Detailed energy and momentum 
considerations are post poned to the article to appear. 

1. The Motion Equations.-The motion equations adopted are formally identical 
with the Lorentz equations of motion, where the transposition into the correspond- 
ing 5-tensors guarantees both the conformal invariance and che extension of the 
domain of definition to the whole of spherespace. In natural co-ordinates they 
run'!* 

d*x* 
~—-A=—y Jay p (4.1) 
de dé 
where f.3 = 4re~'F'43 and « > 0 is some standard charge, hereafter taken to be the 
magnitude of the electron charge, introduced to make the physical dimensions 
correct. If one goes over to a length running parameter s via the noninvariant 
substitution 
ds ? 
= » (4.2) 
dé No 
where the constant length A)» depends on the path, (4.1) take their “special rela- 
tivistic”’ form, 


: dx ad? 7 a : 
= f™. . oo o“A $ + No fom ° (4.3) 
io ae ds? ds 

Consider a “classical” solution of (2.3): Fs, = 0, Fn. = 0. Then, comparing 

(4.3) with the Lorentz motion equations for a point particle P at 2” of charge- 

mass ratio (e/m)p > O and sign of charge Sgn ep in the field F,,,,(2?), we find!’ 


l € ’ ‘ 
No = of" sgn No = sgn Cp. (4.4) 


tr \m/p c* 


The noninertial character of electromagnetic motion has thus been introduced on 
going over to the framework of length geometry via the noninvariant substitution 
(4.2); in the invariant form of the motion equations (4.1) A» does not occur, and the 
electromagnetic forces become inertial. The mathematical root of this is simply 
that, whereas ds is a length, dé is dimensionless, so that, although d@ occurs to the 
—2 power on the left (ef. n. 12) and —1 power on the right, no compensating 
length A» depending on the particle need be inserted in the right member of (4.1). 
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By putting the field (3.2) into (4.3), taking R' = 0 (7 = 1, 2, 3), and, assuming 
\ small, expanding retarded quantities in power series in A, we can study the self- 
force. The self-force is infinity-free here, and there is no need to “smear”’ the source 
and perform integrations as classically; the nonzero \ keeps the particle “at a cer- 
tain distance from itself,”’ thus avoiding the instantaneous action on itself leading 
to infinities. The “electromagnetic mass” is proportional to e?//A!. For |\| & 
'\o| we get agreement with the classical result in sign and order of magnitude 
(identifying \ and a, the 3-sphere radius), but for |A) > |A»o) (A of atomic dimen- 
sions for P an electron), due to the influence of the term Aof"s (dA/ds) in (4.3), the 
self-force vanishes. There is some doubt, however, as to whether the solution 
(3.1), and thus these results, should be taken as valid at such distances. 

For motion in the Coulomb field of a resting source, (4.3) become 


— — 


d?*r r r 
P . —Xo ( ) ~ wr ale 
8 ef (7? + A*)T 


—do(e/e) ig an 
a " r +X ; 
(r? + ?) ds 
d*d hes ( » dx‘ 
> : Ao“A $ No . moet yA 
ds? e/ (r? + A*)” ds 


These have not yet been integrated exactly. The middle equation yields the 
energy first integral'* 


\ " 24. (w/A)? Ao(e/e) BE 
V1 — v?/c? Vr? + i? mec? 


= const., (4.6) 


where the dot means d/ds. Hence the kinetic energy £,;, must be defined as mc? 
times the first term of the left member of (4.6). For free motion (f,3 = 0 in (4.3)) 
we get A? = s* + Ay”, whence VA? + (Ao/A)? = 1 and F,;, goes over into the Ein- 
stein quantity. 

The most important question here-—the physical meaning of A-—will be explored 
in the article to appear. 


' See L. Landau and E. Lifshitz, The Classical Theory of Fields (Cambridge, Mass.: Addison- 
Wesley Press, 1951), chap. viii. 

? P. Dirac, Proc. Roy. Soc. London, A, 167, 148, 1938. 

§ As a short proof, note that in the vanishing divergence equation the Lorentz metric gm» can 
be replaced by the conformal invariant density gm» /(—g fa, 

‘ Index conventions: middle Latin letters /, m, n, p, q, etc., go from | to 4; early Greek letters 
a, 8, y, 6, «, ete., go from 1 to 5. 

'See J. Haantjes, Proc. Ned. Akad. Wetenschap., 43, 1288, 1940; also R. Ingraham, Nuovo 
Cimento, 12, 825, 1954. 

* This, in fact, characterizes the conformal group. 

7 vag is the metric of a Riemannian 5-space of constant curvature + 1, or the surface of a ball 
imbedded in flat 6-space of signature (+++-—+-—). The group of conformal sphere trans- 
formations is just the group of proper rotations of this ball; cf. Ingraham, op. ctt. 

* Note that solution (3.1) is given in nonconformal invariant form and thus refers to a par- 
ticular class of Lorentz frames. The solution in the general frame is obtained by tensorially 


transforming (3.1). 
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® Since (3.1) is homogeneous of degree zero in pm, v™, and a™ in (3.2) can be taken to be dy” /dy', 
d?y™ /(dy*)*, resp. Thenv'=1,a'*=0. 

'© On the nullspheres (A = 0) one gets R = 0, where R? = 0—i.e., just the Coulomb infinity. 

 E.g., \ acts like a sort of gauge, or numerical length assigned to some standard physical rod 


by the observer at 2”. This rod can be assigned any positive value, but not zero. 
'? The left member is the covariant acceleration (d®%x%/de*) + + } (dx® /d@) (dr¥/d6). 
; 


'S ds > O assumed, whence Sgn Ay = Sgn dé. 


'" The identity drt/ds = VU X°* + (4o/A)2/ VY 1 — v2 /c? following from (2.2) has been used. 


MULTIPLICITY THEORY FOR OPERATOR ALGEBRAS 
By Ricuarp VY. Kapison* 
COLUMBIA UNIVERSITY AND THE UNIVERSITY OF COPENHAGEN 


Communicated by John von Neumann, December 3, 1954 


1. Introduction.—We shall outline, in this note, a multiplicity or unitary-invari- 
ants theory for adjoint-preserving homomorphisms of self-adjoint algebras of 
operators closed in the uniform operator topology (briefly, “representations of 
C*-algebras”). A detailed account of this theory containing a historical survey 
together with complete bibliographical references will be published elsewhere. 
We take the liberty of omitting from this note all bibliographical references. 

In its general aspect, the central problem in the study of operator algebras is that 
of associating with each operator algebra a set of entities (invariants) which com- 
pletely determine it, i.e., such that two operator algebras are structurally the same 
if and only if they have identical sets of invariants. The province in which one 
seeks invariants for any structure is indicated by experience with the struc- 
ture and the special theories developed for special cases of ‘the structure. 
Again, experience and applicability to the special cases usually indicates when a 
set of invariants for a structure is “the correct set” in the sense of relevance and 
simplicity. In the case of operator algebras, the central question just mentioned 
has two main parts. On the one hand, we may ask for invariants for the operator 
algebra as an algebraic object (algebraic invariants), and, on the other, we may ask 
for invariants for the operator algebra together with its action on some Hilbert 
space (spatial, multiplicity, or unitary invariants). In an earlier paper we pre- 
sented what we felt were “the correct algebraic invariants” for the general C*- 
algebra (and noted there that, of course, one would insist on more detailed invari- 
ants for special classes of C*-algebras such as the factors of Type II,). In this 
note we give the spatial invariants for a C*-algebra. 

Some remarks concerning the relation of these results to previous work in mul- 
tiplicity theory are in order. The finite-dimensional C*-algebras are completely 
described by the Wedderburn structure theory. The first infinite-dimensional 
results were those of Hahn-Hellinger, describing when two self-adjoint operators 
are unitarily equivalent. Later technical improvements in their description made 
it possible to apply their theory directly to commutative C*-algebras. This last 
theory contains that of the single self-adjoint operator. The next step in the com- 
mutative theory was the study by Nakano-Segal of Abelian rings of operators (a 
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ring of operators is a C*-algebra closed in the weak operator topology). The class 
of Abelian rings of operators is a subclass of the class of commutative C*-algebras 
and admits a very precise set of multiplicity invariants. At first sight, if might 
seem that the theory of Abelian rings of operators contains that of the single oper- 
ator, as does the general theory of commutative C*-algebras, but this is not the 
case. Roughly speaking, the multiplicity theory for a single self-adjoint operator 
(and a commutative C*-algebra) involves its spectrum together with an ideal ot 
Borel sets in the spectrum, the null sets of some measure, while the multiplicity 
theory for an Abelian ring of operators involves the measure algebra of some measure 
(actually, an equivalence class of measures) on its spectrum, i.e., the Boolean alge- 
bra of measurable sets modulo the null sets. In the case of noncommutative rings 
of operators, Murray and von Neumann took the decisive step by investigating 
factors (those rings of operators whose centers consist of scalar multiples of the 
identity) and developing a multiplicity theory for all but the factors of Type IIT. 
Using techniques of Dixmier, which generalize those of Murray and von Neumann 
and apply to general rings of operators (with no portion of Type III), Dve and 
Griffin generalized the Murray-von Neumann multiplicity theory for factors to 
rings of operators (with no portion of Type III). Recently, Griffin developed the 
spatial theory for rings of Type III, and his results show that, in reality, this is the 
simplest case, as far as the spatial theory is concerned. Thus the spatial invariants 
for non-Abelian rings have been determined, carrying the multiplicity theory for 
Abelian rings over to general rings of operators. The theory outlined in this note 
completes self-adjoint multiplicity theory by listing the spatial invariants for non- 
commutative C*-algebras. This is the noncommutative (but self-adjoint) exten- 
sion of the classical spectral multiplicity theory for a single self-adjoint operator. 
The invariants we set down are of the same general character as those obtained 
in the commutative (or single-operator) case and yield the commutative theory 
immediately upon specialization. These invariants are the pure-state space of the 
C*-algebra together with a descending chain of ideals of Borel sets in the pure- 
state space and a distinguished class of such ideals. There is at present no tech- 


nique for determining which systems of the above type are associated with C*- 
algebras, but this defect also appears in the classical spectral multiplicity theory 


of a single self-adjoint operator. In the single-operator situation the question is 
that of which ideals of Borel sets on a bounded subset of the real line can serve as 
the null sets of some Borel measure. 

Since spatial equivalence of operator algebras implies, in particular, algebraic 
equivalence, a set of spatial invariants must include, in some form, the algebraic in- 
variants. We actually set down spatial invariants for representations of abstract 
('*-algebras as concrete C*-algebras, something slightly more general (and useful) 
than locating invariants for concrete C*-algebras (acting on a Hilbert space). By 
so doing, we gain the technical advantage of introducing the assumption of alge- 
braie equivalence without the necessity of employing a bewildering array of ™se- 
less algebraic isomorphisms. We have developed the theory without separability 
restrictions, although the results are new in the separable case, and although our 
natural inclination is to present the somewhat less involved picture permitted 
by the assumption of separability. Our choice was guided by the desire to have 
this theory serve as the final step in self-adjoint multiplicity theory. 
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2. Basie Notions and Definitions.—We recall that an abstract C*-algebra is a 
Banach algebra over the complexes with a *-operation which is a conjugate-linear, 
involutory antiautomorphism satisfying the condition |/A*A) = |A|?, for each 
element A in the Banach algebra, with the norm on the Banach algebra. 
Each such algebra has a faithful (one-one) representation as an algebra of operators 
on a Hilbert space, closed in the uniform (norm) topology (a concrete C*-algebra). 
We shall assume, throughout, that our concrete and abstract C*-algebras contain 
a unit element. The linear functionals which are real and nonnegative on the 


positive operators in the algebra and | at the identity are the states of the (’*- 


algebra. The states form a (compact) convex subset of the space of all continuous 
linear functionals, and the extreme points of this convex set are “the pure states’ 
of the algebra. The closure of the set of pure states is a compact-Hausdorff space 
called “the pure-state space of the algebra.’”’ The obvious map from the elements 
of the C*-algebra to functions on the pure-state space is a linear, order-preserving 
isomorphism of the C*-algebra with a linear subspace of the set of continuous func- 
tions on the pure-state space. (The topology on the pure-state space is the weak- 
est in which the functions of this set are continuous.) We refer to this linear space 


, 


of functions as “the representing function system of the algebra.”’ 

Derinition 1. Jf YM 1s a C*-algebra and ¢ is a representation of A as an algebra 
Yl of operators on the Hilbert space ®, we denote by N, the ideal of Borel subsets of X, 
the pure-state space of A, consisting of those Borel sets S which are null sets of the Borel 
measure induced on X by each positive extension to C(X), the family of continuous 
functions on X, of fed, for each countably additive (equivalently, strongly sequentially 
continuous; equivalently, strongestly continuous) state of %Ao~, the weak closure of Mo, 
where ts the canonical isomorphism of A with its representing function system. We 
call the sets in N, “the permanent null sets of ¢.” 

In view of results of Dixmier, Dye, and Griffin, we could employ, in place of the 
strongly sequentially continuous states f of %~, those states g of Y%)~ arising from 
unit vectors x in 3 (g(A) = (Ag, x)). The strongly sequentially continuous states 
are those strongly continuous on the unit sphere of %->. A projection EF (and its 
range) are said to be cyclic under a C*-algebra Yt when there exists a vector x such 
that [Yr], the closure of the manifold consisting of the vectors Az, with A in %, 
is the range of FE. If this is the case, then F lies in Wl’, the set of operators which 
commute with U. We call this set “the commutant of Yl’ and employ the notation 
Yt’, throughout, for the commutant of %. We say that a projection £ in a ring of 
operators Yl is o-finite when each orthogonal family of projections in A contained in 
i} has at most a countable number of members (we say that F is o-finite relative to 
Yt when the ring to which we are referring is not clearly indicated by the context). 
The ring Mis o-finite when its unit element is o-finite. 

DEFINITION 2. A ring of operators MX ts said to have decomposability character a 
when a is the least cardinal number such that the unit element in A is the sum of a 
family of mutually orthogonal projections with cardinal a, each projection cyclic under 
4’. The ring A is said to be purely decomposable with character a when the ring UP 
has decomposability character a for each central projection P in A. 

The following lemma provides a decomposition of rings of operators which we shall 
employ later. 

LemMaA 1. Jf % is a ring of operators acting on a Hilbert space of dimension a, then, 
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corresponding to each cardinal number b not exceeding a, there is a central projection 
P, in A such that, for each central projection Q which is o-finite relative to the center of 
YA, either QP, = 0 or AQP, has decomposability character b. If (Q»), 6 < a, is a set of 
central projections in A such that Q, has the decomposability property of Py and & » < aQ» 
= ], then Q, = P», for eachb < a. 

We refer to the projection P, constructed in the above lemma as ‘‘the central 
portion of % of character 6.”’ 

DerFinition 3. Jf YU ts a C*-algebra and ¢ ts a representation of A as an algebra 
% of operators on the Hilbert space H, then, fer each unit vector x in #, we call the 
ideal, Nz, of Borel sets, S, in the pure-state space, X, of UA such that S ts a null set of 
the Borel measure induced on X by each positive extension of fed, where d ts the canoni- 
cal isomorphism of A with its representing function system and f is the state of Mo de- 
fined by f(A) = (Az, x), for each A in Mo, a null ideal of g. If the projection in X’ 
with range [Y%x]| is a maximal cyclic projection in A’, we call N, “a characteristic null 
ideal of ¢.”’ 

In addition to the decomposition of the lemma above, we will have use for a part 
of the usual type decomposition of rings of operators. We recall the fact that a 
ring of operators Y contains a central projection P such that MQ is infinite for each 
nonzero, central projection Q contained in P, and % (J — P) is finite (i.e., Q is an 
infinite projection in % and J — P is finite). A projection P in A, with these proper- 
ties, is unique. We call P “the infinite central portion of AX” and/ — P “‘the finite 
central portion.” 

DEFINITION 4. With the notation of Definition 3, if P is the infinite central por- 
tion of %o~, Q the finite central portion of A’, and we let R be PQ, then we call the 
family of characteristic null ideals of the representation ¢, of MX defined by ¢\(A) = 
¢(A)R “the ideal class of ¢.”’ 

The following definitions are intended for that portion of the theory which in- 
volves finite: rings with finite commutants. We could have accomplished their 
effect by employing Griffin’s coupling operator appropriately. This reference, 
however, involves the use of Dixmier’s generalization of the Murray-von Neumann 
trace, a device unnatural to the present theory, which requires, in this connection, 
nothing more than the elementary dimension theory of projections in a finite ring 
of operators. 

Derinition 5. If E is a projection in a finite ring of operators, A, of Type 11, or 
Tm, we call the infimum of all ratios, a/b, such that E is contained in a orthogonal, 
equivalent copies of a projection which has b orthogonal, equivalent copies, “the upper 
dimension of E.”” We denote this number by “u(E)” and define u(@) to be 0. The 
infimum, (Ee), of the numbers u(PE) as P ranges over the nonzero central projections 
in Wt is called “the lower dimension of E.”’ For an arbitrary finite ring A and pro- 
jection E in A, we take as u(E) the supremum of the numbers u(PE) relative to AP, 
where P is the maximal central projection in A such that XP is of Type I], or Im, for 
some m, with (2) defined formally as above. 

The above definition leads to a convenient, trace-free dimension theory for pro- 
jections. 

Derinition 6. Jf Y ts a finite ring of operators with a finite commutant A’ and 
o-finite center, we define “the upper coupling number of A, A’ ” to be 1/1"), where kb’ 
is a maximal cyclic projection in A’, provided that k' # I, and to be u(E), where E is a 
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maximal cyclic projection in UA, if FE’ = I. In case A does not have a a-finite center, 
we define the upper coupling number of A, A’ to be the supremum of the upper coupling 
numbers of AP, XP as P ranges over the nonzero, a-finite, central projections in YA. 
We denote the upper coupling number of A, A’ by “(A, AY.” | 

LemMA 2. Jf WU ts a finite ring of operators with finite commutant WX’, then, corre- 
sponding to each positive real number r, there is a central projection P, such that 
(MQ, X’Q) <r, for each nonzero, central projection Q contained in P,. In addition, 
P, is equal to the unit element in A for r not less than (A, A’). 

We call P, of the above lemma “‘the central portion of 4, A’ with upper coupling 

number r.” , 
3. Principal Results-—The following extension theorem is basic to the theory. 
The detailed account of this theory will contain a stronger result both with respect 
to the generality of maps considered and with regard to the systems being mapped. 
The stronger result will be more in the spirit of a measure-theoretic proposition 
involving a noncommutative measure space with a noncommutative range for the 
integral; for the present, however, the simpler statement will suffice. 

THEOREM |. Jf ¢ ts a representation of the C*-algebra %, acting on the Hilbert 
space %, inlo the C*-algebra Mz acting on Ke, then ¢ has an extension (strongestly 
continuous) mapping %,~ into A.~ if and only if Ny & N,, where Xd is the canonical 
isomorphism of %, with its representing function system. If Ny, = Ny, then the ex- 
tension of g ts an tsomorphism. 

DEFINITION 7. Let A be a C*-algebra, ¢ a representation of A as an algebra Ay 
of operators on the Hilbert space 3, P the infinite central portion of Mo’, Q the infinite 
central portion of %A-(I — P), and R =I — P — Q. Let P,' be the central portion 
of Mo’ P of character a, for each cardinal a not exceeding the dimension of KR. For each 
positive real number or infinite cardinal a, define P, = Q + R, + DY Ps’, where 

bsa@ 


R, ts the central portion of Mo~R, %o'R with upper coupling number a (we use the ob- 
vious convention that R, = R with a an infinile cardinal). We associate with ¢ a 
“multiplicity function,” f,, which assigns to 0 the ideal class of ¢ and to each positive 
real number or infinite cardinal a the ideal of permanent null sets of the representation 
Ga Of A defined by ga(A) = ¢(A)Pa. 

Concerning this definition, we remark that it is appropriate to use the number 0 
for the Q portion of the representation, for this is the portion which is infinite with 
finite commutant and, so, has “coupling” 0. 

THEOREM 2. T'wo representations ¢1, ¢2 of a C*-algebra A as algebras of operators 
%,, U%, acting on Hilbert spaces KH, KH. are unitarily equivalent (1.¢., there exists a uni- 
tary transformation U of %, onto KR, such that Ug(A)U~' = (A), for each A in M) 
if and only tf their associated multiplicity functions are identical. 


* The author is a Fulbright grantee. 









ABELIAN VARIETIES OVER FINITE FIELDS 
By Serce LANG 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF CHICAGO 


Communicated by S. Mac Lane, December 17, 1954 


We shall generalize to Abelian varieties the well-known fact that an elliptic 
curve over a finite field always has a rational point (see Theorem 3). Our first 
theorem is purely algebraic and generalizes parts of Chow’s Theorem 3.' Its proof 
uses techniques of Matsusaka,? to whom I am indebted for valuable discussions 
concerning this matter. Given a separable extension k,/k and a function field 
K,/ky, we give a sufficient condition for the existence of a function field Ak which 
lifts to Ky over ky. 

THEOREM 1. Let k be a field, k,/k a finite separable extension, and k*/k any ex- 
tension linearly disjoint from k, over k. Let Ky/k, and K*/k* be two function fields 
such that K, ¢ K*k;. (We assume, of course, that K, ts free from k* and K* is free 
from k, over k.) If KyKy’ = Kiki’ for every isomorphism o of K*k, over K*, then 
there exists a unique function field K/k such that K ¢ Ky, K ¢ K*, and Kk, = Kky. 

A geometric formulation of Theorem | will be given below. For applications to 
finite constant fields, we shall not need Theorem 1. However, the construction 
used in its proof will be needed to prove the following essential result. 

THeoreM 2. Let K/k be a function field, ki/k a finite separable extension. If 
Kk, is an Abelian function field over k, then there exists a model V of K over k which 
becomes biregularly equivalent to an Abelian variety over ky. 

We now prove Theorem 1, beginning with the uniqueness. If A, K’ are two 
function fields over k satisfying the required conditions, then KA’ ¢ Ky, AK’ ¢ 
K*; hence KK’ is regular over k. If ky = k(a@), then a has the same degree over 
KK’ as over k, and hence KK’ = K. 

To prove the existence of K, we shall use the geometric language. Our hypoth- 
eses can then be stated as follows. We can write K,; = ky(x;) and K* = k*(x*), 
where 2; is the generic point of a variety V;/k; and x* is the generic point of a va- 
riety V*/k*. There is a rational map 7;: V* — V, defined over k,k* such that 
T\(x*) = 2. Without loss of generality, we can assume that V; and V* are in the 
same projective space, that k; is the smallest field of definition for V, containing &, 
and k,k* the smallest for 7, containing &*. Then the condition on the isomorphism 
can be expressed as follows: If k,, V;, 7; (@ = 1,...,m) are the conjugates of 
ki, Vi, 71, respectively, over k, k, k*, there exist birational correspondences R,: 
V,— V, defined over kik, such that commutativity holds in the following diagram: 


\* 
en oh: Neda 
Py \r 


SE pereiee. 2 
R, 
If (x,) are the conjugates of (x;) over k*(x*), then the 0-cycle 5>(z,) is rational over 
k*(x*). Let & be the Chow point of this cycle. The following statements will 
prove that — has a locus W over k and that K = k(é) satisfies our requirements. 


174 
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1. &k(E) © k*(x*). This follows immediately from what has been said above. 

2. ky(ay) © k(t). The eyele > (x,) is rational over k,(), and (z)) is separable 
over /)(&). Let y be a conjugate of (2;) over 4,(&), and let o be an isomorphism of 
the universal domain mapping x2; —~ y over k,(é). Then o@ leaves > (20 invariant 
because the cycle is rational over k,(g). Hence y = 2x; for some j. But the fact 
that o leaves k, fixed implies V," = V;:; hence x, lies in Vy, and that is the case 
only if 7 = 1 (because the conjugates of V; are distinct and so have distinct generic 


points). Having proved that x; is separable over k,(£) and that it is fixed under 


al] isomorphisms leaving /;(£) fixed, we conclude that it is rational over k,(&). 

3. ky(E) © ky(ay). We must prove that >>(2,) is rational over k,(2;), and for 
this it suffices to prove that }°R, is rational over hk; because S°(x,) = YO(R,(2)). 
Hence it suffices to prove that if ¢ is an automorphism of the smallest normal ex- 
tension k’ of k containing k, which leaves k; invariant, then ¢ leaves }°R, invari- 
ant. But o can be extended to an automorphism o* of k’k* leaving kik* fixed, and 
o* leaves 527; invariant because }°7’, is rational over k*. Recalling the commu- 
tativity of the above diagram, we see that ¢ leaves > R, invariant, as desired. 

!. k(t)/k is a regular extension. This follows immediately from statements 
!and3. This concludes the proof of Theorem 1. 

We can now prove Theorem 2. Starting with k* = k in Theorem 1 and V* a 
model of K* = K, we have constructed another model W of K and a birational map 
F:V,— W such that F(27;) = &. It is clear that F~-' has no fundamental points 
on W. Taking the model V; of A, = Ak, to be an Abelian variety, we shall use an 
argument due to Matsusaka? to prove that F is everywhere defined on V;. Let 
k’ be the smallest normal extension of k containing ky. We note that the product 
rm XX... X& tm has a loeus over k’. By Weil’s Theorem 6° and the fact that k’(2,) = 
k’ rn), We conclude that our locus is itself an Abelian variety, biregularly 
equivalent to V, over k’. Hence the map F is everywhere defined on Vi. If we 
now take V to be the normalization of W over /, then the induced correspondence 
between V, and V (over /;) is biregular. This proves Theorem 2. 

Corotiary. Let Kk be a function field. Let k, be the separable closure of k. 
If Kk, isan Abelian function field over k,, then there exists a model V of K over k having 
the following property: If ky/k is finite separable, and Kk, is an Abelian function field 
over ky, then V isan Abelian variety over ky. 

This follows immediately from Theorem 2 and the following lemma. 

Lemma. Let Vk be a variety which becomes biregularly equivalent to an Abelian 
variety over the separable closure k, of k. Then V is an Abelian variety over k if and 
only if V hasa rational point ink. 

Proof: Vf V is an Abelian variety over k, then the unit element is a rational 
point. Conversely, let 2 be a rational point. Let f: Vox V — V be a law of 
composition defined over /, for which x is the zero element. For any automorphism 
o of kh, /k, f? is a law of composition on V and x is again a zero element for f’. Since 
two laws of composition differ by a translation, we must have f’ = f, i.e., fis defined 
over ky. 

From new on, we assume that / is a finite field with g elements. 

"THEOREM 38. Let K /‘k be a function fie ld over the finite fie ldk. If Kkisan Abelian 
function field over k, then K is an Abelian function field over k. In particular, K has 


a rational place. 
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Geometrically, the theorem states that if a variety Uk becomes birationally 
equivalent to an Abelian variety over the algebraic closure of k, then it is bira- 
tionally equivalent to an Abelian variety over k. (If k is a field with a discrete 
valuation, and x is a prime element, then the curve defined by the equation X" + 
rY" + 2°Z" = 0, with n = 3, shows that the conclusion of the theorem does not 
remain valid. In fact, if the curve has a rational point in a separable extension 
k, of k, then n divides [h,:k] because two terms of the equation must have the same 
absolute value in /y.) 

Proof of Theorem 3:4 Let V be the model of A given by Theorem 2. According 
to the lemma, it suffices to prove that V has a rational point. Let 2 be a generic 
point of Vk, and let x* be obtained from x by raising all co-ordinates of 2 to the 
qth power. Then zr‘ is also a generic point of Vk, and (2 X x") has a locus over / 
which is the graph of a rational map ¢: V— V (onto). But over k, V isan Abel- 
ian variety, and hence ¢ = ¢o + c, where go is an endomorphism of V and ¢ is a 
constant. Hence the map x — x — ¢ is an endomorphism of VV. The map s > 
x? — x — cis an endomorphism of |’, which is easilv seen to have finite kernel. — It 
is therefore onto V, and hence & — -¢ —eforsome eV. Hence we have 
é? = £, and £is therefore a rational point. 

CoroLuary. Jf a variely Uk becomes birationally equivalent to an Abelian 
variely over any extension of k (algebrate or not), then U is birationally equivalent to an 
Abelian variety over k. 

Proof: Standard specialization techniques. 

1'W. L. Chow, “Abelian Varieties over Function Fields,” Trans. Am. Math. Soc. (to appear 

> T. Matsusaka, “Some Theorems on Abelian Varieties,” Natural Set. Rept. Ochanomizu U'ni- 
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ON MODULI IN CONFORMAL MAPPING 
By H. EF. Ravcu 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF PENNSYLVANIA 
Communicated by Einar Hille, December 10, 1954 


1. Jntroduction._-In his fundamental papers! on a theory of functions of a com- 
plex variable based on potential theory and conformal mapping, Riemann observed 
that two Riemann surfaces associated with algebraic functions and of the same 
genus, p, cannot necessarily be mapped on one another conformally but that the 
totality of conformally equivalent classes of such surfaces depends in an unspecified 


way on 3p — 3 continuous complex parameters, called moduli. In a note* in these 
PROCEEDINGS (to which I refer the reader for a more detailed discussion of the pre- 
ceding) I established a theorem which definitely exhibits any one of certain subsets 
of the p(p + 1)/2 nontrivial periods of the normal integrals of the first kind as a set 
of numerical moduli, i.e., parameters defined intrinsically for each surface whose 


equality for two surfaces implies their conformal equivalence. 
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Now, following Riemann, Klein,* Schwarz, Schottky,’ and others took up the 
theory of conformal mapping from the new point of view but also with a return to 
its original genesis in mathematical physics, that is to say, the mapping of plane 
domains and, more generaliy, finite surfaces with boundaries. Here, too, the phe- 
nomenon of moduli enters, and, as a counterpart to the count 3p — 3 described in 
the note referred to above,’ one finds that, if one maps, say, a plane domain bounded 
by nm nonintersecting Jordan curves on a canonical slit (say) model, one may pre- 
scribe the left-hand end-point and length of one slit at will, while the end-points 
and lengths of the next n — 2 slits make up 3(n — 2) = 3n — 6 essential real param- 
eters— again a vague and unsatisfactory specification. 

Thus one is faced again with the problem of defining intrinsic numerical moduli 
for finite Riemann surfaces with boundary, and in this note I shall demonstrate how, 


in essences the solution can be reduced to the solution of the problem already de- 


scribed? fér the closed Riemann surfaces, which case is thus the experimentum crucis. 

For didactic purposes I treat the case of a plane domain, D, bounded by the non- 
intersecting Jordan® contours, Cy, k = n, although the proof applies equally 
well to a surface of genus greater than zero. 

2. The Double of D and Differentials——Schottky® had the ingenious and fertile 
idea of associating with D a closed Riemann surface, S, called the double of D, of 
genus p = n — 1, obtained by taking an identical copy, D, of D, identifying corre- 
sponding points on the boundaries of D and D, and assigning the co-ordinate 2 to 
the point in D whose mate in D has the co-ordinate z. The resulting Riemann sur- 
face is one of those called symmetric by Klein.’ Schottky then admits on D only 
those functions f and those differentials which are real on the boundary, C, of D. 
That is because such functions and differentials, after suitable choice of parameter 
near points of (, may be extended to D and hence to S by Schwarz’s reflection prin- 
ciple. 

Thus armed, I can now indicate all the modifications of the previous note? that 
are necessary. First of all, the starting point here as well as there is the important 
observation of Teichmueller,* which serves as a guide but is not used: 

A. The number of quadratic differentials on D which are linearly independent with 
respect to real coefficients is equal to the number of moduli, viz., 3n — 6. 

By a quadratic differential, df? = f dz*, is meant here an everywhere finite, 
analytic quadratic differential on D, i.e., one where f is analytic on D without ex- 
ception and f dz? is realon C. If one can show that the number of real-independent 
dé? on D equals the number of complex-independent df? on its double, S, then A 
will follow from A in the previous note,’ where it is asserted that the latter number 
is3p — 3 = 3(n — 1) — 3 = 3n — 6. In the first place, any member d¢,? of a real 
basis on D can be extended to a dé,? on S. Now any linear combination >>,c, 
dé,? with complex ¢, is not a dé? on D, not being real on C. If now dé? = Ye, dt? = 
0 on S, then the same would be true of df°* = d¢*(z), z being the co-ordinate of a 
point on S and 2 being the co-ordinate of its symmetric image. But d¢,?* = 
de,*() = d¢,?. Therefore, dg? + d¢?* = >°,(Re e,) dé,? would be zero on D. That 
being impossible, the d¢,? form a complex basis on S. 

I recall now the fact stated in B of the previous note? that, if df,, » = 1,..., p 
are the differentials of the normal integrals of first kind on S, then a basis for the 
d¢? of S consists of the products df, dt;, where i runs through a certain set of indices, 
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J, and j runs through a similar set J—-abbreviated to (7, j) « U/, J), all provided 
that S ts nol hyperelliptic when p = n — 1 > 38. 

Lemma l. Whenn — 1 > 3, the double, S, of the n-ply connected plane domain D 
is nol hyperelliptic in the general case.* 

Proof: Ti Cy i = 1,...,n" — 1, are the interior boundary curves of D, then, 
properly oriented, they become the retrosections, 6,, making up half of a canonical 
system of retrosections, the other half, y,, ¢ = n— 1 = p, being constructed 
out of identical cuts, [, and [,, joining C,; to C)+; and C, to C)+1, in D and D, respec- 
tively, the cuts being run through in the order C; > Cj4, = Cj+1 > C;. If S is 
hyperelliptic, it admits a conformal involution, 7 (the exchange of sheets in a two- 
sheeted representation), which leaves one set of retrosections pointwise invariant, 
while each member of the other, conjugate, set admits 7. Letting the y, be the 
former and the 6, the latter, one sees that D must also admit 7. But that is mani- 
festly impossible in the general case, for n — | > 3, as one sees by mapping D on a 
proper canonical domain,’ where 7 becomes a Euclidean motion of D into itself. 
Only the ring, n — 1 = 1, and the double ring, n — 2 = 2, admit a rigid involution 
without exception. 

The last consideration preliminary to formulating the main theorem is to show 
how the normal integrals of the first kind on S are obtained from the harmonic 
measures on D,. The ith harmonic measure on D is the harmonic function w, which 
assumes the values zero on C, and 6;, 7,7 = 1,...,n — 1, on (Cy. Extending 


w, toan analytic function, w, = w, + ty, one finds that w, adds a purely imagi- 


nary quantity, 7P?,,, every time the argument describes C;. The real part of w, 
being constant on (, (Ow,/O0z) dz + (Ow, Oz) dz = O there, i.e., 1(Ow, Oz) dz is real 
on C and can be extended to S. If y, and 6, are as described in the proof of Lemma 


1, then clearly 


Ow, * Ow, 
[ dz = 23... | gm 9? (1) 
. 7 Oz * é Oz 


If, therefore, I set 7f, = '/2(Qw,/0z) dz, then S° dé, will be the ith normal integral 
of the first kind on S whose periods #,, = !/9P,,. I designate the P,, as the periods 
of the harmonic measures. 

3. The Main Theorem.-lf the df, dé;, (7,7) « Ud, J), constitute a basis (complex- 
independent) for the df? on WS, then, according to the observations in the proof of A, 
the products (Ow,/0z) (Ow,/Oz) dz*, (¢, J) « U/, J), will form a real-independent 
basis for the df? of D. Now the main theorem is 

THrokem 1. Let D be a plane domain bounded by nonintersecting Jordan® contours, 

., n, whose double is not hyper-elliptic, having harmonic measures w,, 
n — 1, and associated analytic functions w, = w; + tu; Let (Ow, Oz) 
(Ow, /Oz) dz*, (1, 7) e UI, J), form a real-independent basis for the quadratic differen- 
tials on D.° Let Py, be the period of the ith harmonic measure over C,, j 
n— 1. Ifa similar domain D’ with functions w,' and periods P,, is such that 


Py=P,', (je), 


then D’ can be mapped 1-1 conformally on D. 
In view of formula (1), Theorem 1 is almost a corollary of Theorem | of the pre- 
vious note.” In facet, the latter implies under the hypotheses of the former that 
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the double, S’, of D’ is conformally equivalent to S. However, the resulting con- 
forma! map may not take D into D’. Therefore, | will start afresh but only supply 
those arguments which, while parallel to those of the earlier note, are different. 

The starting point, as before, is to set up a differentiable 1-1 map, ¢, between D + 
Cand D’ + CC’. One then minimizes 


J(e) = "2 SS, (E + @ dx dy (3) 


among all ¢ deformable into ¢,'' where EF and G are obtained from 
du d® = KE dx? + 2F dr dy + G dy’, 


where, in turn, D’ is in the w-plane and w g(x, y), D being in the z 
plane. 
If go is the minimizing map, on rewriting equation (4) as 


dw dt = b dz dz + 2 Re (a dz”), 


where b = (F + G)/2 anda = (EF — G)/2 — iF, one sees as before that a dz?, a not 
necessarily analytic quadratic differential, is analytic when ¢ = ¢o (another proof. 
of this below) and that go is conformal if a dz* = 0. 

I shall demonstrate in Section 4 that the analyticity of a dz* is a consequence of in- 
terior variation of go only, while, and this is the crux of the paper, the variation of the 
boundary correspondence allows the deduction that a dz* is real on C. 

If one grants this for a moment, then the proof of Theorem | is immediately re- 
duced to that of Theorem | of the preceding note;? for all the ¢, including go, can 
be extended to maps of S on S’, while a d2* can now be extended to S, where, thanks 
to formulas (1) and (2) the proof in the earlier note shows that adz?/b dz * must 
be orthogonal to a basis df, d&,, (7, 7) « U, J), on S and, hence, vanish identically. 

4. The First Variation of the Douglas-Dirichlet Functional J(¢).—-To establish 
the preceding italicized assertion in Section 3, I observe that J (go) may be written 


as 


where dz-dzZ is the exterior product, and calculate the first variation where go 1s 
altered by a small amount. 

To alter gp is to map D into itself and follow that by the mapping ¢» into D’. 
Now if ¢ is sufficiently small and w(z, 2) is a continuously differentiable funetion 
which vanishes outside a neighborhood l of z interior to D, then 

zeU 
2-¢D-U 


defines an interior variation gy + d¢9 of go. 
Now 


dw d@® 


where w o(2’, 3) o(z, 2) 4 ; : ~ (°, while 


te} Re (aw;s)}| + Ole? 
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when ze. | observe that d}b(@ dz — w dz)! is the middle term in brackets multi- 
plied by dz-dz (d is the exterior differential). 
Thus 


4 a i a 
aJ = Qi {| Re (aw,) dz-d2 4 : I d}b(w d? — wdz)} + Of). 
2S St 


Since dz-d2 is purely imaginary, the first term on the right may be written 


2) e Im ffi. aw, dz-dé. 


Now, if 0° denotes the boundary of (, one may transform by Stokes’ theorem to 


; 
get 


6] 2+ «Im SS a,w dz-dt 27 e« Im Sox aw dz + 
te Jou O(w d=? — wdz) + Ole). 


In particular, since J (go) is a minimum while w vanishes on Ol” but is otherwise 
arbitrary, one sees that a, G6. i.e, a is analytic. But if Ol” includes a piece, 
Cy, of C, then w = 0 on Ol”, with the exception of Cp, where w \ dz ds for real 
\, s being the are length on C. Hence, keeping in mind that a, 0, the extremal 
condition becomes 


‘ j /dz\*) ? dzdz dz dz\ 
?) a Im a ds + 1 br — is (). 
Je, V\ds/] § Jey dsds ds ds 


Since the last term is zero and J is an arbitrary real function, one concludes that 
Im (a dz*) 0 on C, which completes the proof of the italicized statement in 
Section 3. 


! Gesammelle mathematische Werke 
>On the Transec odental Moduli of Algebraic Riemann Surfaces,’’ 40, 42, 1955. 
> “Riemann’sche Flachen,’’ Autographierte Vorlesungen (Gottingen, 1891-1892 
Gesammelle mathematische Abhandlungen 
“Uber die conforme Abbildung mehrfach—zusammenhingender ebener Flichen,”’ J. f. Wath., 
83, 300, 1877 
*] actually assume the curves to be continuously differentiable, but the restriction can be 
lifted 
On Riemann’s Theory of Algebraic Functions (Cambridge, England, 1893). 
* Abhandl. Preuss. Akad. Wiss., Jahrg. 1939, No. 22, 1940. 
’ The hyper-elliptic exceptions will be dealt with in a forthcoming note in these PROCEEDINGS, 
where one will also find a description of the canonical domain referred to below. 
” The existence of such J and J is assured by what precedes 
‘Tin = 1. three fixed points must go into three fixed points; if nm = 2, one fixed point into one 
fixed point 





POSSIBLE SEPARATION OF INTERTWINED NUCLEIC ACID CHAINS 
BY TRANSFER-TW UIST 


By Joun R. Plarr 
PHYSICS DEPARTMENT, UNIVERSITY OF CHICAGO 


Communicated h / R S Vudlliken December 2 / aw 


To separate the strands of a twisted rope it is not necessary to unravel them from 


the end. It is mechanically simpler and energetically easier to make a “transfer- 
twist,”’ pulling the strands apart in the middle and letting each strand twist about 
itself. This suggests a method of separating two twisted intertwined helical 
molecules without expending excessive energy. It might be an alternative to the 
method suggested by Delbriick.' 

In particular, Figure 1 shows how transfer-twist might apply to intertwined 

duplicate nucleic acid chains like those in the Watson and Crick model of DNA? 
Fhe 1, and A. halves of the two DNA chains are identical and are the complements 
to the other halves, A,’ and Ay’. Thus A, will intertwine perfectly with either 
A,’ or Ae’ ina“ lock-and-key”’ re- 
lationship but not with A», and 
soon. If we now start with Ay 
linked with As’ and A,’ linked 
with A. and make a transfer-twist 
Which pulls the chains apart at 
their centers of inversion, /; and 
/., then an extensive motion in 
the direction of the arrows will 
cost almost nothing in energy or 
entropy. The energy of twist is 
transferred to the self-twisted 
strands, Ay-Ay’ and A.-A»’. The 
bonds formed at the transfer 
point, 7, are identical with those 
broken, the two processes perhaps 
even facilitating each other to give a very low activation energy. 
” The motive force to drive the twisting forward along the arrows in Figure | might 
come from a change of environment. Perhaps hydration or dehydration of the 
chains might start to be favored, leading to more hydrogen bonding, or less, per 
nucleotide. With hydration favored, water molecules might be added to each 
nucleotide as it separates from its usual bonds in passing the transfer point 7' in the 
forward direction. Once started, there would be a continual energy gain in going 
forward in the same direction and none in going backward. Other environmental 
alterations might favor similar changes in the chain structure and energy per unit 
length. 

After separation of A,-A,’, the regeneration of another duplicate from it can pro- 
ceed in from the ends in reverse fashion, with the self-twisted center untwisting as 
fast as suitable new fragments, a; and a3’, grow and retwist around the ends. 

The general features of transfer-twisting may be followed by manipulating a 


18] 
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two-strand piece of rope. A more elaborate model, constructed as follows, will 


show the lock-and-key relationship in duplicate chains. Take two equal lengths of 
round soft rubber rod (gaskets). Twist each tightly the same number of turns 
and fasten them flat side by side on a table. With a wet knife, slice both together 
from one end almost to the other, keeping the knife parallel to the table and roughly 
near the center line of the rods, but making arbitrary hills and valleys. The top 
halves of the two rods (A; and A.) are then alike and are complementary to the 
bottom halves (A,’ and A.’). The uncut ends are the centers of inversion (/; and 
/,). When the rods are released, each will have its complementary halves tightly 
twisted together; but, when the cut ends are separated, each can be made to twist 
in a unique way with its complementary duplicate in the other rod, as in Figure 1. 

It is possible that two intertwined complementary molecular chains may each 
be made of many segments like those in Figure 1. Opposite segments are dupli- 
cates but are oriented oppositely with their neighbors. Each segment may undergo 
transfer-twist about its own center, as suggested in Figure 2, provided viscous re- 

sistance of the medium does 

not prevent the necessary twist 

of the transverse axes around 

the longitudinal axis. There 

would be no such twist of the 

axes if the segments were equal 

and were alternately right- 

handed and left-handed helices. 

It is tempting to speculate 

on whether such a polysegment 

DNA chain might be a princi- 

pal fragment of a chromosome. 

It has the property that regen- 

Fic. 2 eration could begin simultane- 

ously at the numerous ends of 

segments, /. These would also be the last points to separate from a duplicate and 

the easiest points to break. Hundred-segment chains, when almost ready to sepa- 

rate, would appear of the order of a hundred times shorter and a hundred times 
wider than the extended chains and might be microscopically visible. 

A rope of three or more strands may be converted by transfer-twist into two or 
more self-twisted strands or groups of strands, but these are not in general identical 
in pitch or energy with the original rope. Likewise, a one-strand helix may be 
pulled out at one or more points into self-twisted side arms, but these will all be 
double helices. While transfer-twist could be important in the transformations of 
one-strand molecular helices, the special facility of two-strand transformations 
suggests that, for all suspected helices, special attention should be paid to possible 
two-strand interpretations of the evidence. 

Besides its possible role in self-duplication and separation, transfer-twist might 
serve generally to shrink long chains and to pull distant units together and also to 
bring together two reactants if they were associated, say, with opposite ends of the 
same segment in Figures | or 2. It might facilitate creation of duplicates or com- 
plements even of nonhelical chain molecules if they could be associated with a trans- 
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fer-twist helix in a noninterfering way. And it might be a reversible method, 
largely freed from the randomness of thermal meanderings, by which one molecule, 
reversing the arrows in Figures | and 2, may “examine” its alien, neighbors for 
identity or complementarity so that specificity is insured in any reaction with them. 


I am indebted to Dr. Aaron Novick for much helpful criticism and discussion of 
these questions. 

' M. Delbriick, these ProceEDINGs, 40, 783, 1954. 

2 J. D. Watson and F. H. C. Crick, Nature, 171, 737, 964, 1953; F. H. C. Crick, Sei. American, 
191, No. 4, 54, 1954. That each chain is self-complementary like our A;-A,’ was suggested by 
Crick’s remark that “the figure looks exactly the same whichever end is turned up.” With or 
without this feature, a Watson-Crick “zipper” replication is copied by the present model, if one 
of their chains is identified with a fragment such as A, of Figure 1, points /; and /; being broken. 
But the separation process in Figure | and the polysegment replication and separation processes 
indicated in Figure 2 have no counterpart in their model. 


ERRATA: HAMARTIEXERESIS AS APPLIED TO TABLES INVOLVING 
LOGARITHMS 


In the article of the foregoing title appearing in these ProcEEDINGs, 40, 728-731, 
1954, the following corrections should be made: 


P. 729, line 24: Insert subscript 1 below symbol 7m at right. 
P. 729, line 25: For 117724.51130 read 11 77245.51130. 
P. 730, top line: A closing bracket, ], should precede the = sign. 
P. 730, second line below table: For = read = or =. 
P. 731, line 14from bottom: For Duarte read Alliaume. 
P. 731, line 12 from bottom: For 26823155 read 26823 155. 
P. 731, bottom line: For 26090359 read 26090 359. 
Horace S. UHLER 
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